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SUMMARY 


An  extensive  literature  review  on  the  effect  of  interface/ 
interphase  on  composite  properties  was  followed  by  experiments 
which  illustrate  the  significance  of  interphase  in  composite 
performance . 

The  review  covered  adhesion  theory,  wetting  characteristics 
of  carbon  fiber,  load  transfer  mechanisms  in  a  fibrous  composite, 
test  methods  to  measure  and  evaluate  fiber-matrix  interfacial 
bond  strengths,  effect  of  moisture  at  the  interface  and  properties 
of  the  interface/interphase,  including  impact  toughness,  on  com¬ 
posite  performance.  A  critical  discussion  of  the  literature  is 
presented  along  with  appropriate  citations  and  references. 

In  contrast  to  adhesion  in  glass  fiber /epoxy  composites, 
carbon  fiber/epoxy  adhesion  is  not  simply  promoted  by  "coupling 
agents"  but  rather  it  results  from  physical  and/or  chemical  inter¬ 
actions  with  the  matrix.  Much  work  has  been  done  to  characterize 
the  surface  energetics  of  carbon  fiber  but  no  correlation  has 
been  identified  between  wetting  behavior  and  in ter facial  bond 
strength.  Many  techniques  have  been  developed  to  study  fiber/ 
matrix  interfacial  bond  strength  but  few  are  applicable  to  car¬ 
bon  fiber.  Preliminary  data  have  shown  that  single  filament 
critical  length  determination  is  applicable  to  evaluate  carbon 
fiber  -  epoxy  bond  strength.  This  technique  was  adapted  in  this 
study  to  evaluate  finish  variants  on  carbon  fiber  in  epoxy  com¬ 
posites  . 

Moisture  at  the  interface  of  carbon  fiber/epoxy  composites 
has  been  shown  to  have  little  effect  on  composite  properties . 

The  decrease  in  composite  properties  in  humid  environments  is 
primarily  due  to  plasticization  of  the  matrix  by  water. 

Improvement  of  impact  toughness  of  carbon  fiber /epoxy  com¬ 
posites  can  be  achieved  by  modification  at  the  fiber/matrix  inter¬ 
face.  However,  this  is  usually  associated  with  a  significant 
decrease  in  mechanical  properties.  This  is  because  weak  inter¬ 
facial  bonds  allow  more  impact  energy  to  be  absorbed,  thereby 
increasing  the  impact  toughness.  A  few  workers  have  suggested 
that  "tailoring"  interphase  properties  can  improve  toughness 
without  lowering  the  mechanical  properties,  but,  this  concept 
has  not  yet  been  demonstrated  experimentally. 

In  this  study,  experimental  results  were  obtained  which  show 
the  effect  of  interphase  on  composite  performance  and  illustrate 
its  significance.  Various  finish  variants  were  formulated,  based 
on  different  chemical  and  mechanical  properties,  and  applied  to 
Celion®  6000  carbon  fiber.  Wetting  behavior,  studied  by  the 
Wilhelmy  wetting  force  technique,  and  interfacial  bond  strength, 
evaluated  by  single  filament  critical  length  determination  were 
characterized  for  these  fibers.  Composite  properties  at  room  and 
elevated  temperatures  and  impact  toughness  were  also  measured. 
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It  is  concluded  that  good  wetting  is  essential  to  provide 
good  adhesion  but  better  wetting  properties  do  not  necessarily 
translate  to  better  bond  strength.  Also,  it  is  shown  that  with 
a  weak  interfacial  bond  strength,  impact  toughness  is  improved. 
But,  most  significantly,  a  30%  increase  in  impact  toughness  was 
measured  without  loss  of  room  temperature  and  elevated  tempera¬ 
ture  composite  mechanical  properties  by  "tailoring"  the  inter¬ 
phase  properties. 
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INTRODUCTION 


In  the  last  decade,  much  attention  has  been  focused  on  fiber- 
reinforced  polymeric  matrix  composites.  These  composite  materials 
are  both  versatile  and  complex.  Their  versatility  stems  from  a 
wide  choice  of  constituent  materials  available  and  from  the  variety 
of  ways  in  which  composites  can  be  fabricated  to  provide  a  combina¬ 
tion  of  desired  properties.  These  tailor-made  properties  are  fre¬ 
quently  not  available  and  cannot  be  achieved  in  conventional  isotro¬ 
pic  materials.  On  the  other  hand,  they  are  complex  by  virtue  of 
their  chemical  and  mechanical  nature  and  thus  should  not  be  regarded 
as  a  single  material,  but  as  a  material  system.  Inherent  in  this 
multi-phase  material  system  is  the  fiber/matrix  interface  which  is 
the  critical  link  that  provides  the  structural  integrity  for  the 
fibrous  composites. 

Conventionally,  the  interface  is  regarded  as  the  bond  between 
fiber  and  matrix.  This  interfacial  bond  is  considered  to  have  zero 
thickness  and  to  result  from  the  interaction  between  the  fiber  surface 
and  the  matrix  material.  The  interphase  region  of  a  fibrous  composite 
is  the  area  immediately  adjacent  to  and  including  the  interface  and 
extends  a  finite  distance  into  the  bulk  matrix  material.  This  concept 
is  illustrated  in  Figure  1. 

An  extensive  review  of  the  literature  indicates  that  much 
work  has  been  done  on  fiber/matrix  interface  characterization. 

These  include  the  evaluation  of  fiber  surface  energetics  via  wet¬ 
ting  force  measurement;  surface  elemental  analysis  via  ESCA, 

AUGER,  XPS,  etc.  and  surface  topography /morphology  studies  via 
SEM.  However,  little  effort  has  been  spent  to  correlate  this 
information  with  composite  performance.  In  short,  the  role  of 
the  interface  in  determining  composite  properties  is  not  known. 
Furthermore,  little  attention  has  been  given  to  the  understand¬ 
ing  of  the  nature  of  the  interphase  and  how  it  affects  composite 
properties.  In  particular,  the  load  transfer  mechanism  in  the 
interphase  region  is  not  clear.  This  report  is  intended  to 
clarify  these  areas  of  uncertainty  and  to  illustrate  the  influ¬ 
ence  of  the  interphase  region  on  composite  performance. 
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BACKGROUND  AND  LITERATURE  REVIEW 

ROLE  OF  THE  INTERFACE  IN  COMPOSITE  MATERIALS 

Among  many  factors  that  govern  the  characteristics  of  a 
fibrous  composite,  it  is  clear  that  the  adhesion  at  the  fiber/ 
matrix  interface  plays  an  important  role  in  structural  integrity 
of  these  two-phase  materials.  Since  there  are  a  limited  number 
of  ways  to  directly  study  adhesion  at  the  interface,  composite 
failure  mode  and  shear  strength  are  commonly  employed  to  assess 
interfacial  bond  strengths. 

Two  extreme  cases  generally  occur  in  interfacial  failure 
mode  studies.  For  weak  interfacial  bonding,  fiber  "pull  out" 
is  responsible  for  the  rupture.  This  is  because  the  fibers  are 
not  bonded  by  the  matrix  and  they  slide  in  their  cavities. 
Therefore,  load  cannot  be  transferred  in  such  a  discontinuous 
medium.  In  this  case,  the  composite  fails  in  a  shear  mode  and  the 
shear  strength  is  usually  low.  On  the  other  hand,  if  the  inter¬ 
facial  bonding  is  strong,  then  the  matrix  will  transfer  the  load 
to  the  fibers  until  it  reaches  the  breaking  strength  of  the 
fibers.  In  this  case,  the  failure  is  sudden  and  catastrophic. 

A  higher  shear  strength  is  associated  with  this  kind  of  tensile 
failure . 

Virgin  carbon  fibers  exhibit  lower  bondability  with  polymer 
matrixes  (e.g.,  epoxy,  polyester)  and  surface  treatments  are 
employed  to  increase  the  interfacial  bond  strength.  All  commer¬ 
cially  available  carbon  fibers  are  surface  treated.  These  surface 
treatments  are  generally  oxidative  in  nature  and  provide 
functional  groups  on  the  fiber  surface  which  increase  wettability 
and/or  bondability.  Much  of  the  early  work  has  emphasized 
increases  of  composite  shear  strength. 

Goan  et  al(l)  studied  a  sodium  chlorate/sulfuric  acid 
oxidation  process  which  doubled  composite  short  beam  shear 
strength.  However,  this  process  degraded  the  fiber  tensile 
strengths.  A  4%  loss  in  fiber  strength  (impregnated  strand 
tensile  test)  was  observed.  Druin  et  al(2)  studied  gas  phase 
oxidation  processes  where  similar  results  were  obtained. 

Carbon  fiber  strength  is  surface  flaw  sensitive.  Extensive 
oxidation/etching  can  increase  surface  flaws  which  leads  to  a 
reduction  in  properties.  Ehrburger  et  al, (3)  Molleyre  and 
Bastick, (4)  Fitzer  et  al(5)  and  Duffy<6)  studied  various  types 
of  oxidation  processes  to  optimize  the  incorporation  of  functional 
groups  onto  the  fiber  surface  and  reduce  the  possibility  of  exten¬ 
sive  oxidation/etching. 

Ehrburger  et  al(7,8)  found  that,  by  the  insertion  of  weak 
acidic  groups,  interlaminar  shear  strength  (ILSS)  could  be  en¬ 
hanced  just  as  much  as  with  the  incorporation  of  strong  acidic 
groups  on  AC  (type  II),  ex-acrylic  carbon  fiber  surfaces.  However, 
the  mode  of  failure  was  different.  A  shearing  mode  of  fracture  was 
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observed  with  weak  acidic  groups  while  tensile  failure  was  ob¬ 
served  with  strong  acidic  groups.  A  possible  explanation  suggested 
by  the  author  for  this  observation  was  that  the  weak  acidic  groups 
increased  the  surface  energetics  of  the  fibers  and  resulted  in 
better  wetting  and  good  adsorption  which  enhanced  the  coupling. 

The  strong  acidic  groups  -  mainly  carboxylic  or  phenolic  groups  - 
are  coupled  with  the  matrix  via  tighter  and  stronger  primary 
bonds  at  the  interface  which  create  a  brittle  composite. 

Although  little  is  known  about  the  adhesion  mechanisms 
involved  in  a  fibrous  composite  material,  it  is  clear  that  the 
role  of  the  interface  in  a  composite  material  is  significant. 
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EFFECT  OF  INTERPHASE  PROPERTIES  ON  COMPOSITE  PERFORMANCE 

Commercially  available  carbon  fibers  are  usually  sized  with 
an  epoxy  compatible  finish  to  enhance  handling  characteristics. 
Epoxy-based  systems  are  generally  used  since  epoxy  resins  are  com¬ 
monly  used  as  components  of  the  matrix  materials  and  such  mate¬ 
rials  can  be  easily  applied  to  the  fibers.  These  sizes  also  act 
as  lubricants  as  well  as  "loose"  binders  which  prevent  fiber  dam¬ 
age  and  minimize  fiber  fuzz.  The  loose  binder  should  not  hinder 
the  spreading  of  the  fibers  during  the  prepregging  process,  but 
should  allow  the  resin  to  penetrate  bundles  and  wet  individual 
fibers  to  eliminate  dry  spots  or  voids. 

Sizing  can  be  accomplished  by  passing  tows  of  fibers  through 
a  bath  containing  a  dilute  solution  of  the  finish  material.  Fol¬ 
lowing  such  a  coating,  the  fibers  must  be  dried  to  remove  solvent. 
Better  adhesion  properties  were  observed  by  Drzal  et  al(l)  and 
Goan  et  al(2)with  commercial  epoxy-sized  fibers.  This  phenomenon 
can  be  explained  on  the  basis  that  the  fibers  were  sized  with  a 
low  viscosity  fluid,  thus  enhancing  wetting.  Furthermore,  the 
spontaneous  spreading  of  the  sizing  solution  minimizes  trapped 
air  pockets  at  the  interface  between  fiber  and  finish. 

The  overall  effect  of  size  materials  on  composite  properties 
is  not  fully  understood.  This  is  partly  because  the  size  materials 
are  applied  onto  the  fiber  in  very  thin  layers,  normally  1.5%  by 
weight  of  the  fiber,  which  is  equivalent  to  an  approximately  800  A 
thick  coating  (Appendix  A) . 

Few  analytical  techniques  can  be  applied  to  study  the  inter¬ 
phase  of  a  carbon  fiber  epoxy  composite  because  of  its  small 
dimension.  Furthermore,  the  fiber  geometry,  and  the  low  atomic 
number  elements  on  the  surface  further  limit  the  applicable  analyti¬ 
cal  techniques.  These  add  to  the  difficulties  in  understanding  the 
role  of  the  interphase  in  composite  properties.  However,  some  of 
the  effects  of  interphase  can  be  observed  in  mechanical  tests . 

Arai  et  al(3)  showed  an  increase  in  both  interlaminar  shear 
strength  (ILSS)  and  flexural  strength  of  a  carbon  fiber  epoxy  com¬ 
posite  by  electrolytic  polymerization  of  vinyl  monomers  onto  the 
carbon  fibers.  Pinchin  and  Woodhams(4)  pyrolyzed  graphite  onto  car¬ 
bon  fiber  surfaces  and  increased  the  surface  energetics.  This  was 
reflected  in  an  increase  in  shear  (ILSS)  and  flexural  strengths. 
Furthermore,  a  slight  improvement  on  impact  toughness  was  also 
measured.  An  hypothesis,  based  upon  this  phenomenon,  is  that  the 
pyrolytic  graphite  deposit  forms  a  sheath  around  the  fiber  which 
creates  an  additional  interface.  During  the  debonding  and/or  crack 
propagating  processes,  this  additional  interface  absorbs  more  energy 
by  exposing  additional  free  surfaces. 

Broutman  and  Agarwal^  •  7 )  did  a  theoretical  study  on  composite 
properties  with  an  interphase  layer.  Assuming  both  phases  behave 
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elastically  with  perfect  bonding  at  the  interface,  it  was  shown 
that  composite  strength  reached  a  maximum  when  the  interphase 
modulus  reached  a  value  of  69  MPa  (104  psi) .  It  should  be  noted 
that  Broutman  used  the  term  interface  modulus  but,  according  to 
the  definitions  of  interface  and  interphase  described  earlier, 
interphase  modulus  is  more  appropriate  because  interface  is  re¬ 
garded  as  zero  thickness.  Further  increases  in  the  interphase 
modulus  do  not  change  the  composite  strength.  This  is  shown  in 
Figure  2.  However,  composite  modulus  continues  to  increase  with 
increasing  interphase  modulus  (Figure  3) .  This  further  increase 
in  modulus  without  composite  strength  changing  further,  indicates 
that  composite  elongation  can  be  maximized  at  this  critical  inter- 
phase  modulus  (Figure  4) . 

Although  Broutman' s  theoretical  study  was  not  supported  by 
experimental  data,  the  critical  interphase  modulus  determined  in 
his  study  was  significantly  lower  than  the  modulus  of  the  matrix. 

This  is  similar  to  the  flexible  interlayer  concept  discussed  by 
Benedetto  and  Nicolais  (5) .  This  concept  attributes  the  effective¬ 
ness  of  finishes  to  stress  relief  in  the  interphase  region.  Since 
there  is  a  difference  in  thermal  expansion  coefficient  of  the  fiber 
and  the  matrix,  significant  residual  stress  will  accumulate  at  the 
interface  due  to  elevated  temperature  curing  as  well  as  shrinkage 
of  the  matrix  during  polymerization.  However,  by  incorporation  of 
a  flexible  interlayer,  residual  stresses  can  be  relieved.  Although 
these  two  concepts  are  different  in  approach,  they  both  conclude 
that  the  properties  of  the  composite  can  be  improved  through  selec¬ 
tive  tailoring  of  the  interphase. 

(12) 

Marom  and  Arridge  studied  the  stress  patterns  of  soft  inter¬ 

layers  on  stainless  steel  inclusions.  Reduction  of  stress  concen¬ 
tration  at  the  interface  was  observed  and  an  improvement  in  the  ul¬ 
timate  tensile  strength  was  measured.  However,  Marom  stated  that 
optimization  of  this  soft  interlayer  is  necessary  to  prevent  poor 
transverse  properties.  This  effect  was  reported  by  Gatward  and 
Hull, (13 ) when  they  coated  a  flexibilized  polyvinyl  butyral  polymer 
onto  glass  fiber  and  a  decrease  in  transverse  tensile  properties 
was  measured. 

Taniguchi  et  al (9) incorporated  a  flexibilized  epoxy  at  the  in¬ 
terphase  of  a  carbon  fiber/epoxy  composite  and  preliminary  data 
showed  improvement  in  flexural  and  shear  strength.  Kardoes  et  al(10/H) 
and  Lavengood  et  al(14)  COated  boron  fibers  and  glass  fibers  with  a 
low  cross-linking  density  epoxy  system  to  a  thick,  flexible  inner- 
layer.  Increase  in  mechanical  properties,  both  moisture  conditioned 
and  unconditioned,  were  measured.  A  1000%  increase  in  torsional 
fatigue  life  was  measured  and  Figure  5  shows  the  improvement  of 
transverse  properties  of  glass  fiber/epoxy  composites  by  incorpora¬ 
tion  of  a  soft  interphase. 
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FIGURE  2. 


FIGURE  3. 
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*From  Broutman  and  Agarwal  (Ref.  7) 
Reproduced  by  permission. 
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FIGURE  4 .  COMPOSITE  STRAIN  ENERGY  ABSORBED  AS 
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*From  Broutman  and  Agarwal  (Ref.  7). 
Reproduced  by  permission. 
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*From  Lavengood  and  Micheno  (Ref.  14) . 
Reproduced  by  permission. 


Riess  et  al ^  developed  an  interphase  material  with  poly- 
[ isoprene-b- (styrene-alt-maleic  anhydride)]  block  copolymer: 


CCH2 


1 3 

C  =  CH  -  CH0]  UCH., 
Z  n  Z 


-  CH  -  CH  -  CH] 


m 


C  C. 

/  v  \ 


(Isoprene) 


(SAM) 


According  to  the  compatibility  of  the  polymers,  the  polyiso 
prene  will  agglomerate  to  form  micellar  structures  at  the  inter¬ 
phase  while  the  SAM  block  is  mixed  with  the  matrix  resin  to  form 
a  continuous  phase.  This  inclusion  of  a  second,  elastomeric  phase/ 
is  to  relieve  residual  stress  and  to  damp  the  elastic  energy  as 
discussed  earlier.  A  proposed  interphase  structure  is  shown  in 
Figure  6 . 


FIGURE  6.  MICELLE  STRUCTURE  OF  THE  COPOLYMER 
MICELLAR  ADSORPTION* 

*From  Riess  et  al  (Ref.  8). 

Reproduced  by  permission 


This  PI -SAM  interphase  material  can  couple  to  the  carbon 
fiber  surface  through  the  dipole— dipole  interaction  sketched  in 
Figure  7 . 
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FIGURE  7.  FIXATION  OF  BLOCK  COPOLYMERS  ON  CARBON 
FIBRES  BY  DIPOLE-DIPOLE  BONDING* 


The  introduction  of  this  "soft"  interphase  not  only  increased 
the  shear  strength  of  the  composite,  both  unconditioned  and  con¬ 
ditioned  in  water  at  50 °C,  but  also  improved  the  impact  toughness. 
Some  of  these  results  are  shown  in  Figures  8  and  9. 


FIGURE  8.  SHEAR  STRENGTH  OF  COMPOSITES  AFTER  WATER  TREATMENT 
AT  50 °C* 
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*From  Riess  et  al.  (Ref.  8) . 
Reproduced  by  permission. 


FIGURE  9.  SHEAR  STRENGTH  AND  RESILIENCE  OF: 
RESIN  (EPOXY) 

RESIN  +  UNTREATED  FIBRES 
RESIN  +  GRAFTED  FIBRES* 


*From  Riess  et  al  (Ref.  8). 
Reproduced  by  permission. 


In  reviewing  the  theoretical  concepts  together  with  the 
limited  experimental  data,  a  soft  interphase  between  fiber  and 
matrix  clearly  improves  some  of  the  composite  properties.  The 
objective  of  this  report  is  to  provide  further  understanding  of 
the  role  of  interphase  in  composite  properties  and  to  define 
guidelines  and  principles  for  the  selection  of  optimal  interphase 
materials  to  improve  overall  performance. 
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ROLE  OF  INTERFACE/INTERPHASE  IN  COMPOSITE  TOUGHNESS 

Toughness  in  composite  materials  is  not  a  well  defined 
property.  This  is  due  in  part  to  the  difficulties  in  designing 
test  methods  to  study  the  toughness  of  a  multi-phased  aniso¬ 
tropic  material.  In  general,  toughness  is  measured  by  the 
resistance  of  the  material  to  crack  growth  and  propagation. 
However,  for  multi-phase  anisotropic  materials,  this  is  not 
totally  representative. 

An  increasingly  common  method  in  toughness  measurement  is 
instrumented  impact.  This  method  utilizes  a  load  cell  mounted 
onto  an  impacter.  The  load  cell  is  capable  of  recording  load 
during  impact  in  microseconds.  A  microprocessor  then  outputs  a 
load  deflection  curve  together  with  the  energy  absorption 
profile  during  the  course  of  impact.  The  maximum  force  on  the 
load  deflection  curve  is  taken  as  the  required  load  for  penetra¬ 
tion. 


Due  to  the  multi-phase  characteristic  of  a  composite,  its 
toughness  can  be  improved  in  several  ways.  A  toughened  matrix 
and/or  a  higher  strain  fiber  are  both  capable  of  increasing  the 
toughness  of  a  composite.  Although  little  is  known  about  the 
interface  and  interphase,  it  has  the  potential  to  significantly 
influence  the  toughness  of  composites. 


The  interface  can  affect  the  toughness  of  a  composite  by 
providing  various  energy  absorbing  mechanisms  like  debonding, 
fiber  stress  relaxation  and  fiber  pullout  during  fracture. 

In  the  case  of  weak  interfacial  bonds,  debonding  takes  place  in 
advance  of  the  crack  tip.  The  degree  of  debonding,  as  well  as 
the  energy  absorbed  due  to  the  debonding  process,  depends  largely 
upon  the  interfacial  bond  strength.  In  the  opposite  case  of 
strong  inter facial  bonds,  fibers  will  break  prior  to  debonding 
and  energy  is  absorbed  due  to  stress  relaxation  at  the  ends  of  the 
broken  fibers.  Furthermore,  energy  is  also  absorbed  due  to  some 
fiber  pullout  which  is  related  to  the  strength  of  interfacial 
bonds.  When  all  of  these  energies  are  accounted  for,  it  is 
expected  that  the  weak  interfacial  bond  will  provide  very  high 
fracture  toughness  in  composites.  This  phenomenon  has  been 
observed  by  many  workers. (1“4) 


However,  by  tailoring  the  interphase  properties,  a  signifi¬ 
cant  improvement  in  composite  toughness  can  be  obtained  without 
a  severe  loss  of  mechanical  properties  due  to  weak  interfacial 
bond  strength.  Broutman  et  al(5-7)  have  analytically  demonstrated 
that  interphase  properties  can  be  tailored  to  obtain  both  toughness 
and  strength.  In  general,  a  softer  material  at  the  interphase  is 
preferred . 

Plueddemann (8)  incorporated  a  rubber  interphase  in  glass/ 
epoxy  composites.  Improved  toughness  together  with  relatively 
good  adhesion  strength  was  obtained.  Hancox  and  Wells coated 
different  amounts  of  silicone  rubber  at  the  surface  of  carbon 
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fibers.  A  107%  increase  in  work  of  fracture  without  measurable 
decrease  in  flexural  strength  was  determined  with  1%  by  weight 
silicone  rubber  at  the  interphase  of  a  carbon  fiber/epoxy  com¬ 
posite.  Williams (10)  polymerized  polyvinylalcohol  and  polysul- 
fone  onto  carbon  fiber  and  enhancement  in  shear  strength,  notched 
fracture  toughness  and  low-velocity  impact  toughness  were  achieved. 
Subramanian  and  Jakubowski (11)  coated  carbon  fibers  with  various 
kinds  of  polymers  by  electropolymerization  and  demonstrated  that 
impact  toughness  and  shear  strength  can  simultaneously  increase 
when  optimum  interphase  properties  were  incorporated. 

In  addition  to  the  mechanical  properties  of  the  interphase, 
the  physical  thickness  of  the  interphase  is  another  critical 
parameter  that  needs  to  be  optimized  in  order  to  achieve  the 
ultimate  properties  of  composite  materials.  As  discussed  by  Han- 
cox  and  Wells, (9)  when  silicone  rubber  exceeded  the  1%  by  weight 
level,  a  dramatic  decrease  in  flexural  strength  was  observed. 

Similar  observations  were  also  reported  by  Cavano  and  Winters (12) 
and  Peiffer. (13) 

Another  approach  was  taken  by  Atkins (14)  where  intermittent 
bonding  along  the  fiber  surface  was  suggested.  The  weakly  bonded 
sections  allow  debonding  and  fiber  pullout  during  fracture  while 
the  strongly  bonded  areas  enable  load  transfer.  The  combination 
of  these  could  then  simultaneously  provide  high  fracture  toughness 
and  good  tensile  strength.  Furthermore,  Reiss  et  al(15)  incorpor¬ 
ated  a  two-phase  block  copolymer,  polyisoprene-polystyrene-maleic 
anhydride,  at  the  interphase.  Where  the  polyisoprene  is  providing 
the  energy  absorbing  characteristic  and  the  stryene-maleic  anhy¬ 
dride  block  (SAM)  is  enhancing  the  compatibility  and  joint  strength 
of  the  fiber  with  the  matrix,  increases  in  shear  strength  and 
toughness  were  both  measured.  In  addition,  this  tailored  interphase 
enhanced  the  moisture-conditioned  shear  strength. 

In  view  of  the  results  discussed  above,  composite  toughness 
can  be  increased  without  accompanying  degradation  of  the  remain¬ 
ing  mechanical  properties.  This  can  be  achieved  by  proper  selec¬ 
tion  of  interphase  properties,  such  as  transition  temperature, 
modulus,  strength,  etc.  and  the  physical  characteristics  such  as 
the  quantity  and  thickness  of  this  interphase. 
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EFFECT  OF  MOISTURE  ON  THE  CARBON  FIBER/EPOXY  INTERFACE  AND  INTERPHASE 

Moisture  is  usually  found  on  all  solid  surfaces  and  water 
molecules  are  usually  the  inhibitors/contaminants  for  adhesive 
bonding.  Water  is  generally  attracted  and  adsorbed  onto  the  sur¬ 
face  of  solids  through  hydrogen  bonding  to  the  surface  hydroxyl 
groups .  The  water  molecules  are  normally  one  or  two  monolayers 
deep  on  the  solid  surface  but  this  is  enough  to  induce  stress 
corrosion  and  assist  crack  propagation. 

Evidence  of  adsorbed  water  molecules  at  glass  fiber/epoxy 
interfaces  has  been  shown  and  its  detrimental  effect  on  composite 
properties  were  discussed  by  Bascom. (1)  According  to  Bascom, 
water  molecules  are  attracted  by  the  highly  concentrated  surface 
hydroxyl  groups  on  the  glass  fibers  and  the  reduction  in  composite 
properties  are  due  to  stress  corrosion  at  the  interface  as  well 
as  in  the  bulk  matrix. 

Although  carbon  fibers  also  contain  many  functional  groups 
on  their  surfaces,  moisture  has  a  smaller  effect  on  carbon/epoxy 
composites.  According  to  Gauchel  and  Nash, (2)  only  minor  reduc¬ 
tions  in  composite  properties  are  observed  after  long-term  exposure 
to  water  (12  months) .  Furthermore,  these  property  reductions 
were  primarily  due  to  the  plasticization  effect  of  water  on  the 
epoxy  matrix.  Springer  and  Loos, (3)  and  Augl(4,5)  showed  that 
by  characterizing  effects  of  moisture, on  epoxy  matrix  properties, 
the  moisture  effect  on  carbon/epoxy  composites  performance  can  be 
predicted  with  a  high  degree  of  confidence.  In  other  words,  com¬ 
posite  performance  reduction  is  primarily  due  to  bulk  matrix 
plasticization . 

In  addition,  McKague  et  al(^)  at  General  Dynamics  have  studied 
the  moisture  effect  on  carbon/epoxy  composites  as  a  function  of 
the  sizing  materials.  McKague  concluded  that  the  size  variants 
studied  showed  no  effect  on  total  moisture  absorptivity  (Figures  10 
through  15)  or  glass  transition  temperature  (Table  I),  nor  did  it 
affect  the  resistance  to  thermal  cycling  between  -54 °C  and  124 °C 
( -65 °F  and  255°F) ,  (Table  II). 

In  conclusion,  for  carbon  fiber/epoxy  composites,  the  decrease 
in  shear  strength  after  moisture  conditioning  is  primarily  due  to 
plasticization  of  the  matrix  by  water.  This  is  contrary  to  glass 
fiber/epoxy  composites,  wherein  irreversible  damage  has  been  re¬ 
ported  due  to  moisture  degradation  of  chemical  bonds  at  the  inter¬ 
face.  This  has  not  been  shown  for  carbon  fiber/epoxy  composites. 
Although  moisture  may  have  some  effect  on  the  interphase  region  in 
carbon  fiber  composites,  experimental  techniques  are  not  available 
which  can  separate  the  effects  of  plasticization  of  the  interphase 
from  plasticization  of  the  bulk  of  the  matrix  itself.  Furthermore, 
there  is  no  evidence  of  chemical  degradation  of  the  interphase  region 
by  moisture  at  ordinary  use  temperatures.  Therefore,  no  experimental 
effort  will  be  devoted  in  this  study  to  demonstrate  the  significance 
of  moisture  at  the  interface  and  interphase. 
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CEL  I ON  E000  FINISH 
<y  WITH  NOMINRL  I . I  %  FINISH 
V  WITH  NO  FINISH 


CEL  I DM  E000  FINISH 
❖  WITH  NDMINHL  1.1  %  FINISH 
▼  WITH  0.5  Ji  5210  EPOXY 
(FROM  DRBHNIC  SOLVENT) 


CEL I DM  H000  FINISH 
<>  HITH  NDMINHL  l.l  X  FINISH 
▼  WITH  Z.0  %  £210  EPDXY 
(FROM  DREHNIC  SOLVENT) 


CEL  I ON  E000  FINISH 
V  HUH  NO  FINISH 
<>  WITH  NOMINRL  l.l  X  FINISH 
0  WITH  2.0  %  5210  EPOXY 
(FROM  OREHNIC  SOLVENT) 


TABLE  I.  GLASS  TRANSITION  TEMPERATURES 
ON  INTERPHASE  VARIANTS* 


Sample 

Dry  Tg  °C 

( °F) 

No  Finish 

232 

(450) 

Normal  1.1% 
Finish 

232 

(450) 

0.5%  5210 

Epoxy 

232 

(450) 

1%  5210 

Epoxy 

232 

(450) 

2%  5210 

Epoxy 

232 

(450) 

1%  5210 

Epoxy 

(emulsion)  232  (450) 


Moisture (%) 

Wet  Tg  °C 

( °F) 

0.70 

166 

(330) 

0.67 

166 

(330) 

0.80 

160 

(320) 

0.76 

168 

(335) 

0.78 

166 

(330) 

0.73 

160 

(320) 

Note:  Each  value  is  the  average  of  two  specimens. 

Maximum  pair  difference  was  43.3°C  (110°F) . 


*  From  McKague  et  al  (Ref.  6). 
Reproduced  by  permission. 
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TABLE  II 


THERMAL  CYCLING  DATA* 


0°  Flexural  Properties 
Strength _  _ Modulus 


Uncycled 

Cycled 

Uncycled 

Cycled 

Sample 

MPa 

ksi 

MPa 

ksi 

GPa 

Msi 

GPa 

Msi 

No  Finish 

1931 

280 

1924 

279 

142.1 

20.6 

144.8 

21.0 

Normal  1.1% 

1869 

271 

1952 

283 

140.7 

20.4 

144.8 

21.0 

Finish 

0.5%  5210 

1876 

272 

1897 

275 

138.6 

20.1 

136.6 

19.8 

Epoxy 

1%  5210 

1897 

275 

1938 

281 

135.2 

19.6 

138.6 

20.1 

Epoxy 

2%  5210 

1841 

267 

1814 

263 

131.0 

19.0 

132.4 

19.2 

Epoxy 

1%  5210 

1938 

281 

1986 

288 

142.1 

20.6 

139.3 

20.2 

Epoxy 

(Emulsion) 


*From  McKague  et  al..  (Ref .  6). 
Reproduced  by  permission. 
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ADHESION  THEORY 


Adhesion,  in  physical  chemistry,  is  the  force  of  attraction 
between  a  solid  surface  and  a  second  phase.  The  second  phase  can 
be  a  solid  or  a  liquid.  Adhesion  leads  to  sorption,  which  may  be 
adsorption  on  a  surface  or  adsorption  into  a  surface  layer.  Ad¬ 
hesion  can  be  due  to  electrostatic  forces,  to  van  der  Waals  forces 
or  to  "chemical"  valence  forces.  In  the  latter  case,  one  speaks  of 
"chemisorption."  The  magnitude  of  the  adhesive  forces  is  derived 
from  equilibrium  measurements  of  adsorption  energies  by  thermodynamic 
methods . 

In  adhesion  technology,  only  the  interaction  between  a  solid 
surface  and  a  second  liquid  or  solid  phase  is  termed  adhesion.  The 
technical  process  of  producing  adhesion  between  two  solids  is  called 
adhesive  bonding.  The  bonded  system  is  not  necessary  in  a  thermo¬ 
dynamic  equilibrium;  frequently  the  bonding  process  is  irreversible. 
Therefore,  surface  interactions  cannot  be  separated  experimentally 
from  bulk  properties  of  the  system. 

In  addition,  the  state  of  fundamental  knowledge  concerning  the 
nature  of  the  forces  of  attraction  between  bodies,  determination  of 
the  magnitude  of  such  intrinsic  forces,  and  their  relation  to  measured 
adhesive  joint  strengths  undoubtedly  lags  behind  state-of-the-art 
applied  technology.  One  of  the  main  reasons  why  developing  theory 
has  followed  behind  technology  is  that  the  science  of  adhesion  is  a 
truly  multi-disciplined  subject.  The  adhesion  scientist  often  needs 
to  consider  aspects  of  surface  chemistry  and  physics,  chemistry, 
polymer  chemistry,  polymer  physics,  rheology,  stress  analysis  and 
fracture  phenomena  to  interpret  his  data  fully. 

The  mechanisms  of  adhesion  are  still  not  fully  understood  and 
many  theories  are  to  be  found  in  the  current  literature.  Much  of 
present  confusion  results  from  the  fact  that  test  methods  commonly 
employed  to  measure  the  strengths  of  adhesive  joints  are  not  well 
suited  to  theoretical  analysis.  They  introduce  geometrical  factors 
and  loading  factors  which  are  difficult  to  analyze,  and  the  measured 
joint  strength  includes  indeterminate  contributions  from  rheological 
energy  losses  in  the  adhesives  and  substrate.  Thus,  although  the 
intrinsic  adhesion  forces  acting  across  the  adhesive/substrate  inter¬ 
face  may  affect  joint  strength  they  are  usually  completely  obscured 
by  other  contributions,  and  information  concerning  the  magnitude 
of  such  forces  may  only  be  indirectly  obtained.  This  inability  to 
measure  the  interfacial  interactions  has  been  the  main  obstacle  to 
the  development  of  a  comprehensive  theory  of  adhesion. 
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The  four  principal  mechanisms  of  adhesion  which  have  been  pro¬ 
posed  and  widely  discussed  are: 

1.  Mechanical  Interlocking 

2.  Diffusion  Theory 

3.  Electronic  Theory 

4 .  Adsorption  Theory 

These  are  discussed  in  the  following  paragraphs. 

Mechanical  Interlocking 


This  theory  proposes  that  mechanical  keying  or  interlocking  of 
the  adhesive  into  the  irregularities  of  the  substrate  surface  is 
the  major  source  of  intrinsic  adhesion.  Much  work'  '•  in  this  area 
has  shown  the  importance  of  surface  topography /rugosity  in  intrinsic 
adhesion  forces.  However,  the  observed  increase  in  measured  joint 
strength  with  increasing  surface  rugosity  may  be  attributable  to 
other  mechanisms. 

Diffusion  Theory 


The  diffusion  theory  of  adhesion  is  based  upon  the  intrinsic 
adhesion  of  high  polymers  to  themselves  (autohesion) ,  and  to  each 
other,  due  to  the  mutual  diffusion  of  polymer  molecules  across  the 
interface.  This  requires  that  the  macromolecules  or  chain  segments 
of  the  polymers  (adhesive  and  substrate)  possess  sufficient  mobility 
and  are  mutually  soluble,  that  is,  they  possess  similar  values  of 
the  solubility  parameter.  The  solubility  parameter,  <5,  is  defined 
by 


5 


AH 


V 


RT 


V 


(1) 


where 

AH  =  molar  heat  of  vaporization 

R  =  gas  constant 

T  =  temperature  (°K) 

V  =  molar  volume 

Voyutskii^8)  and  others (9-13)  have  presented  the  evidence  of 
molecular  diffusion  and  their  significant  contributions  to  the  in¬ 
trinsic  adhesion.  However,  with  polymer  systems  which _ are  highly 
crosslinked  and  have  a  high  degree  of  crystallinity,  diffusion 
theory  is  not  likely  to  be  the  principal  mechanism  of  adhesion. 
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Electronic  Theory 


The  electronic  theory  of  adhesion  arises  from  the  different 
electronic  bond  structures  of  the  adhesive  and  the  substrate  which 
cause  some  electron  transfer  on  contact  to  balance  Fermi  levels. 

This  will  result  in  the  formation  of  a  double  layer  of  electrical 
charge  at  the  interface  which  contributes  to  the  intrinsic  adhesion. 
This  theory  was  proposed  by  Deryaguin (l4-!6 )  and  supported  by  other 
workers . (I/-24)  However,  the  arguments  were  not  strong  enough  to 
show  its  significant  contribution  to  the  intrinsic  adhesion.  None¬ 
theless,  surface  electronic  states  should  not  be  overlooked. 

Adsorption  Theory 


The  adsorption  theory  of  adhesion  is  by  far  the  most  generally 
accepted  theory  and  has  been  discussed  in  depth  by  many  authors i2 5 ~2 8 ' 
This  theory  proposes,  provided  that  sufficiently  intimate  inter- 
molecular  contact  is  achieved  at  the  interface,  that  materials  will 
adhere  because  of  the  surface  forces  acting  between  the  atoms  in  the 
two  surfaces .  The  most  common  of  such  forces  are  the  van  der  Waals 
forces  (secondary  bonds).  In  addition,  chemisorption  may  well  occur 
and  thus  ionic,  covalent  and  metallic  bonds  may  operate  across  the 
interface  (primary  bonds) . 

The  thermodynamic  work  of  adhesion  required  to  separate  a  unit 
area  of  two  phases  forming  an  interface,  WA,  may  be  related  to  the 
surface  free  energies  by  the  Young-Dupre  equation . (29 )  In  the  absence 
of  chemisorption  and  inter-diffusion,  the  reversible  work  of  adhesion, 
WA,  in  an  inert  medium  may  be  expressed  by 

WA  =  YX  +  Y2  -  Y12  (2) 

where 

=  surface  free  energy  of  phase  1 
y2  =  surface  free  energy  of  phase  2 
Y12,=  interfacial  free  energy  between  the  phases 

This  equation  only  applies  to  a  solid/liquid  or  liquid/liquid  inter¬ 
face.  However,  by  assuming  that  the  surface  energy  of  a  liquid 
does  not  change  significantly  where  it  solidifies  isothermally  and 
ignoring  the  shrinking  stresses,  it  may  be  applied  to  solid  adhesives/ 
substrate  interfaces . 

Dukes  (20)  an£  others (31-34)  have  been  trying  to  correlate 
measured  strengths  of  adhesive  joints  with  the  thermodynamic  work 
of  adhesion,  WA.  However,  only  some  of  the  authors  have  confirmed 
correlations  while  others  have  found  no  correlations  at  all.  The 
major  cause  of  this  conflict  arises  from  the  use  of  test  methods  for 
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measuring  the  strengths  of  adhesive  joints  which  are  not  well 
suited  to  theoretical  analysis.  Moreover,  correlation  can  only 
be  expected  when  the  locus  of  joint  failure  is  purely  interfacial. 
Therefore,  even  though  the  thermodynamic  work  of  adhesion  affects 
joint  strength,  it  is  often  obscured  by  other  contributions. 

Although  it  has  been  shown  that  secondary  bonding  forces  alone 
can  give  rise  to  acceptable  joint  strengths  at  the  interface,  many 
adhesion  scientists  believe  that  the  additional  presence  of  primary 
bonding  can  increase  the  joint  strength  and  is  certainly  a  must 
for  establishing  environmentally  stable  interfaces.  Evidence  of 
the  presence  of  primary  interfacialy  bonding  has  been  shown  via 
the  use  of  sophisticated,  surface-specific^  analytical  tech¬ 
niques  such  as  Laser -Raman  spectroscopy, (35 )  x-ray  photoelectron 

spectroscopy ,  (36 )  secondary-ion  mass  spectroscopy,  (37 , 38)  ancj 

attenuated  total  reflectance  infra-red  spectroscopy. (39) 

It  was  recognized  by  adhesion  scientists  that  the  establish¬ 
ment  of  an  intimate  molecular  contact  at  the  interface  is  a  necessary, 
though  sometimes  insufficient,  requirement  for  developing  strong  ad¬ 
hesive  joints.  In  other  words,  the  adhesive  needs  to  be  able  to  be 
spread  over  the  solid  substrate  surface  and  needs  to  displace  air 
and  any  other  contaminants  that  may  be  present  on  the  surface. 

Ideally,  an  adhesive  should  conform  to  the  following  conditions:  (40) 

1.  Adhesives  in  the  liquid  state  must  exhibit  a  zero  or 
near-zero  contact  angle  with  the  solid  substrate . 

2.  The  adhesive  should  have  a  relatively  low  viscosity  (a 
few  centipoises)  at  some  time  during  the  bonding  opera¬ 
tion. 

3 .  The  adhesive  should  be  brought  together  with  the  sub¬ 
strate  at  a  rate  and  in  a  manner  that  will  promote  the 
displacement  of  any  trapped  air  at  the  interface. 

In  order  to  assess  the  ability  of  a  given  adhesive/substrate 
combination  to  meet  these  criteria,  it  is  necessary  to  consider 
wetting  equilibria,  to  ascertain  values  of  the  surface  free  energies 
of  the  adhesive  and  substrate  and  the  free  energy  of  the  adhesive/ 
substrate  interface.  In  addition,  the  kinetics  of  the  wetting 
process  should  also  be  examined. 
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WETTING  PHENOMENA 


Surface  tension  is  a  direct  measurement  of  intermolecular 
forces.  The  tension  in  surface  layers  is  due  to  the  attraction 
forces  exerted  on  the  surface  molecules  by  the  bulk  material.  This 
internal  attraction  force  leads  to  a  reduced  number  of  molecules  at 
the  surface  which  result  in  an  increased  intermolecular  distance. 

This  increase  in  intermolecular  distance  requires  external  work  to 
be  done  which  explains  why  surface  tension/surface  free  energy  exist. 

The  most  common  type  of  physical  surface  attractive  forces  are 
van  der  Waals  forces  which  can  be  attributed  to:  (a)  dispersion 

(London-d)  forces  arising  from  internal  electron  motions  which  are 
independent  of  dipole  moments  and;  (b)  polar  (Keesom-p)  forces  aris¬ 
ing  from  the  orientation  of  permanent  electric  dipoles  and  the  in¬ 
duction  effect  of  permanent  dipoles  or  polarizable  molecules.  The 
dispersion  forces  are  usually  weaker  than  the  polar  forces  but  they 
are  universal  and  all  materials  exhibit  them. 

Another  type  of  force  that  may  operate  is  the  hydrogen  bond. 
Hydrogen  bonding  is  the  attraction  between  a  hydrogen  atom  and  a 
second,  small  and  strongly  electronegative  atom  such  as  fluorine, 
oxygen  or  nitrogen.  A  secondary  bond  is  actually  formed  between 
the  hydrogen  atom  and  the  electronegative  (host)  atom  through  an 
electron-pair  on  the  host  species. 

Equilibrium  wetting  occurs  when  the  liquid  and  solid  phases 
are  in  contact  and  can  no  longer  be  perturbed  by  external  forces . 

The  equilibrium  condition  for  the  wetting  of  a  solid  surface  by  a 
liquid  is  a  three-phase  equilibrium  of  solid,  liquid  and  vapor.  As 
shown  in  Figure  16,  one  can  describe  the  equilibrium  point  of 


Vapor  Y 

JjV 


Y 


sv 


FIGURE  16.  EQUILIBRIUM  STATE  OF  A  LIQUID 
DROP  ON  A  SOLID  SURFACE 
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contact  as  the  intersection  of  three  interfaces — solid-liquid  (SL) , 
liquid-vapor  (LV) ,  and  solid-vapor  (SV) .  The  balance  of  these 
surface  forces/surface  tensions  on  this  three-phase  intersection 
is  best  described  by  the  Young-Dupre  equation; (1) 


Ysv  "  ysl  =  ylv  cos  0  (3) 

where 


Ysv 

=  surface 

tension 

at 

ysl 

=  surface 

tension 

at 

ylv 

=  surface 

tension 

at 

0 

=  contact  angle 

solid-vapor  interface 
solid-liquid  interface 
liquid-vapor  interface 


When  0  >90°,  the  liquid  is  non -wetting  and  for  90°>e>0°, 
the  liquid  is  wetting  but  it  is  nonspreading.  When  e  =  0°,  the 
liquid  wets  the  solid  completely  and  spontaneously  and  spreads 
freely  over  the  surface  at  a  rate  depending  on  the  liquid  viscos¬ 
ity  and  solid  surface  roughness.  Nevertheless,  it  is  also  possible 
for  a  liquid  to  spread  and  wet  a  solid  with  contact  angle  e>0°  by 
the  application  of  an  external  pressure  or  force  to  the  liquid. 

To  apply  these  concepts  to  adhesive  technology,  the  surface 
and  interfacial  free  energies  need  to  be  evaluated.  However,  there 
are  no  direct  methods  available  to  measure  surface  free  energy,  ysv/ 
which  complicates  the  analysis.  Nonetheless,  the  contact  angle  0 
is  a  physically  measurable  parameter  and  it  serves  as  a  convenient 
means  for  describing  the  wettability  or  degree  of  wetting  on  solid 
surfaces. 

When  a  liquid  has  a  zero  contact  angle  with  a  solid  substrate, 
surface  tension  gradients  sometimes  exist  at  the  spreading  front 
which  may  either  assist  or  hinder  spreading  depending  upon  their 
directions . (2 ,3)  These  surface  tension  gradients  arise  from  thermal 
gradients ,  or,  in  the  case  of  liquids,  from  the  presence  of  a  more 
volatile  component  of  different  surface  free  energy  (e.g.,  a  trace 
of  impurity).  The  effect  of  the  surface  tension  gradients  on  the 
spreading  rates  are  relatively  small.  However,  it  might  have  an 
effect  on  the  redistribution  of  the  wetting  liquid  after  it  wets  the 
solid  surface. 

. The . topography  of  the  substrate  surface  also  has  an  effect  on 
the  kinetics  of  wetting.  A  liquid  can  spread  along  fine  pores, 
scratches  and  other  inhomogeneities  via  capillary  action.  This 
spontaneous  spreading  occurs  even  if  the  liquid  is  non-wetting  on  a 
smooth,  planar  surface.  A  fifty  per  cent  increase  in  spreading  rate 
was  reported  by  Bascom  et  al'2)  and  Cottington  et  al .  (4)  Again, 
topography  plays  a  minor  rola  in  affecting  the  rate  of  spreading, 
especially  if  the  liquid  is  forcibly  spread  over  the  substrate  sur¬ 
face.  Nonetheless,  topography  does  affect  the  redistribution  of 
the  wetting  liquid  after  its  initial  "wet-out." 
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Much  work (5-12)  has  been  done  to  obtain  quantitative  ex¬ 
pressions  for  wetting  upon  various  types  of  substrate  surface 
topography.  Fundamental  understanding  of  the  importance  of 
wetting  kinetics  upon  wetting  processes  (e.g.,  a  process  to  elim¬ 
inate  trapped  air  at  the  interface  during  wetting)  and  interfacial 
bond  strength  have  yet  to  be  studied. 

Dynamic  wetting  occurs  when  the  liquid  or  the  solid  or  both 
are  kept  in  motion  relative  to  the  other  phase  throughout  the 
wetting  process.  This  process  is  best  described  by  a  sessile  drop 
placed  on  an  inclined  plane  where  two  contact  angles  are  present. 
As  shown  in  Figure  17,  the  advancing  angle  is  the  larger,  frontier 
angle  as  the  liquid  is  progressing  to  "wet"  the  solid  surface. 

The  advancing  angle  reflects  the  tendency  of  the  liquid  to  wet  a 
virgin  surface.  The  receding  angle  is  the  smaller,  tail  end 
angle  where  the  liquid  is  retreating  from  the  "wetted"  surface. 

The  receding  angle  is  a  measure  of  the  ability  of  the  solid  sur¬ 
face  to  remain  wetted. 


FIGURE  17.  A  SESSILE  DROP  ON  AN  INCLINED  PLANE 

(GRADIENT  =  tan<J> )  SHOWING  THE  ADVANCING 

6  AND  RECEDING  6  ANGLES 
^  r 


The  dynamic  wetting  angles  best  characterize  the  condition 
where  the  wetting  liquid  is  forcibly  spread  onto  the  substrate 
surface.  Moreover,  the  advancing  angle  is  very  sensitive  to 
surface  contamination  and  surface  heterogeneity.  This  property 
enables  the  characterization  of  the  substrate  surface  and  provides 
insight  on  the  wetting  process.  The  receding  angle  can  be  related 
to  the  uniformity  of  the  coating  after  the  liquid  has  wetted  the 
substrate  surface. 


41 


In  addition,  wetting  hysteresis  can  be  obtained  by  the  differ¬ 
ence  between  the  advancing  angles  and  the  receding  angles.  The 
hysteresis  of  wetting  is  due  to  the  presence  of  a  large  number  of 
metastable  states  which  differ  slightly  in  energy  at  the  surface 
and  are  separated  from  each  other  by  small  energy  barriers.  These 
metastable  states  are  generally  attributed  to  either  surface  rough¬ 
ness  or  surface  chemical  heterogeneity,  or  both.  Surface  hetero¬ 
geneity  can  also  arise  from  impurities  concentrated  at  the  surface, 
from  crystal  imperfections  or  from  differences  in  the  properties 
of  different  crystal  faces.  In  general,  surfaces  are  both  rough 
and  heterogeneous.  It  is  possible  to  compare  surface  roughness  and 
heterogeneity  by  measuring  the  hysteresis  of  wetting. 

Experimental  Techniques 


The  surface  energetics  of  fibers  can  critically  affect  the 
wetting  behavior  and  the  end-use  performance  of  fibrous  materials, 
e.g.,  in  composite  structures.  The  measurement  of  spontaneous  sur¬ 
face  wetting  is  one  of  the  few  experimental  techniques  available 
for  the  study  of  the  solid-surface  energetics. 

Assessing  contact  angle  measurements  on  filamentous  material 
is  a  very  difficult  experimental  task  due  to  the  fiber  geometry. 

A  preferred  technique  in  direct  contact-angle  measurement  for  fibers 
is  a  derivative  of  the  tilted  plate  concept. (iJ)  This  involves  a 
small  reservoir  of  liquid  which  is  pierced  by  the  fiber  and  an 
arrangement  whereby  the  fiber  can  be  tilted  relative  to  the  liquid 
surface  until  the  meniscus  disappears.  The  contact  angle  is  then 
measured  from  the  fiber  to  the  liquid  surfaces.  However,  considerable 
precaution  must  be  taken  to  obtain  the  true  contact  angle.  This  ex¬ 
perimental  technique  is  extremely  difficult  to  apply  on  fibers  with 
small-cross  sectional  areas. 

An  alternate  experimental  approach  which  does  not  require  the 
direct  measurement  of  contact  angles  is  the  wetting  force  measure¬ 
ment  based  on  the  Wilhelmy  balance  principle . (14 )  This  technique 
is  deemed  to  be  promising  with  filamentous  materials,  especially 
with  small  denier  (fine)  fibers  because  of  the  development  of  high 
sensitivity  microbalances  and  the  ability  to  assess  both  equilibrium 
and  dynamic  wetting  data.  For  fibrous  composites,  dynamic  wetting 
data  is  very  significant  because  of  its  ability  to  prescribe  the 
conditions  for  prepregging  processes. 
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The  wetting  force  measurement,  according  to  Wilhelmy , (14) 
is  the  force  exerted  on  a  vertical  rod  (fiber)  inserted  into  a 
liquid.  This  force  is  expressed  by: 

FT  =  tt  D  yT„  cos  0  (4) 

tyj  L  V 

where 

Fw  =  wetting  force 

D  =  diameter  of  the  fiber 

y  =  surface  tension 
1  LV 

9  =  contact  angle 


A  mathematical  derivation  of  this  expression  was  performed 
by  A.  J.  G.  Allan. (15)  The  balance  of  forces  acting  on  the  fiber 
before  and  after  contact  with  liquid  is  illustrated  in  Figure  18  . 

An  important  feature  of  this  technique  is  its  applicability  to  both 
equilibrium  and  dynamic  conditions.  The  advancing  wetting  force 
can  be  obtained  during  the  fiber  penetrating  into  the  wetting  liquid. 
The  equilibrium  wetting  force  is  obtained  during  the  stable  condition 
when  the  fiber  is  not  in  motion  with  respect  to  the  liquid.  The  re¬ 
ceding  wetting  force  can  be  measured  during  the  pulling  out  of  the 
fiber  from  the  wetting  liquid.  The  corresponding  wetting  force  ex¬ 
pressions  are: 


Ftt  =  t r 
W 

D 

ylv 

cos 

0 

(5) 

F  =  ir 
a 

D 

ylv 

cos 

e 

a 

(6) 

F  =  ir 
r 

D 

ylv 

cos 

0 

r 

(7) 

where 

F„  and  0  are  the  equilibrium  wetting  force 
W  and  contact  angle  respectively 

F  and  0  are  the  advancing  wetting  force 
a  a  and  contact  angle  respectively 

F  and  0  are  the  receding  wetting  force 
r  r  and  contact  angle  respectively 

By  measuring  the  wetting  force,  F^,  and  the  surface  tension 

of  the  wetting  liquid,  y  ,  work  of  adhesion,  Wa  can  be  determined. 

LV  ^ 


wA  = 


LV 


+  tlv  c°s 


(8) 


or 


WA  -  ^LV  +  IH 
ttD 


(9) 
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(a) 


Filament  suspended  (b)  Filament  partially 

freely  immersed 


where : 


Fb  ~  buoyancy  force  ylv  =  surface  tension  of 

M^g  =  gravitational  force  wetting  liquid 

FVj  =  wetting  force  0  =  contact  angle 

D  =  diameter  of  fiber 

assuming  Fb  is  negligible 

then 


7T  D  Yt„  COS  6 

LV 


FIGURE  18 .  FORCES  ACTING  ON  A  SUSPENDED  FILAMENT  BEFORE  AND 
AFTER  CONTACT  WITH  WETTING  LIQUID 
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Furthermore,  by  attaching  a  microprocessor  to  the  micro¬ 
balance  and  controlling  the  fiber  immersion/emersion  operation, 
a  continuous  monitoring  of  the  advancing  wetting  force  and  receding 
wetting  force  along  the  fiber  can  be  achieved.  This  would  enable 
the  determination  of  fiber  surface  heterogeneity  and  provide  in¬ 
sight  into  the  uniformity  of  the  liquid  film  adhered  to  the  fiber 
surface.  Moreover,  this  information  also  shows  the  frequency  of 
occurence  of  critical  surface  contamination  which  prevents  wetting/ 

bonding. 

Critical  surface  contamination  is  present  in  those  areas  of 
the  surface  that  has  a  surface  energy  lower  than  the  wetting  liquid. 
This  lower  surface  energetics  create  a  thermodynamically  unstable 
condition  for  the  higher  surface  energy  liquid  to  "wet"  the  solid 
surface.  Therefore,  a  localized  rejection  of  the  wetting  liquid 
occurs  and  this  leads  to  the  formation  of  voids  during  composite 
fabrication . 

In  addition,  the  hysteresis  of  wetting  can  also  be  studied  by 
this  technique.  A  plot  of  the  immersion  f orce/emersion  force  along 
the  fiber  immersion  length  can  be  used  to  map  the  wetting  hysteresis. 
Figure  19  shows  a  typical  plot  and  identifies  the  hysteresis  of 
wetting. 


FIGURE  19.  HYSTERESIS  OF  WETTING  MEASURED  BY  PLOTTING  ADVANCING 
AND  RECEDING  WETTING  FORCE  VERSUS  IMMERSION 
FIBER  LENGTH 
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Although  the  Wilhelmy  wetting  force  measurement  can  provide 
much  information,  one  experimental  parameter  can  affect  the  accuracy 
of  all  these  data. (1®)  The  rate  of  the  fiber  penetrating  into  and 
receding  out  of  the  wetting  liquid  critically  affects  the  measur- 
ables.  As  shown  in  Figure  20,  a  90°  contact  angle  (total  disappear¬ 
ance  of  the  meniscus)  is  created  by  this  dynamic  phenomena  at  the 
"critical  penetration  velocity."  This  phenomenon  cannot  be  totally 
eliminated.  However,  the  degree  of  meniscus  depression/elevation  can 
be  minimized  and  controlled  by  operating  at  a  much  slower  fiber 
immersion/emersion  rate  than  the  critical  velocity  and  maintaining 
this  rate  throughout  all  studies.  This  technique  has  been  used  by 
many  authors  (17-23 )  an(j  information  on  small  denier  fibers,  such  as, 
carbon  fibers  and  polyester  fibers,  have  been  reported. 


v<v  v=v  v>v 

c  .  c  C 

I  I  i 


where 


V  =  fiber  penetration  velocity 

Vc  =  critical  velocity  at  which  meniscus 
disappears  as  in  (c) 


FIGURE  20.  EFFECT  OF  FIBER  PENETRATION  VELOCITY  ON  CONTACT  ANGLES 


Carbon  Fiber  Surface  Energetics 

Hammer  and  Drzal(l^)  studied  the  surface  energetics  of  Hercules 
Type  A  carbon  fibers  via  the  micro-Wilhelmy  wetting  force  technique. 
Correlations  were  found  between  the  oxygen  concentration  measured  by 
x-ray  photoelectron  spectroscopy  (XPS)  and  the  change  in  polar/dis¬ 
persive  ratios  of  the  carbon  fibers.  This  indicated  that  surface 
treatments  of  carbon  fibers  which  promote  better  fiber  matrix  adhesion 
resulted  from  an  increase  in  the  surface  free  energy.  This  increase 
in  surface  free  energy  was  due  to  the  addition  of  surface  oxygen  con¬ 
centration  during  the  surface  treatment  process.  Some  of  the  contact 
angle  data  are  shown  in  Table  III  while  Table  IV  lists  the  measured 
surface  free  energies  of  carbon  fibers  with  their  polar  and  dispersive 
components . 


TABLE  III. 


CONTACT  ANGLES  FOR  CARBON  FIBERS* 


o 

Liquid  (mJ/m  ) 


Water 

72.8 

Glycerol 

64.0 

Formamide 

58.3 

Methylene  Iodide 

50.8 

Ethylene  Glycol 

48.3 

Br omonaphtha lene 

44.6 

Polypropylene 
Glycol  PG  1200 

31.3 

n-Hexadecane 

27.6 

_  AU  _ 

As  Received  30Q°C  V.T. 


o 

1— 1 

± 

4° 

44° 

+ 

4° 

33° 

± 

2° 

40° 

± 

1° 

37° 

24° 

± 

7° 

23° 

± 

5° 

41° 

+ 

11 

35° 

± 

6° 

30° 

+ 

6° 

29° 

± 

3° 

20° 

± 

4° 

3° 

± 

3° 

25° 

± 

1° 

CO 

o 

± 

8° 

27° 

± 

7° 

_ AS _ 

As  Received  3Q0°C  V.T. 


29° 

± 

4° 

33° 

± 

4° 

25° 

+ 

1° 

42° 

+ 

9° 

22° 

± 

5° 

27° 

± 

3° 

23° 

± 

5° 

30° 

+ 

12' 

31° 

± 

7° 

29° 

± 

2° 

26° 

± 

3° 

21° 

± 

3° 

3° 

± 

3° 

25° 

± 

1° 

22° 

± 

7° 

23° 

± 

7° 

y  =  surface  tension  of  liquid 
LV 

AU  =  Hercules  Type  A  high  strength  carbon 
fiber  non-surface  treated 

AS  =  Hercules  Type  A  high  strength  carbon 
fiber  surface  treated 

V.T.  =  vacuum  treated 


*From  Hammer  and  Drzal  (ref.  19). 
Reproduced  by  permission. 
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TABLE  IV. 


SURFACE  FREE  ENERGIES  OF  CARBON  FIBER* 


Fibers 

Surface  Free  Energies 

As  Received 

P 

Ysv (mJ/m2) 

Ygv (mJ/m2) 

Ysv(mJ/m2) 

AU 

23.6  +  2.6 

27.4  ±  0.3 

51.0  ±  2.6 

AS 

30.0  ±  1.7 

26.4  +  0.1 

56.4  +  1.7 

300°C  V.T. 

AU 

24.1  +  1.4 

26.3  ±  1.1 

50.4  ±  1.4 

AS 

26.8  ±  1.4 

26.0  +  1.2 

52.8  +  1.4 

Ysv  -  Keesom-polar  surface  free  energy  of  solid. 

-  London -dispersion  surface  free  energy  of  solid. 

Ysv  “  Surface  free  energy  of  solid,  ysv  =  YgV  +  YgV 
AU  -  Hercules  Type  A  carbon  fiber — non -surface  treated. 
AS  -  Hercules  Type  A  carbon  fiber — surface  treated. 

V.T.  -  Vacuum  treated. 


*From  Hammer  and  Drzal  (ref.  19). 
Reproduced  by  permission. 
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Sell  ^  studied  the  wetting  properties  of  carbon  fibers  and 
measured  the  effect  of  surface  treatments  on  wetting  properties. 
Figure  21  summarizes  his  reported  results . 


FIGURE  21.  ANGLE  OF  WETTING  0  OF  VARIOUS  LIQUIDS  WITH 
RESPECT  TO  CARBON  FIBERS* 


*From  Sell  (ref.  18) . 
Reproduced  by  permission. 


Kaelble  et  al ^1/ 24-25)  have  intensively  studied  carbon  fiber 
surface  energetics  via  the  micro-Wilhelmy  balance  technique.  Some 
of  his  data  on  the  contact  angles  of  carbon  fibers  are  shown  in 
Table  V.  Kaelble  also  extended  the  information  on  the  Keesom- 
polar  and  London-dispersive  surface  free  energy  with  the  modified 
Griffith  fracture  energy.  This  enabled  the  prediction  of  the  inter* 
facial  conditions  to  promote  good  bond  strength  and  environmental 
stability.  This  expression  was  done  through  the  extension  of  the 
adsorption  theory  of  interfacial  interactions . '  The  critical 
value  of  the  Griffith  fracture  stress,  oc  under  normal  plane  stress 
loading  can  be  expressed  by  the  following  relationships: 


a  = 
c 


/2Ey  \h  /2E\15 

(if)  ■  fc) 


2  2  , 

(R  -R  v-  >0 


2  2 

R  "  Ro 


0.25C  (ax  -  a  3 )  2  +  (l 


(ct2  -  H)2  +  (e2  -  K)2 


H  =  0 . 5  (a  +  a3) 


0.5  (0  +  B3) 


,  d  v  h 


(y  p  )** 
V  1  L  v 


d 


(y  ) 
v  T  svy 


(YP  )% 
T  SVy 


a  2  and  32 


surface  properties  of  the  immersion  phase 
at  the  crack  tip 


For  air:  a2  =  3^  -  0 


For  water:  a2  =  4.67  (mJ/m2) 


=  7.14  (mJ/m2) 
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TABLE  V.  WILHELMY  BALANCE  MEASUREMENTS  OF  HERCULES 
TYPE  AU  AND  TYPE  AS  CARBON  FIBER  AT  23°C* 


AU  AS 


Liquid 

ylv 

(mJ/in2) 

F 

a 

(yg) 

COS0 

a 

6 

a 

w 

a 

(mJ/m2 ) 

Fa 

(yg) 

cos0 

a 

0 

a 

w 

a 

(mJ/m2 

Water 

72.8 

181 

0.8299 

33.9° 

133.22 

172 

0.8966 

26.3° 

138.07 

Glycerol 

64.0 

182 

0.9492 

18.3° 

124.75 

161 

0.9547 

17.3° 

125.10 

Formamide 

58.3 

131 

0.7500 

41.4° 

102.03 

120 

0.7811 

38.6° 

103.84 

Diiodomethane 

50.8 

115 

0.7556 

40.9° 

89.18 

102 

0.7620 

40.4° 

89.51 

Ethylene  Glycol 

48.3 

137 

0.9467 

18.8° 

94.03 

120 

0.9429 

19.5° 

93.84 

1-Bromonapthalene 

44.6 

126 

0.9430 

19.4° 

86 . 66 

111 

0.9445 

19.2° 

86.72 

Glycol  PG-E-1200 

43.5 

13  0^1.  000 

0° 

— 

114 

1.000 

0° 

— 

Glycol  PG-15-200 

36.6 

106 

0.9667 

14.8° 

71.98 

94 

0.9747 

12.9° 

72.27 

Glycol  PG-P-1200 

31.3 

94(1)1.000 

0° 

— 

83(1)1.000 

0° 

— 

n-Dodecane 

25.4 

71 

0.933 

21.1° 

49.10 

63 

0.9413 

19.7° 

49.31 

Ethanol  (abs) 

22.4 

— 

— 

— 

— 

56.3 

0.9413 

19.7° 

43.49 

Y  -  surface  tension  of  liquid 
LV 

Fa  -  advancing  wetting  force 

0a  -  advancing  contact  angle 

W  -  work  of  adhesion 
a 

AU  -  Hercules  type  A  non-surface  treated  carbon  fiber 
AS  -  Hercules  type  A  surface  treated  carbon  fiber 


(1)  Wetting  conditions  with  cosine  0  =  1.0. 

*From  Kaelble  et  al  (ref  21) 

Reproduced  by  permission. 
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Through  this  thermodynamic  analysis,  Kaelble  predicted  that 
in  order  to  have  a  moisture  insensitive  interface,  a  non-polar  sur¬ 
face  for  the  fiber  and  the  matrix  [£■]_,  S3  ^1 . 0  (mJ/m2  )  %]  was  neces¬ 
sary.  To  assure  spontaneous  wetting  to  enhance  interfacial  bond¬ 
ing,  an  increase  in  London -dispersive  contribution  to  the  fiber 
surface  free  energy  >6 . 0 (mJ/m2 ) is  required. 

Penn  et  al(22)  studied  the  advancing  and  receding  contact 
angles  of  carbon  fibers  and  the  effect  of  wetting  hysteresis. 

The  advancing  angle  is  used  to  demonstrate  how  readily  the  liquid 
wets  the  fiber  surface;  the  receding  angle  is  to  illustrate 
how  well  the  fiber  surface  remains  wetted.  The  data  is  presented 
in  Table  VI.  According  to  Penn,  hysteresis  of  wetting  would 
disappear  below  a  critical  surface  tension  of  the  wetting  liquid. 

In  addition,  Penn  also  tried  to  correlate  the  wetting  properties 
of  carbon  fibers  to  composite  shear  strength.  (22)  Penn  sized  the 
carbon  fiber  with  a  silicone  fluid  to  create  a  hydrophobic  surface 
and  also  had  a  hydrophilic  size  to  compare  the  effect  of  sizings 
to  composite  properties  both  unconditioned  and  conditioned. 

Some  of  these  data  are  shown  in  Table  VII.  However,  the  re¬ 
ported  composite  data  do  not  correlate  with  the  wetting  angles 
measured.  The  hydrophilic  sized  fiber  showed  the  smallest  con¬ 
tact  angles  (best  wetting)  but  it  had  the  lowest  short  beam  shear 
strength  both  conditioned  and  unconditioned.  Furthermore,  silicone 
fluid  has  a  poor  compatability  with  epoxy  and  during  prepregging, 
the  silicone  fluid  could  diffuse  away  from  the  fiber  surface  to 
the  prepreg  surface.  This  effect  would  be  more  pronounced  with 
filament  winding  procedures.  This  is  because  the  filament  winding 
resin  systems  are  more  fluid  and  molecule  mobility  is  enhanced. 

The  similarity  of  the  composite  short  beam  shear  strength  of  the 
silicone  fluid  sized  fiber  and  the  unsized  fiber  could  be  the 
reflection  of  the  disappearance  of  the  silicone  size  from  the  fiber 
surface.  Moreover,  Penn  stated  that  silicone  fluid  was  found  on 
the  surface  of  the  laminates.* 

Overall,  much  work  has  been  done  in  measuring  contact  angles 
of  carbon  fibers  and  theoretical  analyses  have  been  made.  However, 
reported  values  are  generally  scattered  and  the  conditions  at 
which  these  measurements  were  made  are  not  generally  clear .  Above 
all,  fiber  origins,  types  and  properties  have  not  been  generally 
characterized  and  specified.  Therefore,  attempting  to  correlate 
the  reported  data  is  difficult  and  unclear.  This  report  is  in¬ 
tended  to  clarify  some  of  the  gray  areas  and  to  try  to  correlate 
wetting  properties  of  carbon  fibers  to  fiber  matrix  adhesion  and 
composite  properties. 


*Private  communication  with  Dr.  Lynn  Penn. 
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TABLE  VI.  ADVANCING  AND  RECEDING  CONTACT 
ANGLES  OF  CARBON  FIBERS* 


ytti  Carbon  Fiber  ^ 

LV  '  "  ' 


(mJ/iu2 ) 

COS0 

0 

cose 

e 

Liquid 

a 

a 

r 

T 

Water 

72.4 

-0.067 

93.8° 

0.713 

44.5° 

Formamide 

55.6 

0.445 

62.9° 

0.562 

55.8° 

Methylene  Iodide 

48.6 

0.824 

34.5° 

0.917 

23.5° 

Ethylene  Glycol 

47.5 

0.538 

57.5° 

0.895 

26.5° 

Br -Naphthalene 

45.3 

0.957 

16.9° 

0.957 

16.9° 

Dimethyl  Acetamide 

36.7 

0.993 

• 

00 

o 

0.993 

6.78 

Hexadecane 

28.4 

0.964 

15.4° 

0.964 

15.4° 

(1)  Origin  and  type 

of  carbon 

fiber 

unknown . 

-  surface  tension  of  liquid. 

LV 

0  -  advancing  contact  angle. 

-  receding  contact  angle. 


*From  Penn  and  Miller  (Ref.  22)  . 
Reproduced  by  permission. 
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STRESS  TRANSFER  FROM  MATRIX  TO  FIBER 


When  a  tensile  load  is  applied  to  an  elastic  matrix,  the 
matrix  will  elongate  to  a  value  proportional  to  the  load.  How¬ 
ever,  if  fiber  with  a  higher  elastic  modulus  is  embedded  in  the 
matrix,  the  fiber  will  restrict  the  matrix  from  freely  elongating 
in  the  region  of  the  fiber.  Shear  stresses  are  produced  on  planes 
parallel  to  the  axis  of  the  fibers  in  the  direction  of  this  axis. 
The  tensile  stresses  in  the  fiber  are  introduced  by  these  inter¬ 
facial  shear  stresses.  These  stress  patterns  can  be  identified 
and  evaluated  by  birefringence  techniques . I1'2' 


A  number  of  models  have  been  proposed  for  calculating  the 
longitudinal  tensile  stress  distribution  in  an  embedded  fiber. 

The  models  of  Cox, <3)  Dow(4)  and  Rosen(5)  considered  the  case  of 
an  elastic  fiber  in  an  elastic  matrix.  If  the  matrix  deforms  elas¬ 
tically,  the  interfacial  shear  stresses  peak  at  a  maximum  value  near 
the  fiber  ends  and  then  decay  rapidly  toward  the  center  length  of 
the  fiber.  Conversely,  the  tensile  stresses  in  the  fiber  will  build 
to  a  maximum  value  provided  that  the  fiber  is  sufficiently  long. 

The  tensile  strain  of  the  fiber  and  the  tensile  strain  of  the  matrix 
are  nearly  equal  in  the  mid-length. 

When  the  load  is  applied  to  the  fiber  by  shear  stress  at  the 
interface,  physically,  at  equilibrium,  one  must  always  have: 


dP 

dx 


TT 


Tx,r  =  35<3f 


(20) 


where 

x  =  distance  from  end  of  fiber 
P  =  tensile  load  in  the  fiber 
d  =  diameter  of  the  fiber 

=  the  shear  stress  evaluated 
Tx ,  r=3sd  at  the  fiber  surface 

Generally,  Tx  r  in  the  matrix  is  a  function  of  both  x  and  r. 
The  magnitude  of  the  interfacial  shear  stress  is  largely  dependent 
on  the  fiber  content,  on  the  magnitude  of  the  applied  load,  and  on 
the  physical  properties  of  the  filament,  matrix  and  inter facial 
bonds . 

In  order  to  obtain  an  expression  for  the  shear  stress,  the 
following  assumptions  are  made.  (3-5) 
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1.  Fiber  and  matrix  both  behave  elastically. 

2.  Fiber  is  surrounded  by  a  cylinder  of  matrix  material. 

3.  The  bond  between  fiber  and  matrix  is  perfect  throughout 
any  deformation. 

4.  Straight  radial  lines  in  the  fiber  and  matrix  remain 
straight  lines  after  deformation. 

The  shear  stress  between  fiber  and  matrix  could  then  be  ex¬ 
pressed  as  follows : (3-5) 


T  x 


Peff 

sinh^— “  [j  “ 

A  E£ 

/XL  \ 

m  f  +  A 

cosh /  - —  ) 

Af 

E  J 

vV  J 

L 

0  <  X  <  - 
2 


(21) 


m 


where : 


X  = 


(22) 


and 


t  =  interfacial  shear  stress  between  fiber  and  matrix 
in  psi 

A  =  cross-sectional  area  in  in2 
E  =  Young's  modulus  in  psi 
G  =  shear  modulus  in  psi 

X  =  distance  from  end  of  fiber,  in  inches 
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df  =  diameter  of  fiber,  in  inches 

L  =  effective  length  for  load  differential  between 
fiber  and  matrix 


aa  =  applied  field  stress  in  psi 

Subscripts : 

f  =  fiber 
m  =  matrix 

Hence,  the  tensile  load  in  the  fiber  at  distance  X  is  given  by:  (3-5) 


Pf  =  ?f  0  +  »  /  Tx  dx 


Therefore,  the  tensile  stress  in  the  fiber  at  position  X  is 
obtained  by  integrating  as  follows : (3-5) 


where : 


The  transfer  length  required  to  build  up  a  maximum  stress 
equal  to  a  fraction  0  of  the  stress  in  an  infinitely  long  fiber 
can  also  be  calculated  by  the  equation: (3-5) 


^  -  ±  oosh  -1  (  -i-  ) 
df  A  \  1-0  / 


(27) 


and  the  load  transfer  pattern  is  demonstrated  in  Figure  22. 
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°f  Hh 
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a.  Elastic  Fiber 
Elastic  Matrix 


b.  Elastic  Fiber 
Plastic  Matrix 


FIGURE  22.  INTERFACIAL  SHEAR  STRESS  PATTERNS  FOR  THE 
ELASTIC  AND  PLASTIC  MATRIX  CASES 
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From  equation  (21),  one  can  see  that  a  maximum  shearing  stress 
occurs  at  the  ends  of  the  fiber.  When  a  high  tensile  load  is 
applied,  the  shear  stress  is  likely  to  exceed  the  yield  stress  of 
the  matrix.  If  the  matrix  is  unable  to  plastically  deform,  either 
the  matrix  fails  or  an  inter facial  slip  occurs  between  the  fiber 
and  the  matrix.  If  the  matrix  is  able  to  deform,  then  the  inter¬ 
facial  shear  stress  never  rises  above  the  yield  shear  strength  of 
the  matrix  and  equation  (21)  becomes  irrelevant.  Therefore,  the 
plastic  flow  model  of  Kelly  and  Tyson*6)  appears  to  be  more  approp¬ 
riate  . 

In  the  case  of  a  plastic  matrix,  a  yield  shear  strength  Ty  is 
constant  throughout  the  deformation.  Then  one  can  obtain  a  tensile 
stress  distribution  in  the  fiber  by  integrating  equation  (20):  (6) 


P  = 


ird^T  dx 
f  y 


=  TTd-X  X 

f  y 


(28) 


Equation  (28)  shows  that  the  tensile  load  builds  up  linearly 
and  a  transfer  length  of  the  fiber  is  obtained.  (.6). 

(29) 


(30) 


a  =  tensile  strength  of  the  fiber  at  transfer  length 

df  =  fiber  diameter 

Ty  =  interfacial  shear  strength 

L  =  critical  length  of  the  fiber  required  for  load  transfer 
from  a  fiber  end 

£  =  critical  fiber  length  required  for  load  transfer, 

c  £  =  2l 

c  c 

and  the  load  transfer  mechanism  for  the  plastic  deformable  .matrix 
is  shown  in  Figure  22 (b) . 


f 
4  T 


or 


T  = 


gfdf 
2  £ 


where 
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METHODS  FOR  MEASURING  INTERFACIAL  BOND  STRENGTH 
Fiber  Pullout  Test 

Broutman^  developed  the  single  filament  fiber  pullout  test 
to  measure  the  interfacial  joint  strength  and  interfacial  frictional 
strength.  A  typical  load -displacement  curve  resulting  from  such 
a  test  is  shown  in  Figure  23.  The  geometry  of  the  specimen  for  the 
fiber  pullout  test  can  either  be  a  button-type  sample  as  shown  in 
Figure  24,  or  a  fiber  end  of  controlled  length  embedded  into  a 
bulk  matrix  as  in  Figure  25.  The  relationship  which  assumes  the 
shear  stress  to  be  uniformly  distributed  along  the  interface  is: 


T 


ct  r 


m 


2irr£  2  Z 


(31) 


where 

t  =  average  shear  strength  of  joint 
=  maximum  load  applied  to  fiber 
r  =  radius  of  fiber 
5,  =  embedded  fiber  length 

=  maximum  stress  applied  to  fiber 


The  embedded  fiber  length  is  influenced  by  the  fiber  strength 
so  that  the  maximum  embedded  length  to  be  used  is  determined  by: 


where 


SL 

max 


°  i  t  r 

ult 


2  T 


(32) 


a  ,  =  ultimated  tensile  strength  of  the  fiber 

ult 

(2) 

The  theoretical  analysis  of  this  test  was  performed  by  Lawrence 
and  Bartoi§(3)  who  concluded  that  this  test  was  sensitive  to  interfacial 
adhesion  properties. 

However,  this  test  is  not  applicable  to  carbon  fiber/epoxy 
joint  strength  analysis  because  of  the  small  diameter  of  the  carbon 
fiber  and  the  relatively  high  interfacial  bond  strength.  Using  the 
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FIGURE  23 .  A  TYPICAL  LOAD-DISPLACEMENT  CURVE 
FOR  FIBER  PULLOUT  TEST* 


*From  Broutraan  (ref.  1) 
Reproduced  by  permission. 
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FIGURE  24.  BUTTON-TYPE  SPECIMEN  FOR  FIBER  PULLOUT  TEST 


FIGURE  25.  BLOCK-TYPE  SPECIMEN  FOR  FIBER  PULLOUT  TEST 


typical  properties  of  carbon  fiber  shown  in  Table  VIII,  the  maximum 
embedded  length  for  carbon  fiber/epoxy  determined  from  equation 
(32)  had  to  be  less  than  lOOum  (3.5  mils). 


TABLE  VIII.  SOME  TYPICAL  PHYSICAL  PROPERTIES  OF  A 
SINGLE  FILAMENT  CARBON  FIBER 


Tensile  strength  of  carbon  fiber  at  2.54  x  10"2m  gauge  length 


ult 


=  3241  MPa  (470  ksi) 


Diameter  of  carbon  fiber,  d  =  7  x  10  6m  (0.14  mil) 
Interfacial  shear  strength,  t  %  70  MPa  (10  ksi) 


Favre  et  al(4)  attempted  to  use  this  technique  to  study  the 
joint  strength  of  carbon  fiber  in  various  thermoset  matrix  materials 
(epoxy,  polyester  and  polyimide) .  A  special  fabrication  technique 
was  designed  to  provide  a  short  embedded  length.  However,  the  thick¬ 
ness  of  the  resin  lozenge  was  difficult  to  control  and  large  varia¬ 
tions  in  the  data  were  observed.  Moreover,  the  critical  embedment 
lengths  of  the  fiber  for  each  matrix  were  not  determined.  The 
critical  embedment  length  could  be  determined  by  testing  specimens 
with  various  embedded  lengths  and  relating  the  fiber  pullout  force 
to  the  embedded  length  as  shown  in  Figure  26.  According  to  the 
authors,  the  joint  strength  of  equation  (31)  should  be  calculated 


FIGURE  26.  RELATIONSHIP  OF  RESIN  LOZENGE  THICKNESS  TO  FIBER 
PULLOUT  FORCE 
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based  on  P 


c  and  £c  as  in  equation  (33) . 


P 

m,c 

ird  l 

c 


(33) 


where : 

Pmc=  fiber  pullout  force  at  maximum  embedded  length 
&c  =  maximum  embedded  length  which  allows  fiber  pullout 

This  relation  was  not  determined  by  Favre  possibly  because  of  the 
difficulties  in  controlling  and  fabricating  ultrathin  resin  lozen¬ 
ges.  Some  of  his  data  is  shown  in  Tables  IX  and  X. 


TABLE  IX.  VALUES  OF  THE  ADHESION  STRENGTH  a 
FOR  AS -RECEIVED  FIBERS*  A 


Fibre 

Type 

Modulus  (GNm-*) 

Matrix  * 

o\  Nmm  -* 

Courtaulds 

HT 

270 

Polyester  A 

5.4 

(long  staple) 

Polyester  B 

10.4 

Epoxy 

9.7 

Polyimide 

17.1 

HTS 

283 

Epoxy 

57.3 

(continuous  tow) 

HM 

330  (1st  batch) 

Epoxy 

5.5 

(continuous  tow) 

Polyimide 

12.8 

357  (2nd  batch) 

Epoxy 

1.8 

HMS 

350 

Epoxy 

31.4 

(continuous  tow) 

Morgan 

II  treated  (sample) 

240 

Epoxy 

41.9 

Carbone-Lorraine 

AG  (short  staple) 

440 

Epoxy 

1.0 

AGT  (continuous  tow) 

258 

Epoxy 

41.6 

Columbia 

CMR  (sample) 

42t 

Epoxy 

8.3 

*Poly ester 

A  Rhodester  1108  ( Rh6ne -Pou 1 enc ) 

Polyester 

B  Stratyl  A.30(St  Gobain) 

Epoxy 

Araldite  LY 

556  +  HT  972 

(CIBA) 

Poly imide 

Kerimid  601 

( Rh^ne -Poulenc ) 

f Manufacturer ' s  data 

From  Favre  and  Perrin  (Ref . 

4) 

Reproduced  by  permission. 
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F 

aA  ird£ 


In  Table  IX, 


where 

F  =  the  force  required  for  fiber  pullout 
d  =  diameter  of  the  fiber 
i  =  thickness  of  the  pastille 
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Single  Filament  Shear  Fracture  Energy  Test 


Outwater  and  Murphy (5“7)  developed  a  new  test  method  to  meas¬ 
ure  the  debonding  fracture  energy  between  epoxy  resin  and  single 
filament  glass  fibers.  The  authors  considered  an  infinitely  long 
fiber  of  diameter  "df"  and  cross-sectional  area  "Af"  embedded 
in  a  semi -inf inite  solid  as  in  Figure  27.  Let  P  be  the  tensile 
load  at  the  end  of  the  fiber  protruding  a  distance  L  from  the  sur¬ 
face  of  the"  resin  and  let  the  fiber  be  debonded  a  distance  "x" 
into  the  matrix  and  GI;[  be  the  energy  required  to  debond  the  fiber . 

The  strain  energy,  in  the  filament  is: 

P2L  /~x(P-nrd,x)2  (34  s 

ft£  =  -  +  /  - - -  dx  1 

2AEf  JO  2AEf 


where 


t  =  interfacial  shear  strength 
Ef  =  Young's  modulus  of  the  fiber 


FIGURE  27 .  EMBEDDED  SINGLE  FIBER  IN  A  SEMI-INFINITE 
SOLID  UNDER  TENSILE  LOAD* 

*From  Outwater  and  Murphy  (Ref.  5) . 

Reproduced  by  permission. 
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By  equating  the  energy  required  to  debond  an  incremental 
length,  dx,  Gjjfrdfdx  and  the  strain  energy  released  by  an  incre¬ 
ment  dx  of  debonding, 


(P-TTrdfX)  2 


where  A  =  cross-sectional  area,  the  following  equation  is 
obtained: 


Gii"df  - 


(p_xird  X) 


Substituting  for  P  =  a* A,  then 


where 


V  d,  /  \  d ./ 


a  =  stress  at  the  fiber  end  not  embedded  in  resin 

Directly  pulling  the  fiber  makes  this  measurement  experimen¬ 
tally  difficult.  A  modified  technique  is  to  embed  the  fiber  in  a 
resin  block,  severing  the  fiber  at  the  center  of  the  block  and 
then  loading  the  specimen  in  compression  and  observing  the  stress 
required  to  debond  the  fiber.  A  diagram  of  this  method  is  shown 
in  Figure  28.  The  G  value  was  determined  from: 


where 


ar  =  applied  stress 
x  =  debonded  length 

and,  for  initial  debonding,  where  x  =0,  then 


e  2E^d 

r  f  f 


where 
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=  compressive  strain  in  the  resin 


FIGURE  28 .  SPECIMEN  FOR  MEASURING  THE  DEBONDING  ENERGY 
OF  THE  FIBER  FROM  THE  SURROUNDING  MATRIX* 


*From  Outwater  and  Murphy  (Ref.  7). 
Reproduced  by  permission. 


The  data  obtained  from  this  test  for  glass  fibers  in  an  epoxy 
matrix  was  conclusive.  However,  due  to  the  small  diameter  and  the 
brittle  nature  of  the  carbon  fiber,  this  test  is  not  applicable 
for  a  study  of  carbon  fiber/epoxy  adhesion  properties.  This  test 
requires  the  observation  of  initial  debonding  during  compression 
of  the  specimen  and  the  small  diameter  of  a  carbon  fiber  makes  the 
observation  difficult  if  not  impossible. 

Single  Filament  Interfacial  Shear  and  Transverse  Tensile  Debonding 
Strength  Determination 


Broutman (1/ 8, 9) has  also  developed  tests  to  characterize  the 
interfacial  shear  debonding  strength  and  interfacial  transverse 
tensile  strength.  Interfacial  shear  strength  is  determined  by 
embedding  a  single  filament  in  the  center  of  an  epoxy  casting  as 
shown  in  Figure  29.  Upon  subjection  to  compression,  the  fiber 


Fiber 

1.27cm 


A-1  .27cm — *| 


3 . 81cm 


it 


FIGURE  29.  INTERFACIAL  SHEAR  STRENGTH  SPECIMEN 


debond  at  the  interface  by  shear.  The  shear  stress  can  be 
determined  from: 


x  =  2.5a 


(39) 


where 
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t  =  interfacial  shear  stress 
a  =  axial  compression  stress  in  specimen 


Interfacial  transverse  tensile  strength  was  determined  in  a 
similar  manner.  The  specimen  for  this  test  is  shown  in  Figure  30. 


FIGURE  30.  INTERFACIAL  TRANSVERSE  TENSILE  STRENGTH  SPECIMEN 


Due  to  the  curvature  of  the  matrix  and  the  difference  in  the 
Poisson's  ratios  of  the  fiber  and  the  matrix,  the  fiber  will 
debond  at  the  interface  through  transverse  tensile  stress.  This 
transverse  tensile  stress  can  be  determined  from: 

_ 1  Ef  (40) 

(1+VEf  + 

where 

S  =  transverse  tensile  stress 
am  =  axial  compression  stress  on  specimen 
y  =  Poisson's  ratio 
E  =  elastic  modulus 

Subscripts:  f  =  fiber 

m  =  matrix 

These  tests  have  been  used  to  evaluate  carbon  fiber/epoxy  ad¬ 
hesion  properties  by  Hawthorne and  Mullin.  However,  they 

both  concluded  that  these  tests  were  not  applicable  to  carbon  fibers 
due  to  the  low  compression  strain  of  the  carbon  fibers.  Compression 
fracture  and  local  buckling  occured  prior  to  the  debonding  of  the 
fiber . 
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Single  Filament  Critical  Length  Determination 


The  single  filament  critical  length  determination  is  an  experi¬ 
mental  technique  derived  from  the  load  transfer  analysis  discussed 
under  "Stress  Transfer  from  Matrix  to  Fiber."  Kelly (16)  illustrated 
the  multiple  fracture  phenomenon  in  a  system  consisting  of  tungsten 
wires  embedded  in  a  copper  matrix.  After  deformation,  the  copper 
matrix  was  dissolved  and  the  wire  fragments  extracted.  On  analysis 
of  the  fragment  length  distribution,  it  was  found  that  its  bounds 
were  approximately  Zc  and  £c/2;  where  £c  was  the  critical  fiber 
length  and  was  determined  from: 


where 


Z 

a 


afd  f 


(41) 


of  =  tensile  strength  of  fiber  at  Zc  gauge  length 

df  =  diameter  of  fiber 

Ty  =  interfacial  shear  strength 

Equation  (41)  is  a  rearrangement  of  equation  (30) .  Recalling  the 
definition  of  critical  fiber  length,  this  result  is  not  unexpected. 
The  longest  fibers  capable  of  surviving  in  this  load  situation  will 
have  a  length  just  under  Z  ;  since  longer  fibers,  by  definition, 
must  fracture. 

Analogous  experiments  were  run  by  Ongchin  et  al, (18)  McGarry(19) 
and  Fraser  et  al(12)  for  glass  fibers  in  both  thermoplastic  and 
thermoset  (epoxy)  matrixes.  Conclusive  results  were  obtained. 
Wardsworth  and  Spilling,  (17)  Ishikawa  et  al(H)  and  Drzal  et  al  (14,15) 
applied  this  technique  to  carbon  fiber  in  an  epoxy  matrix. 

Drzal (14)  embedded  a  single  filament  of  carbon  fiber  into  an 
ASTM  two-and-one-half  inch  length  epoxy  casting  as  in  Figure  31. 

These  specimens  were  then  subjected  to  a  tensile  deformation.  The 
schematic  of  the  fiber  fracture  and  stress  distribution  is  presented 
in  Figure  32  (taken  from  Ref.  14) .  Figures  33  and  34  show  the 
transmitted  polarized  light  micrographs  of  a  non-surface  treated 
(fiber  A)  and  a  surface  treated  (fiber  B)  fiber  during  tensile  de¬ 
formation  and  their  corresponding  shear  strength  data  are  presented 
in  Table  XI.  The  effect  of  carbon  fiber  surface  treatments  on  ad¬ 
hesion  is  shown  clearly  by  the  histogram  of  the  aspect  ratios  of 
the  two  types  of  fibers  in  Figure  35. 


In  addition,  the  failure  mechanisms  during  debonding  can  also 
be  observed  and  characterized. (18)  The  epoxy  sized  carbon  fiber 
promoted  a  higher  interfacial  bond  strength.  The  higher  joint 
strength  created  matrix  cracking  after  fiber  fracture  whereas  no 
matrix  cracking  was  observed  with  the  non-surface  treated,  unsized 
fiber  and  the  surface  treated,  unsized  fiber. 
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FIGURE  31.  SCHEMATIC  DIAGRAM  OF  SINGLE  FIBER  INTERFACIAL 
SHEAR  STRENGTH  SPECIMEN* 


*From  Drzal  et  al  (Ref.  14) . 
Reproduced  by  permission. 


a 


b. 


c . 


a3><T2>CTl 


1TWA 


a.  At  low  initial  stress  values 

b.  After  fiber  flaws  and  defects  have  broken. 

c.  After  critical  length  has  been  reached. 


FIGURE  32.  FIBER  FRACTURE  AND  STRESS  DISTRIBUTION* 

*From  Drzal  et  al  (Ref.  14) . 

Reproduced  by  permission 


FIGURE  33.  TRANSMITTED  POLARIZED  LIGHT  MICROGRAPH  OF  A 
TYPICAL  FIBER  A  FRACTURE  AS  A  FUNCTION  OF 
INCREASING  STRAIN* 


*From  Drzal  et  al  (Ref.  14) . 
Reproduced  by  permission. 
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FIBER  B/EPON  828  -  mPDA 


IOO 


FIGURE  34.  TRANSMITTED  POLARIZED  LIGHT  MICROGRAPH  OF  A 
TYPICAL  FIBER  B  FRACTURE  AS  A  FUNCTION  OF 
INCREASING  STRAIN*. 


*From  Drzal  et  al  (Ref.  14) . 
Reproduced  by  permission. 
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TABLE  XI.  SINGLE  FILAMENT  INTERFACIAL  SHEAR  STRENGTH 
OF  FIBER  A(1)  AND  FIBER  B(2)* 


Fiber  A 


(1) 


Fiber  B 


(2) 


Total  Specimens 

Total  Fragments 

Weibull  a  (eqn.  3) 

Weibull  3  (eqn.  4) 

Fiber  Tensile  Strength 
(1"  gage  length) 

E828/mPDA  Yield  Strength 


15 

23 

669 

485 

3.074 

3  .308 

48.50 

109.05 

2959 

(429,000 

MPa 

psi) 

2614 

(379,000 

MPa 

psi) 

40 

(5,800 

MPa 

psi) 

40 

(5,800 

MPa 

psi) 

Interfacial  Shear  Strength  (t)  40.86±26.21  MPa 
(eqn.  1)  (5,925±3800  psi) 


15 . 6±8 . 89  MPa 
(2 , 262±1290  psi) 


(1)  As— 1 — surface  treated  carbon  fiber,  unsized. 

(2)  AU— 1 — virgin  carbon  fiber  with  no  surface  treatment  and  unsized. 


*From  Drzal  et  al  (Ref  14) . 
Reproduced  by  permission. 
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Number  breaks 


FIbER  R  „t.  FIbER  b 


1/ d  Aspect  Ratio 


Total  number  of  breaks  as  a  function  of  length 
to  diameter  ratio. 


FIGURE  35.  HISTOGRAMS  OF  ASPECT  RATIOS  FOR  FIBER  A  AND  B* 


*From  Drzal  et  al  (Ref.  14) . 
Reproduced  by  permission. 


Overall,  this  test  has  shown  some  conclusive  results  and  pro¬ 
vided  much  information  on  the  characterization  of  fiber/matrix  ad¬ 
hesion.  This  technique  has  been  adapted  by  the  author  to  further 
study  the  effect  of  interphase  properties  on  carbon  fibers  epoxy  ad¬ 
hesion. 


Microdebonding  Test 


Mandell  and  McGarry(20)  developed  a  microdebonding  test  for 
in-situ  measurement  of  fiber/matrix  bond  strength  in  fiber  com¬ 
posites.  The  test  involved  the  compressive  loading  of  a  fiber  on 
a  polished  specimen  surface  to  produce  debonding  as  shown  in  Fig¬ 
ure  36.  Glass,  aramid  and  carbon  fibers  in  epoxy  composites  have 
been  tested.  Initial  results  of  these  evaluations,  especially 
glass,  on  unconditioned  and  moisture  conditioned  samples  were 
promising.  For  carbon  fiber,  experimental  difficulties  were  encoun¬ 
tered  due  to  the  skin-core  properties  of  the  fiber.  Refinements  to 
simplify  interpretation  of  the  debonding  force  in  terms  of  inter¬ 
facial  shear  strength  and  improvement  on  the  apparatus  were  required 
in  order  to  make  the  test  more  reproducible  for  carbon  fibers . 


FIGURE  36.  MICRODEBONDING  TEST* 

*From  Mandell  and  McGarry  (Ref.  20) 

Reproduced  by  permission. 
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EXPERIMENTAL 


INTRODUCTION 

After  reviewing  the  information  available  in  the  literature, 
it  is  apparent  that  little  effort  has  been  spent  on  the  correlation 
of  theoretical  analyses  with  experimentally  determined  composite 
properties;  and  there  is  limited  evidence  that  the  nature  of  the 
interphase  can  be  used  to  optimize  composite  properties.  In  order 
to  clarify  some  of  these  uncertain  areas,  the  following  experiments 
were  designed  to  examine  the  relationship  between  wetting  and  in^ 
terfacial  shear  strength.  Also,  various  finish  variants  were 
formulated  to  demonstrate  the  effect  of  interphase  properties  on 
composite  performance. 


MATRIX  RESIN 

For  carbon  fiber  composites,  diglycidyl  ether  of  bisphenol  A 
( DGEBA )  epoxies  are  commonly  used  in  the  matrix  resin.  In  order 
to  have  a  typical  matrix  that  is  suitable. for  both  critical 
length  determination  and  composite  fabrication,  system  "F"  was 
formulated  for  the  evaluations  reported  herein.  This  resin  is 
prepared  by  mixing  100  parts  of  Epi-Res  508*  and  22.5  parts  of 
Epi-Pure  841.*  A  unique  feature  of  this  system  is  its  high  ten¬ 
sile  elongation-to-break  value  of  8%.  A  cure  schedule  of  two 
hours  at  93°C  was  used  combined  with  a  post-cure  at  150°C  for  two 
hours.  The  tensile  properties  of  system  "F"  castings,  evaluated 
at  different  cure  schedules,  are  shown  in  Table  XII.  This  high 
elongation-to-break  property  of  system  "F"  allows  multiple  frac¬ 
ture  of  the  embedded  fiber  without  matrix  failure  which  enables 
critical  length  determination.  Some  of  the  other  properties  ob¬ 
tained  from  castings  of  system  "F"  are  presented  in  Table  XIII. 

Finish  Variants 


To  study  the  effect  of  fiber/matrix  interphase  on  composite 
performance,  four  finish  variants  were  developed  and  applied  to 
Celion®  6000  (C-6K)  unsized  fiber.  These  four  variants  are: 

1.  C-6K/SR  -  RTV  615**silicone  rubber 

2.  C-6K/FRE  25  -An  organophosphazene***elastomer  blended  with 

Epon  828  at  25  parts  elastomer  to  100  parts 
epoxy 


*  Product  of  Celanese  Plastics  and  Specialties  Co.,  Louisville,  Ky. 

**  Product  of  General  Electric  Co. 

***Product  of  Firestone  Co. 
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TABLE  XII.  TENSILE  PROPERTIES  OF  SYSTEM  "F"  AT 
VARIOUS  CURE  SCHEDULES 


Composition: 

EPI-RES  508*  100  parts  by  weight 
EPI-CURE  841*  22.5  parts  by  weight 


Cure  Post  Elonga- 

Schedule  Cure  Tensile  tion-to 

Time  Temp  Time  Temp  Strength-  ""Modulus  break 


Sample 

w 

TcT 

w 

T*cT 

MPa 

ksi 

GPa 

Msi 

% 

1 

2 

85 

2 

15  0 

82.1 

11.9 

3.02 

0.44 

5.56 

2 

2 

85 

3 

150 

79.3 

11.5 

2.92 

0.42 

5.22 

3 

2 

85 

4 

150 

82.1 

11.9 

2.90 

0.42 

5.77 

4 

2 

85 

2 

160 

84.8 

12.3 

2.97 

0.43 

6.78 

5 

2 

85 

2 

170 

79.3 

11.5 

2.87 

0.42 

5.36 

6 

2 

85 

2 

180 

77.2 

11.2 

2.88 

0.42 

3.20 

7 

2 

93 

2 

150 

92.4 

13.4 

3.15 

0.46 

7.99 

8 

2 

93 

3 

150 

89.0 

12.9 

3.01 

0.44 

5.89 

9 

2 

93 

4 

150 

91.7 

13.3 

2.93 

0.43 

7.60 

10 

2 

93 

2 

160 

80.0 

11.6 

2.93 

0.43 

5.75 

11 

2 

93 

2 

180 

80.0 

11.6 

2.88 

0.42 

5.55 

*Celanese  Plastics  and  Specialties  Co. 
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TABLE  XIII.  RESIN  CASTING  PROPERTIES  OF  SYSTEM  "F"* 


Tensile 


Strength 

Modulus 

MPa 

ksi 

GPa 

Msi 

RT 

92.4 

13.4 

3.15 

0.456 

121°C 

39.3 

5.7 

1.93 

0.280 

3  Day  I^O  boil, 
tested  at  RT 


Elonga¬ 
tion  -to  Flexural 

break  Strength  Modulus 


o. 

o 

MPa 

ksi 

GPa 

Msi 

7.99 

102.8 

14.9 

2.90 

0.420 

9.71 

72.4 

10.5 

2.01 

0.292 

97.2 

14.1 

2.60 

0.377 

Transition  Temperature  Determined  by  TMA : 

Initial  Tg  =  176°C 
After  3  day  H~0  Boil  =  142 °C 
After  7  day  H^O  Boil  =  143 °C 


Heat  Deflection  Temperature  =  155 °C 
[@  1.82  MPa  (264  psi)] 


Notched  Izod  =  58.7  J/m  of  notch 

(1.10  ft-lbs/in  of  notch) 

*Composition :  Epi-Res  508,  100  parts  by  weight 

Epi-Cure  841,  22.5  parts  by  weight 

Cure  Schedule:  2  hours  @  93 °C 

2  hours  @  150 °C 
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3.  C-6K/ETP  10  -  An  epoxy-thermoplastic  blend  (in-house  devel¬ 

oped  material)  with  10  parts  thermoplastic  in 
100  parts  epoxy. 

4.  C-6K/ETP  50  -  Same  as  item  3  except  this  consisted  of  50 

parts  thermoplastic  in  100  parts  epoxy. 

Silicone  rubber  (a  prepolymer)  was  dissolved  in  n-hexane, 
then  the  curing  agent  was  added  to  the  solution  before  sizing. 
The  silicone  rubber  was  mixed  at  a  ratio  of  10  parts  silicone 
rubber  prepolymer  to  1  part  curing  agent. 


The  FRE  25  sizing  solution  was  prepared  by  dissolving  the 
organophosphazene  prepolymer  in  acetone.  Then,  Epon  828  was 
added  in  the  ratio  of  25  parts  elastomer  to  100  parts  epoxy. 

Just  before  sizing,  1%  by  weight  of  the  elastomer  of  Esperox  10 
(t-butyl  perbenzoate)  was  added  to  the  solution.  This  organoper- 
oxide  was  used  as  a  curing  agent  for  the  organophosphazine  elas¬ 
tomer  . 

ETP  50  and  ETP  10  were  dissolved  in  acetone  and  the  concen¬ 
tration  of  all  the  sizing  solutions  were  adjusted  together  with 
the  sizing  speed  to  allow  approximately  1%  by  weight  pickup  by 
the  fibers . 

Silicone  rubber  was  selected  to  demonstrate  the  effect  of  a 
very  soft  interphase  together  with  a  weak  interfacial  joint 
strength.  FRE  25  was  used  to  illustrate  the  effect  of  a  second, 
soft  and  dispersed  phase  at  the  interphase.  Organophosphazene 
was  selected  because  of  its  high  temperature  performance  combined 
with  a  relatively  high  strength  for  an  elastomeric  material;  it 
is  also  tough  and  is  compatible  with  epoxy  resins.  The  purpose 
of  the  polyblend  material  is  to  maintain  the  strength  at  the  inter¬ 
face  while  the  elastomeric  material  is  damping  the  impact  and/or 
blunting  crack  tips  to  enhance  toughness. 

The  same  principles  were  applied  in  the  selection  of  the  thermo- 
plastic-epoxy  blend,  ETP  10  and  ETP  50.  The  only  difference  is 
that  the  thermoplastic  epoxy  blends  provide  higher  modulus  as  well 
as  strength.  The  difference  in  concentrations  is  to  provide  step¬ 
wise  increases  in  modulus  and  strength  at  the  interphase. 

Although  the  properties  of  the  finish  variants  were  not  charac¬ 
terized  due  to  many  unsuccessful  attempts  to  make  castings  from 
these  materials,  it  is  estimated  that  the  moduli  of  the  interphases 
should  be  in  the  following  order: 

ETP  10  >  ETP  50  >  FRE  25  >  SR 
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EXPERIMENTAL  RESULTS  AND  DISCUSSION 
Wetting  Force  Measurement 


The  Wilhelmy  micro-balance  wetting  force  technique  was  used 
to  characterize  the  surface  of  carbon  fibers.  A  micro-computer 
was  used  to  control  the  wetting  experiment  as  well  as  to  carry  out 
the  data  processing.  The  details  of  the  apparatus  construction 
are  discussed  in  Appendix  B  and  the  program  developed  by  the  authors 
for  an  HP-85A  micro-computer  is  attached  in  Appendix  C. 

To  characterize  the  surface  energetics  of  the  carbon  fiber, 
double  distilled  water  was  used  initially.  A  very  careful  cleaning 
procedure  for  the  glassware  which  contains  the  wetting  liquid  was 
required  so  that  the  surface  tension  of  the  contained  water  approaches 
72.4  x  10-3  N/m.  This  vigorous  cleaning  process  is  detailed  in 
Appendix  D.  By  employing  this  cleaning  process,  a  consistent  72.3  x 
10"3  N/m  ±  1  x  10“4  N/m  surface  tension  for  water  was  measured.  Con¬ 
tact  angles  for  unsized  fibers  were  measured  successfully  with  water 
as  the  wetting  liquid.  However,  when  the  fibers  contained  finish, 
misleading  wetting  forces  were  obtained.  This  is  primarily  due  to 
the  finish  slowly  dissolving  into  and  contaminating  the  water.  As 
a  consequence,  a  significant  drop  in  surface  tension  of  the  water 
took  place.  This  caused  the  wetting  force  to  drop  during  the  ex¬ 
periment.  Some  of  these  results  are  listed  in  Table  XIV.  The  large 
scatter  of  data  within  each  region  with  little  or  no  hysteresis 
effects  indicates  that  the  surface  tension  of  the  wetting  liquid 
was  changing  during  the  experiment.  The  surface  tension  of  the 
water  dropped  as  much  as  5  x  10-3  N/m  after  experiments  were  conducted. 

Under  this  condition,  water,  as  well  as  other  commonly  used 
liquids  with  known  surface  tension  are  not  applicable  for  this  study 
because  of  the  contamination /dilution  process.  Nevertheless,  in 
order  to  understand  the  wetting  behavior  and  to  prescribe  the  optimum 
prepregging  condition,  a  major  component  of  the  matrix  resin  should 
be  used  as  the  wetting  liquid.  Although  contamination  by  finish 
from  the  sized  fiber  still  continues,  this  is  representative  of  the 
prepregging  process.  A  DGEBA  epoxy  resin,  Epi-Res  508,  was  selected 
to  be  the  wetting  liquid.  Epi-Res  508  was  chosen  because  of  its 
relatively  low  viscosity  and  it  is  the  major  component  of  the 
matrix  used  in  this  study.  A  viscosity-temperature  profile  of  Epi- 
Res  508  was  determined  in  a  Brookfield  viscometer  and  is  shown  in 
Figure  37.  An  insulated  clean  beaker  was  used  to  contain  and  to 
maintain  the  resin  at  approximately  70°C.  At  this  temperature,  a 
reasonable  low  viscosity,  0.16  Pa  (160  cps)  was  obtained  and  no  sig¬ 
nificant  dynamic  wetting  phenomenon  was  observed.  The  wetting  force 
data  obtained  for  a  series  of  fibers  and  finishes  are  shown  in 
Table  XV. 


Comparing  the  data  in  Table  XV,  the  effect  of  finish  variants 
on  wetting  behavior  was  demonstrated.  The  receding  angles  of  all 
the  samples  are  all  similarly  high  suggesting  that,  once  the  fibers 
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TABLE  XIV.  WETTING  FORCE  MEASUREMENTS  OF  CARBON  FIBERS 
WITH  WATER  AT  23 °C 


Advancing,  yg  Equilibrium,  yg  Receding,  yg 

Sample  (contact  angle,  e  )  (contact  angle,  9)  (contact  angle,  9_r) 


C-6KU(1) 

82.3 

± 

5.5 

89.4 

± 

7.2 

97.5 

± 

9.7 

(60° 

± 

2.3°) 

(57° 

± 

3°) 

(53° 

± 

4.3°) 

as-i(2) 

81.7 

± 

7.5 

91.3 

± 

9.8 

108.0 

± 

8.0 

(61° 

+ 

3°) 

(57° 

± 

4°) 

(50° 

± 

3.3°) 

C-6KE(3) 4 

44.2 

± 

14.6 

41.8 

± 

24.0 

63.5 

± 

10.0 

T-6300  ^ 

90.9 

± 

18.3 

91.3 

± 

23.3 

100.4 

± 

21.4 

(1)  Celion®  6000  unsized  carbon  fiber. 

(2)  Hercules  AS-1  unsized  carbon  fiber. 

(3)  Celion®  6000  standard  epoxy  compatible  sized  fiber, 
commercial  product. 

(4)  Thornel  T-6300  standard  commercial  sized  fiber. 
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(xlO3  cps) 


Temperature  (° 


FIGURE  37.  VISCOSITY-TEMPERATURE  PRO] 
EPOXY  RESIN  MEASURED  BY  A 
WITH  SPINDLES  NO.  1  AND  2 


TABLE  XV.  WETTING  FORCE  MEASUREMENTS  OF  CARBON  FIBERS 
WITH  EPI-RES  508  RESIN  AT  70 °C 


Sample 
T-6300  (1) 

as-i(2) 3 

C-6KV^ 
C-6KU(4) 
C-6KE(5) 6 
C-6K/ETP  50 (6^ 
C-6K/ETP  10(7) 8 9 
C-6K/FRE  25  ^ 
C-6K/SR^9) 


Advancing,  yg 


72.8 

± 

27.6 

57.2 

± 

10.6 

67.0 

± 

11.3 

75.6 

± 

19.4 

55.8 

+ 

19.5 

107.4 

± 

28.4 

86.0 

± 

40.0 

71.7 

± 

10.7 

-1.6 

± 

18.2 

Equilibrium,  yg 


95.6 

± 

18.5 

74.4 

± 

5.4 

80.2 

± 

9.9 

83.2 

± 

19.5 

91.4 

+ 

13.3 

127.2 

± 

29.4 

103.0 

± 

39.0 

82.0 

± 

12.6 

50.6 

+ 

4.1 

Receding,  yg 

103.8  ±  20.2 

87.2  ±  7.8 

92.2  ±  14.9 
88.6  ±  20.9 

102.4  ±  12.9 
135.0  ±  29.1 
109.0  ±  43.0 

90.2  ±  14.9 
99.8  ±  48.2 

C-6KV. 


Note:  All  carbon  fibers  were  surface  treated  except 


(1)  Thornel  standard  sized  6000  filament  per  tow  carbon  fiber. 

(2)  Hercules  unsized  carbon  fiber. 

(3)  Celion®  6000  virgin  fibers — no  surface  treatment  and  unsized. 

(4)  Celion®  6000  unsized  carbon  fiber. 

(5)  Celion®  6000  standard  epoxy  compatible  size,  commercial  product. 

(6)  Celion®  6000  sized  with  finish  variant  ETP  50. 

(7)  Celion®  6000  sized  with  finish  variant  ETP  10. 

(8)  Celion®  6000  sized  with  finish  variant  FRE  25. 

(9)  Celion®  6000  sized  with  finish  variant  SR. 
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are  "wetted,"  there  is  little  or  no  rearrangement  of  the  resin 
film  on  the  surface  of  the  fiber.  However,  in  the  advancing 
measurements,  large  variations  are  observed.  This  suggests 
that  some  finish  variants  have  better  compatibility  with  the 
matrix  resin  and  can  be  wetted  easier.  For  example,  silicone 
rubber  coated  Celion®  6000  showed  the  lowest  wetting  force 
(negative)  which  indicates  that,  in  order  to  wet  this  surface, 
additional  external  forces  are  needed  to  forcibly  wet  it.  This 
will  most  likely  also  result  in  poor  adhesion. 

On  the  other  hand,  the  finish  variant  ETP50  showed  the  highest 
advancing  wetting  force  which  suggests  that  the  resin  is  highly 
compatible  with  this  surface  and  it  will  probably  translate  to 
good  interfacial  adhesion.  Furthermore,  with  a  high  advancing 
wetting  force  fiber,  better  quality  prepreg  is  likely,  due  to 
the  minimization  of  trapped  air  at  the  interface  by  rapid  advance¬ 
ment  of  the  wetting  liquid/resin.  In  other  words,  improved  lam¬ 
inate  properties  are  anticipated. 

Single  Filament  Adhesion  Test 


As  previously  discussed,  the  single  filament  critical  length 
determination  was  identified  to  be  one  of  the  better  methods  for 
characterization  of  carbon  fiber/epoxy  interfacial  bond  strength. 
Recall  the  load  transfer  mechanism,  wherein  the  critical  length, 
i  is  determined  by  equation  (41)  : 

C 


a 

c 


(41) 


where 

a  is  the  tensile  strength  of  the  fiber  at  critical 
length 

ci 

f  is  the  diameter  of  the  fiber 

is  the  interfacial  shear  strength 


By  determining  Of  and  measuring  ac  and  df,  the  interfacial  shear 
strength,  Ty  can  be  obtained  as  already  stated  by  rearranging 
equation  (41)  to  the  form  shown  for  equation  (30) : 


(30) 
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The  critical  length  and  the  diameter  of  the  fiber  can  be  measured 
optically  by  embedding  a  single  filament  in  an  epoxy  matrix  and 
stressing  the  specimen  to  a  strain  value  such  that  no  additional 
fracture  can  be  obtained  with  additional  stress.  However,  to  ob¬ 
tain  the  tensile  strength  of  carbon  fiber  at  the  critical  length 
is  itself  an  involved  task.  Carbon  fiber  strength  is  sensitive 
to  gauge  length  and  the  critical  length  of  carbon  fiber  is  too  small 
to  experimentally  test  in  a  reliable  manner.  Nevertheless,  it  was 
shown (1)  that  a  linear  relationship  can  be  obtained  by  plotting 
single  filament  tensile  strength  versus  the  logarithm  of  gauge 
lengths.  The  test  method  is  described  in  Appendix  E. 

Gauge  lengths  of  7.62  cm  (3 -inch) ,  2.54  cm  (1-inch) ,  1.27  cm 
(%-inch) ,  0.635  cm  (J?-inch)  and  0.254  cm  (1/10-inch)  were  tested 
for  each  fiber  type  and  finish  variant  and  best  fit  straight  lines 
were  obtained  from  these  five  data  points  for  each  variant.  These 
straight  lines  were  then  extrapolated  to  the  critical  length  and 
the  tensile  strengths  were  obtained.  The  results  of  these  gauge 
length  studies  are  presented  in  Appendix  F . 

To  perform  the  adhesion  test,  a  single  filament  of  carbon  fiber 
was  aligned  in  the  center  of  an  ASTM  6.35  cm  (2h  inch)  dogbone  epoxy 
casting.  The  dimensions  of  the  dogbone  specimens  are  0.32  cm 
(1/8  inch)  wide  x  0.16  cm  (1/16  inch)  deep  with  a  2.54  cm  (1  inch) 
long  gauge  section.  Fabrication  of  the  specimen  and  alignment  of 
the  fiber  was  accomplished  with  the  aid  of  an  RTV-630*  silicone 
rubber  mold.  A  silicone  rubber  mold  was  used  because  of  the  ease 
in  making  a  flexible  mold  cavity  by  casting  the  rubber  around  stand¬ 
ard  ASTM  injection-molded  dogbone  specimens.  Moreover,  the  use  of 
silicone  rubber  molds  also  prevents  fiber  surface  contamination  be¬ 
cause  they  do  not  require  the  application  of  mold  release  agents. 

Sprue  spots  were  cut  precisely  in  the  center  of  each  end  for 
each  dogbone  to  a  depth  of  0.08  cm  (1/32  inch)  to  align  and  secure 
each  fiber  to  be  embedded.  Once  the  fiber  was  aligned,  a  drop  of 
rubber  cement  was  used  to  hold  it  in  place.  During  the  fiber  selec¬ 
tion,  mounting  and  the  application  of  rubber  cement,  the  section  of 
the  fibers  placed  in  the  mold  cavity  were  not  handled  nor  contamin¬ 
ated. 


After  the  rubber  cement  was  set,  system  "F"  resin  was  intro¬ 
duced  with  an  eye  dropper  into  the  cavity.  Extreme  care  was  exer¬ 
cised  to  prevent  damaging  and/or  mis-aligning  the  fiber.  Then  the 
whole  was  degassed  in  a  vacuum  oven  and  cured  according  to  the 
schedule  previously  described. 

Since  system  "F"  is  an  amber  transparent  solid,  visual  and 
microscopic  evaluation  of  the  fiber  and  the  critical  length  can  be 
performed  with  ease.  Furthermore,  system  "F"  also  exhibits  photo¬ 
elastic  properties  which  can  be  used  to  evaluate  the  stress  patterns 
at  the  in ter phase. 

*General  Electric  Co.,  Silicon  Rubber  Dept.,  Waterford,  NY.  ~ 
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To  apply  this  test  to  the  above  mentioned  specimen,  a  special 
tensile  fixture  was  designed  and  mounted  onto  a  stage  micrometer 
for  a  Leitz  ortholux  optical  microscope.  Long  working  distance  ob¬ 
jectives  and  polarizing  condensers  were  used.  The  details  of  the 
design  for  the  tensile  fixture  are  presented  in  Appendix  G. 

To  carry  out  this  study,  the  specimen  was  initially  examined 
for  fiber  alignment,  breakage  and  residual  stress.  Due  to  the  low 
cure  and  post  cure  temperatures  employed  with  this  system,  together 
with  a  slow  cool  down  rate,  no  apparent  residual  stresses  were  ob¬ 
served  in  the  samples.  The  specimen  was  then  stressed  to  different 
strain  values  at  approximately  0.5%  intervals  and  at  each  strain 
level,  the  number  of  the  broken  segments  was  counted.  When  this 
number  remains  constant  with  increasing  strain,  the  critical  length 
of  the  fiber  is  reached.  This  test  was  repeated  on  at  least  ten 

specimens  for  each  fiber  type/finish  variant  and  the  cumulative 
fiber  lengths  were  used  to  determine  the  critical  length. 

Once  the  segment  lengths  of  the  fibers  were  obtained,  a 
statistical  approach  to  determine  the  critical  length  of  the  fiber 
was  used.  According  to  Drzal  et  al,(2)  these  data  sets  fit  the 
two  parameter  Weibull  distribution  auite  well. 

F(X)  =  l-exp-|-(^)a|  for  X>o  (42) 

where 

X  is  the  aspect  ratio,  £  /d^ 

a  is  the  shape  parameter 

3  is  the  scale  parameter 


The  maximum  likelihood  methods  were  employed  to  estimate  the 
parameters  of  the  Weibull  distribution.  These  parameters  are 
solutions  to 


(44) 
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A  computer  program  was  developed  to  evaluate  these  two 
parameters,  a,  3,  numerically,  theso  were  then  used  to  calculate 
the  shear  strength  and  the  listing  is  attached  in  Appendix  H. 

The  mean  and  the  variance  of  the  interfacial  shear  strength  were 
then  obtained  from 

?  T  ( 1  -  ±  ) 

23  y  a  ) 


and 

2  , 

VAR (x )  =  )  r(1_  2 

432  \ 


where 


F  is  the  Gamma  function 


(45) 


)  -  (i- 


(46) 


Another  program  was  also  developed  to  plot  the  cumulative 
function  and  the  population  density  function  (histogram)  of  the 
critical  aspect  ratio  for  each  fiber  type/finish  variant.  This 
program  is  also  listed  in  Appendix  I. 

The  results  of  the  interfacial  bond  strength  determinations 
for  the  different  fiber  types/finish  variants  are  listed  in 
Table  XVI  and  the  histograms  of  the  critical  aspect  ratios  are 
shown  in  Figures  38  through  46.  The  results  clearly  indicate 
that  the  surface  treatment  improved  the  interfacial  adhesion 
significantly  and  that  coating  the  carbon  fibers  with  silicone 
rubber  inhibits  interfacial  bonding.  The  remainder  of  the  sized 
and  unsized  fibers  together  with  the  other  finish  variants  did 
not  show  meaningful  difference. 

An  underlying  reason  for  the  lack  of  additional  distinctions 
is  the  large  variation  in  the  shear  strengths  and  the  1:1  rela¬ 
tionship  between  extrapolated  tensile  strength  and  shear  strength. 
This  effect  is  seen  in  the  AS-1  data  shown  in  Appendix  F  whereby 
an  abnormally  high  extrapolated  tensile  strength  translates 
directly  to  a  proportionally  high  shear  strength. 

In  conclusion,  this  test  can  be  used  to  identify  strong  and 
weak  interfacial  bond  strengths  successfully.  However,  due  to 
the  nature  of  this  test,  it  cannot  reliably  differentiate  small 
differences . 
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TABLE  XVI.  SINGLE  FILAMENT  INTERFACIAL  SHEAR  STRENGTH 
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38.  HISTOGRAPH  OF  CRITICAL  ASPECT  RATIO  OF 
T-6300  STANDARD  SIZED  CARBON  FIBERS 
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39.  HISTOGRAPH  OF  THE  CRITICAL  ASPECT  RATIO  OF 
AS-1  UNSIZED  CARBON  FIBER 
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FIGURE  40.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  RATIO  OF  CELION® 
6000  WITH  NO  SURFACE  TREATMENT  AND  UNSIZED 
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FIGURE  41.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  RATIO  OF  CELION® 
6000  WITH  SURFACE  TREATMENT  BUT  UNSIZED 
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FIGURE  42.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  RATIO  OF  CELION® 
6000  STANDARD  EPOXY  COMPATIBLE  SIZE 
(COMMERCIAL  PRODUCT) 


FIGURE  43.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  RATIO  OF  CELION® 
6000  WITH  ETP-50  SIZE 
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FIGURE  44.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  RATIO  OF 
CEL I ON®  6000  WITH  ETP-10  SIZE 
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FIGURE  45.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  OF 
CELION®  6000  WITH  FRE-25  SIZE 
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FIGURE  46.  HISTOGRAM  OF  THE  CRITICAL  ASPECT  RATIO  OF 
CELION®  6000  WITH  SILICONE  RUBBER  SIZE 
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Although  this  test  cannot  differentiate  small  differences 
numerically  in  interfacial  bond  strength,  it  does  provide  some 
information  on  fracture  behavior  of  different  fiber  types/finish 
variants.  Three  distinctly  different  types  of  fracture  behaviors 
were  observed.  The  optical  photomicrographs  of  these  fracture 
behavior  patterns  are  shown  in  Figures  47  through  49  and  a 
schematic  diagram  of  these  three  fracture  modes  is  presented  in 
Figure  50 . 

From  this  information,  it  is  clear  that  the  commercial  epoxy- 
based  sized  fibers  provide  a  higher  interfacial  shear  strength  and 
that  brittle  matrix  cracks  were  initiated  at  the  fiber  fracture 
ends.  However,  ETP-50  sized  Celion®  6000  produced  a  ductile, 
energy  absorbing  type  of  crack  at  the  interphase.  Further  in¬ 
crease  in  stress  on  the  latter  fibers  created  a  sharp,  brittle 
crack  extending  into  the  bulk  matrix  material.  These  observa¬ 
tions  suggest  that  it  may  be  possible  to  achieve  both  high  shear 
strength  and  increased  fracture  toughness  simultaneously. 
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CROSS  POLAR  TRANSMISSION  PHOTOMICROGRAPH  100X 


TRANSMISSION  PHOTOMICROGRAPH  320X 


FIGURE  47.  FRACTURE  MODE  AT  FIBER  END  OF  LOW  INTERFACIAL 
BOND  STRENGTH  (SILICONE  RUBBER  SIZED  AND 
UNSIZED  CARBON  FIBERS) 


TRANSMISSION  PHOTOMICROGRAPH  3 2 OX 


CROSS  POLAR  PHOTOMICROGRAPH  20 OX 


FIGURE  48.  FRACTURE  MODE  AT  FIBER  END  OF  HIGH  INTERFACIAL  BOND 
STRENGTH  WITH  STANDARD  EPOXY  SIZED  CARBON  FIBER. 
BRITTLE  FAILURE  INTO  MATRIX  AT  FIBER  ENDS. 


CROSS  POLAR  PHOTOMICROGRAPH  320X 


TRANSMISSION  PHOTOMICROGRAPH  3 2 OX 


FIGURE  49.  FRACTURE  MODE  AT  FIBER  END  OF  LOWER  MODULUS  FINISH 
(ETP  50)  CARBON  FIBER.  A  DUCTILE  FRACTURE  MODE  AT 
INTERPHASE  INDICATES  A  TOUGH  FRACTURE  BEHAVIOR. 

THE  INNERLAYER  IS  IDENTIFIED  BY  THE  SHARP  BRITTLE 
CRACK  EXTENDED  INTO  THE  BULK  MATRIX  MATERIAL. 
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a)  Typical  unsized  fiber  or  weak  interfacial 
bonding  -  no  matrix  crack. 


b)  Celion®6000  standard  epoxy  sized  fiber 
brittle  cracks  initiated  at  fiber  ends. 


different  crack  geometry  immediately  adjacent 
to  fiber.  Brittle  matrix  crack  developed 
outside  the  interphase  region. 


FIGURE  50.  A  SCHEMATIC  DIAGRAM  OF  THE  THREE  TYPES  OF  FAILURE 
MECHANISMS  AT  FIBER  ENDS  WITH  DIFFERENT  FINISHES 
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Effect  of  Sizings  on  Composite  Properties 


The  finish  variants  were  formulated  and  applied  to  Celion®  6000 
unsized  carbon  fiber.  A  three-inch  wide  prepreg  tape  with  a  thick¬ 
ness  of  0.127  mm  (5  mils)  was  prepared  from  each  variant.  Laminates 
were  then  fabricated  using  a  match-metal-die  molding  technique. 
0°-tensile  properties  together  with  0°-flexural,  90°-flexural  and 
short  beam  shear  (SBS)  at  both  room  temperature  and  93.3°C  (200°F) 
were  measured.  An  ETI-Dynatup  drop  dart  instrumented  impact  tester 
(Appendix  J)  was  used  to  measure  the  impact  toughness  of  additional 
laminates.  An  approximately  quasi-isotropic ,  [+45/90/-45/0/  45/90]  s 
lay-up  was  used  for  the  impact  toughness  measurements.  To  character¬ 
ize  the  toughess  of  these  laminates,  the  maximum  penetration  force, 

Pp  recorded  on  the  tester  is  reported  along  with  Pj ,  the  force  assoc¬ 
iated  with  initial  damage.  These  composite  properties  together  with 
the  wetting  force  measurements  and  the  single  filament  adhesion  strength 
are  presented  in  Table  XVII. 

To  further  evaluate  this  toughness  enhancement  obtained  with 
resin  system  "F,"  a  commercial  resin  system  520S  is  used.  Finish 
variant  ETP-50  and  a  standard  commercial  size  (serving  as  a  control) 
were  applied  to  Celion®  6000  unsized  fiber  and  their  properties  are 
presented  in  Table  XVIII. 

Comparing  the  wetting  force  measurements  with  both  the  single 
filament  adhesion  strength  and  the  composite  shear  strength,  it  can 
be  seen  that  good  wetting  is  required  for  good  adhesion.  However, 
better  wetting  as  measured  herein  does  not  necessarily  provide  stronger 
interfacial  bonding.  This  also  suggests  that  once  the  fibers  are 
"wetted,"  regardless  of  how  the  wetting  is  achieved,  either  forcibly 
wet  or  due  to  differences  in  surface  energetics,  then  some  adhesion 
can  be  achieved.  This  is  identified  in  the  silicone  rubber  sized  car¬ 
bon  fiber  composites  where  structural  integrity  is  maintained  in  spite 
of  chemical  incompatibility.  However,  stronger  interfacial  bond 
strengths  are  clearly  obtained  with  the  other  finish  variants  and  the 
control.  This  probably  arises  from  more  favorable  chemical  bonding 
at  the  interface. 

Although  a  1:1  correlation  between  the  single  filament  adhesion 
strength  and  the  composite  shear  strength  was  not  demonstrated,  the 
single  filament  test  is  still  very  useful.  This  test  method  provides 
valuable  insight  into  adhesion  fundamentals,  the  underlying  composite 
shear  strength  and  the  failure  mechanisms  involved  in  corresponding 
composites . 

It  is  noteworthy  that  three  of  the  finish  variants  selected  in 
this  study,  ETP-50,  ETP-10  and  FRE-25  all  show  significant  increases 
in  impact  toughness  with  no  corresponding  decreases  in  other  mechani¬ 
cal  properties  either  at  room  or  elevated  temperature.  This  impact 
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TABLE  XVIII.  CEL I ON®  6000/5208  LAMINATE,  PROPERTIES 


TEST 

SAMPLE 

PROPERTY 

CONDITION 

TEMPERATURE 

C-6KE 

C-6K/EPT50 

Impact  ^ 

PF  (N) 

RT 

1348 

1744 

et(j) 

7.46 

9.53 

(2) 

Tensile 

STR (MPa) 
MOD (GPa) 

RT 

1669 

140.7 

1834 

142.1 

STR (MPa) 
MOD (GPa) 

132  °C 

1738 

154.5 

1814 

153.1 

STR (MPa) 
MOD (GPa) 

WET 

132  °C 

1745 

156.6 

1786 

153.1 

0°  Flex1 (2) 

STR (MPa) 
MOD (GPa) 

RT 

1972 

118.6 

2041 

121.4 

STR (MPa) 
MOD (GPa) 

132°C 

1731 

119.3 

1986 

119.3 

STR (MPa) 

WET  (3^ 

132°C 

1076 

869 

90°C  Flex 

STR (MPa) 
MOD (GPa) 

RT 

87.6 

9.31 

82.8 

9.03 

STR (MPa) 

132°C 

63.5 

56.6 

STR (MPa) 

WET  ^3) 

132°C 

31.0 

26.9 

SBS 

STR (MPa) 
STR (MPa) 

RT 

132°C 

117.2 

79.3 

111.7 

80.7 

STR (MPa) 

WET  ^ 

132°C 

53.8 

49.0 

±45°  Tension 

STR  (MPa) 
MOD (GPa) 

RT 

161.4 

4.76 

159.3 

4  .35 

STR (MPa) 
MOD (GPa) 

132  °C 

144.8 

3.72 

131.7 

3.24 

STR (MPa) 
MOD (GPa) 

WET  (3) 

132  °C 

129.7 

2.69 

112.4 

2.28 

(1)  Normalized  to  1.5  xlO-3  meter  thickness. 

(2)  Normalized  to  62%  vol.  fraction. 

(3)  Immersed  in  water  at  71.1°C  for  14  days. 

Note:  All  data  except  impact  are  average  of  5  specimens. 

An  average  of  6  specimens  is  used  for  impact  data. 


110 


toughness  improvement  is  also  shown  in  the  5208  resin  system.  The 
mechanisms  involved  are  not  fully  understood.  It  is  postulated  that 
the  softer  innerlayer  at  the  interphase  diminishes  some  of  the  in¬ 
ternal  residual  stresses  and  thereby  eliminates  stress  concentrations. 
The  significant  improvements  in  impact  toughness  may  arise  from  damp¬ 
ing  effects  of  shock  waves  at  each  fiber  interphase  where  an  energy 
dissipating  material  is  not  present.  In  any  case,  this  study  has 
demonstrated  that  tailoring  of  the  interphase  properties  can  enable 
significant  improvement  of  composite  performance. 


MODEL  DEVELOPMENT 


The  analytical  model  discussed  under  the  section,  "Stress 
Transfer  from  Matrix  to  Fiber, "  is  valid  only  if  there  is  no 
interphase  region.  However,  it  was  shown  by  Drzal  et  al^1)  and 
by  the  authors  that  an  interphase  does  exist  and  its  effect  on 
composite  properties  is  significant.  Therefore,  a  modification 
of  the  existing  analytical  model  is  essential. 

From  the  observations  and  the  results  of  this  study,  several 
recommendations  can  be  made  by  the  authors  to  support  the  develop¬ 
ment  of  a  new  analytical  model  which  includes  a  third  phase  -  the 
interphase  region. 

I 

Experimentally,  the  authors  have  shown  that  incorporation  of 
a  lower  modulus  interphase  region  results  in  the  observation  of 
a  different  fracture  mechanism.  This  energy  absorbent  inner  layer 
also  improves  the  impact  toughness  of  the  composite.  Further, 
there  is  some  indication  that  the  transverse  matrix  crack  which 
occurs  at  a  fiber  fracture  is  blunted  and  that  the  tensile  strength 
of  a  laminate  is  slightly  enhanced.  However,  it  has  to  be  made 
clear  that  not  only  the  physical  properties  of  the  interphase  are 
important,  but  the  chemical  nature  of  the  interphase  also  has  to 
be  considered.  This  is  because  some  chemical  bonding  is  required 
to  provide  strong  interfacial  shear  strength  to  maximize  the  load 
transfer  across  the  interface.  The  thickness  of  this  interphase 
may  itself  be  crucial.  An  annulor  layer  which  is  too  thin  may 
not  be  effective  in  improving  impact  resistance,  while  too  thick 
a  layer  may  lower  other  laminate  properties  such  as  transverse 
tension  and  compression,  etc. 

A  schematic  diagram  of  the  modified  model  with  the  incorporated 
inner  layers  is  shown  in  Figure  51.  A  micromechanistic  approach 
can  be  taken  to  analyze  the  stresses  in  each  region.  Two  approaches 
with  different  assumptions  can  be  taken. 

First,  consider  the  fiber  to  be  very  rigid,  when  under  stress 
and  the  true  deformation  of  the  fiber  to  be  negligible  compared 
to  the  matrix  and  the  inner  layer.  In  this  case,  a  ductile  inner 
layer  of  low  shear  modulus  combines  with  the  matrix  to  undergo 
large  deformations  in  transferring  tensile  load  to  the  fiber. 

Or  second,  consider  the  fiber  to  be  an  inclusion  in  the  inner 
layer  material,  then  treat  the  fiber/inner  layer  composite  with 
its  inherent  characteristics  as  an  inclusion  in  the  bulk  matrix. 

After  the  stress  states  are  determined,  then  classical  micro- 
mechanistic  approach  can  be  used  to  study  the  composite  properties. 
The  resulting  analytical  model  should  include  a  third  phase  intro¬ 
duced  by  either  (or  both)  of  the  approaches  illustrated  in  Figure  51. 

To  incorporate  the  experimental  results  reported  here,  some 
refining  experiments  are  necessary: 
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y  is  the  shear  modulus 
k  is  the  bulk  modulus 
and  subscript 

m  =  matrix 
I  =  interphase 
f  =  fiber 


FIGURE  51.  A  SCHEMATIC  DIAGRAM  OF  AN  INNERLAYERED  INCLUSION 


What  is  the  modulus  of  the  interphase? 

What  is  the  composition  gradient  (and  property  gradient) 
across  the  interphase? 

What  is  the  effect  of  interphase  thickness? 

What  is  the  degree  of  bonding  required  at  the  fiber 
surface?  (Is  the  perfect  bonding  assumption  a  valid  one?) 

What  is  the  effect  of  finish  on  the  load  transfer  to  an 
adjacent  fiber?  (What  is  the  loaded  length  of  the  ad¬ 
jacent  fiber  in  the  vicinity  of  a  fracture?) 

The  answers  to  these  questions  and  additional  experimental 
observations  would  be  amenable  to  analytical  modelling  to  predict 
and  account  for  the  observed  phenomenon.  A  valid  model  must  be 
consistent  with  the  observations. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


An  extensive  literature  review  on  the  role  of  interface  and 
interphase  in  composite  properties  has  been  completed  and  ex¬ 
periments  were  designed  to  clarify  some  of  the  uncertainties 
encountered  in  this  survey. 

Wetting  angle  measurement  on  fiber  surfaces  is  the  most  commonly 
used  technique  to  predict  adhesion  strength  due  to  wetting.  However, 
no  direct  correlation  has  been  shown  between  wetting  angle  measure¬ 
ments  and  fiber- to-matrix  bond  strength.  In  this  study,  wetting 
behavior  has  been  characterized  by  the  Wilhelmy  wetting  force 
measurement  and  interfacial  shear  strengths  have  been  studied  via 
the  single  filament  critical  length  determination.  The  results 
indicate  that  intimate  contact  is  necessary  for  adhesion,  but 
chemical  bonding  may  be  responsible  for  higher  interfacial  bond 
strength . 

Further,  the  results  also  indicate  that  by  tailoring  the  inter¬ 
phase,  different  failure  mechanisms  can  be  achieved  without  signif¬ 
icant  loss  in  interfacial  bond  strength.  This  approach  was  in¬ 
vestigated  more  thoroughly  by  an  in  depth  evaluation  of  laminate 
properties.  A  significant  improvement  in  impact  toughness  was 
obtained  without  loss  of  static  mechanical  properties.  The  con¬ 
clusions  of  this  study  can  be  summarized  as  the  following: 

1.  Good  wetting  properties  are  essential  to  provide  good  bond- 
ability,  but  better  wetting  properties  do  not  necessarily 
provide  stronger  adhesive  joint  strengths. 

2.  Single  filament  critical  length  determination  is  a  good 
analytical  tool  to  characterize  carbon  fiber/epoxy  adhesion 
properties  and  inherent  failure  mechanisms. 

3.  Single  filament  adhesion  strength  and  composite  shear  strength 
are  only  loosely  correlated. 

4.  Different  failure  mechanisms  are  identified  in  testing  of 
single  filaments  of  carbon  fiber  with  different  finish  variants 
embedded  in  an  epoxy  matrix. 

5.  A  significant  30%  increase  in  impact  toughness  was  measured 
without  loss  at  room,  elevated  temperatures  and  hot -wet 
mechanical  properties  by  incorporation  of  a  "tailored- 
interphase"  into  the  composite. 
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6.  A  modified  analytical  model  is  necessary  to  better  understand 
the  effect  of  the  interface  and  interphase  on  composite  per¬ 
formance  . 

7 .  A  better  understanding  of  the  load  transfer  mechanisms  in¬ 
volved  in  this  more  complex  3-component  composite  environ¬ 
ment  is  necessary  to  prescribe  the  optimum  interphase  proper¬ 
ties  . 


APPENDIX  A.  CALCULATIONS  FOR  THE  THICKNESS  OF  A  1.5%  BY 
WEIGHT  EPOXY  SIZE  ON  CARBON  FIBER 


Density  of  carbon  fiber,  Pf  =  1.76  g/cc 
Diameter  of  carbon  fiber,  r  =  7  x  10“4  cm 
Density  of  epoxy  size,  Pg  =  1.182  g/cc 
Assume  unit  length,  SL  =  1  cm 

Volume  of  the  fiber,  V^  =  irr^fl,  =  1.539  x  10  8  cc 
Weight  of  the  fiber,  Wf  =  VfPf  =  2.709  x  10"6  g 
.*.  Weight  of  the  eooxy  size,  Ws  =  1.5%  Wf 

=  4.064  x  10"8  g 


Volume  of 

But,  Vs  = 
R  = 

Thickness 


the  epoxy  size,  Vs  =  Ws/P s 

=  3.438  x  10“8 
Tr(R2-r^)£  =  3.438  x  10-8  Cc 
7.08  x  10-4  cm 

o 

of  eooxy  size,  A  =  R-r  =  800  A 


cc 
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APPENDIX  B.  APPARATUS  FOR  WETTING  FORCE  MEASUREMENT 


A  Cahn  26  automatic  electrobalance  was  used  to  measure  the 
wetting  force  of  a  single  filament  carbon  fiber.  This  electro¬ 
balance  was  interfaced  through  a  BCD  interface  module  to  an 
HP-85A  microcomputer.  This  arrangement  allowed  the  data  to  be 
stored,  plotted  and  manipulated  during  the  experiment. 

In  order  to  assess  the  advancing,  receding  and  equilibrium  wetting 
force,  a  constant  speed  moving  stage  for  the  wetting  liquid . 
was  required.  This  stage  was  made  from  a  modified  micromanipulator 
and  a  high  torque,  variable  speed  electric  motor.  The  electric 
motor  was  controlled  automatically  by  the  microcomputer  to  allow 
automation  of  this  experiment.  The  program  developed  for  this 
process  is  listed  in  Appendix  C  and  a  schematic  diagram  of  this 
set-up  is  presented  in  Figure  B-l. 

A  single  filament  of  carbon  fiber  was  selected  and  attached  to. 
the  end  of  a  thin  gauge  hang-down  wire  with  rubber . cement .  This 
fiber  was  then  connected  to  the  balance  and  the  weight  of  the 
fiber  was  tared.  The  wetting  liquid  contained  in  a  perfectly, 
cleaned  beaker  was  raised  by' the  micromanipulator  until  the  fiber 
end  was  just  immersed  in  the  liquid.  Then  the  program  would 
automatically  accumulate  the  data  from  the  microbalance . and  control 
the  stage  to  obtain  the  wetting  force  during  the  advancing,  station¬ 
ary  and  receding  positions. 
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Closed  Chamber  to 
Eliminate  Mr  Current  and  Dust 


FIGURE  B-l.  APPARATUS  ASSEMBLY  FOR  WETTING  FORCE  MEASUREMENT 


APPENDIX  C.  WETTING  FORCE  MEASUREMENT  PROGRAM 


10  REM  - >AUT0ST 

15  OPTION  BASE  1 
40  ON  ERROR  GOTO  80 
50  LOADBIN  ’REDZER’ 

60  OFF  ERROR 
70  CHAIN  •WETBAL* 

80  IF  ERRN=25  THEN  100 

90  PRINT  ’Error  number  ’ t ERRN » ’  was  encountered  at  line  ’ ? ERRL » ’  .  ’  @  GOT 
0  120 

100  PRINT  ’A  binary  program  is  already  in  memory  causing  an  error.  Plea 

se*  a  ~  , 

105  PRINT  'eliminate  this  binary  (i.e.  by  turning  machine  off  or  execut 

ins' 

110  PRINT  * 'SCRATCH' )♦  Then  you  may  beSin  the  program  again.’ 

120  END 
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Introduction 


REDZER  is  a  binary  program  designed  to  enhance  array  manipulation  capa¬ 
bilities  of  the  HP-85.  An  array  can  be  reorganized  (redimensioned)  or 
initialized  to  zero  with  use  of  one  of  the  two  statements  offered.  This 
binary  program  can  be  loaded  into  the  HP-85  by  the  following  command: 

LOADBIN  "REDZER" 

For  a  further  discussion  of  binary  programs,  see  page  193  of  the  HP-85 
Owner's  Manual  and  Programing  Guide. 

ARRAY  (RE  ORGANIZATION 

You  can  reorganize  an  array  into  a  more  useful  configuration  by  redim- 
sioning  that  array.  This  changes  the  working  size  of  the  array;  subse¬ 
quent  statements  affect  only  the  elements  included  in  the  new  working 
size.  However,  elements  not  included  in  the  new  working  size  are  still 
associated  with  the  array.  The  values  of  these  elements  are  not  changed, 
and  they  can  be  accessed  if  the  array  is  again  redimensioned. 


REDIM  array  (redim  subscripts)  [.array  (redim  subscripts)...] 


The  redimensioning  subscripts  are  numeric  expressions,  variables,  or  con¬ 
stants  that  specify  a  new  upper  bound  for  each  dimension.  The  number 
of  subscripts  must  be  the  same  as  the  number  specified  in  the  original 
DIM,  REAL,  SHORT,  or  INTEGER  statements.  Furthermore,  the  total  number 
of  elements  in  the  new  working  size  cannot  exceed  the  number  originally 
dimensioned. 
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Examples: 

DIM  A(4) ,B(2,4),C(1,9)  OPTION  BASE  0  assumed. 

REDIM  A(3)  Redimensions  working  size  from 

five  to  four  elements. 

REDIM  B (3 ,2 )  Redimensions  B  from  3x5  matrix 

(15  elements)  into  4x3  matrix 
(12  elements). 

REDIM  A(4),B(2,4)  Redimensions  A  and  B  back  to 

original  sizes. 

X=2  @  REDIM  C(2*X-1,10/X-1)  Redimensions  C  from  2  x  10  matrix 

into  4x5  matrix. 

When  a  matrix  is  redimensioned,  the  values  of  its  elements  are  reassigned 
to  different  positions  within  the  matrix.  Values  of  matrix  elements  are 
stored  in  order  from  left  to  right  along  each  row,  from  the  first  row 
to  the  last.  The  redimensioning  takes  the  elements  out  of  the  matrix 
in  that  order,  and  reassigns  them  in  accordance  with  the  new  working 
size  of  the  matrix. 

The  following  example  shows  how  values  of  matrix  elements  are  reassigned 
when  a  matrix  originally  declared  to  be  3  x  3  is  redimensioned  into  a 
2x2  matrix.  The  values  of  the  original  matrix  elements  are  integers 
that  indicate  the  order  in  which  the  elements  are  stored  before  the  re¬ 
dimensioning. 


Example: 

10  OPTION  BASE  1 
20  DIM  A(3,3) 

30  DATA  1,2, 3, 4, 5, 6, 7 ,8, 9 
40  FOR  1-1  TO  3 
50  READ  A(1 ,1)  ,A(1 ,2)  ,A(1 ,3) 
60  NEXT  I 

70  K=3  @  G0SUB  130 
80  REDIM  A (2 ,2) 


Array  A  is  originally  dimensioned  as 
3x3. 


Assigns  values  1  through  9  to  the 
elements  of  matrix  A. 


Displays  original  3x3  matrix. 

Redimensions  A  down  to  a  2  x  2 
matrix. 


03 


90  K=2  0  GOSUB  130 
100  REDIM  A (3 , 3) 

110  K=3  @  GOSUB  130 
120  END 

130  FOR  1=1  TO  X 
140  FOR  J=1  TO  X 
150DISP  A (I ,J); 
160  NEXT  J 
170  DISP 
180  NEXT  I 
190  DISP 
200  RETURN 


Displays  redimensioned  2x2 
matrix. 

Redimensions  A  back  up  to  a 
3x3  matrix. 

Displays  redimensioned  3x3 
matrix. 


Subroutine  that  displays  matrix  A. 


um 

1  2  3 
4  5  6 
7  8  9 


Nine  elements  of  original  3x3 
matrix. 


1  2 

3  4 

12  3 

4  5  6 
7  8  9 


Four  elements  of  redimensioned 
2x2  matrix.  Values  of  elements 
have  been  reassigned  sequentially 
within  new  working  size. 

Nine  elements  of  3  x  3  matrix  with 
original  values  still  assigned. 


Note  that  redimensioning  a  matrix  does  not  isolate  a  submatrix.  In  other 
words,  if  you  redimension  a  3  x  3  matrix  into  a  2  x  2  matrix,  the  resulting 
matrix  is  not  the  2x2  submatrix  from  the  upper  left  corner  of  the  original 
matrix. 
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5  !  WILHELMY  BALANCE  DATA  ACQUISITION . . .PROGRAM  WETBAL. ♦ .G.V.  NELSON  5/7 
/81  ♦  .  ♦ VERSION  5.12.81 

8  OPTION  BASE  1 

9  DIM  D(l»1400) tS2(10) »S1$C603»V1*C84] 


WILHELMY  BALANCE*  0  DISP  *  D 

0  INPUT  Y$ 

data  collection  interval  (1~3  seconds 


collected  for  Advancing  Region  (300  m 


collected  for  Still  Region  (200  max)* 


27  INPUT  N2 

28  CLEAR  0  ON  KEY*  1>*START*  GOTO  50 

29  ON  KEY*  6 »  *  UP  *  GOSUB  6200 

30  ON  KEY*  7f*D0WN*  GOSUB  6250 

31  ON  KEY*  8  r  *  STOP  *  GOSUB  6300 

32  OFF  KEY*  3  0  OFF  KEY*  4  0  OFF  KEY*  2 

37  ! 

38  GCLEAR  0  CLEAR  0  KEY  LABEL 

39  GOSUB  6000  !  INITIALIZE  INTERFACE 

40  GOTO  40 

50  DISP  'Enter  Sample  Description  (20  Chars  maxf  no  commas)* 

51  INPUT  A$ 

52  PRINT  USING  *7/*  »  @  PRINT  *  WILHELMY  BALANCE*  0  PRINT 

DATA  ACQUISITION*  0  PRINT 

53  PRINT  A$ 

54  PRINT  0  PRINT  'DATA  INTERVAL*  *  J  S2»  *  SECONDS/POINT* 

55  PRINT  'ADVANCING  REGION** ?N1 i *  POINTS* 

56  PRINT  'STILL  REGION*  *  »N2 »  *  POINTS* 

57  PRINT  'RECEEDING  REGION**  »N1 »  *  POINTS* 

58  PRINT  0  PRINT 

80  Y 1 =0  0  Y2=0  0  M=0 

85  GOSUB  7000  0  D1=D/100  0  D2=D1  !  START  TIMER 

86  LOCATE  OtlZOrOtlOO 

87  FRAME 

90  SCALE  0r2*Nl+N2»0.Dl 
100  MOVE  NlrO  0  DRAW  N1.D1 

102  MOVE  N1 +N2 1 0  0  DRAW  N1+N2»D1 

103  LORG  6  0  MOVE  N1/2.D1  0  LABEL  “ADVANCING* 

104  MOVE  N1 +N2/2  »D1  0  LABEL  'STILL* 

105  MOVE  N1+N2+N1/2»D1  0  LABEL  ’RECEDING* 

120  MOVE  0.0 

135  GOSUB  6230  !  MOTOR  ON-FWD 

140  1*0  0  J=I 

150  IF  I=J  THEN  150 

155  GOSUB  5030 

160  Y(I) » Y=D 

170  J=I 

180  GOSUB  770 

182  PEN  1 

185  MOVE  I » Y ( I ) 


10  DIM  Y(700) fM(3) »S(3) f A$C203 

15  1=0  0  02=1400 

16  DISP  USING  *7/*  i  0  DISP  * 

ATA  ACQUISITION* 

17  DISP  'Data  from  tape  file(Y/N)* 

18  IF  UPC*(Y$Cldl)  =  *Y*  THEN  2000 

21  DISP  USING  *3/*  i  0  DISP  'Enter 
/point ) " 

22  INPUT  S2 

23  DISP  'Enter  no.  of  points  to  be 
ax)  * 

24  INPUT  N1 

26  DISP  'Enter  no.  of  points  to  be 
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190  PLOT  IrY(I) 

200  IF  I =N 1  THEN  GOTO  220 
210  GOTO  150 
220  ! 

225  GOSUB  6330  !  MOTOR  OFF 
230  PRINT 

250  OFF  TIMER*  1  0  PRINT  'ADVANCING! * 

254  PRINT  USING  ■ 10A > 6D ♦ 3D t / t 10A > 6D ♦ 3D '  r  '  MEAN  =  '»M»'  STDEV  =  *?SQ 
R < SI ) 

260  PRINT 

270  M  < 1 ) =M  0  M=0  0  Y1  =  0  0  Y2=0 

280  !  TURN  ON  TIMER  FOR  STILL  REGION 

290  GOSUB  7000  !  START  TIMER 

310  1=0  0  J=I 

320  IF  I=J  THEN  320 

325  GOSUB  5030 

330  Y<N1+I)»Y=D  0  J=I 

340  GOSUB  770 

342  MOVE  N1+I»Y(N1+I) 

344  DRAW  Nl+IrY(Nl+I> 

346  IF  I=N2  THEN  390 
350  GOTO  320 

390  OFF  TIMER*  1  0  PRINT  'STILL! ' 

400  PRINT  USING  ' 10A > 6D . 3D » / *  1 OA » 6D *  3D '  i  *  MEAN  =  '»M»*  STDEV  =  '»SQ 
R  (  S 1  ) 

410  PRINT 

420  M ( 2 ) =M  0  M=0  0  Y1=0  0  Y2=0 

430  !  TURNON  TIMER  FOR  RECEDING  REGION 

450  GOSUB  7000  !  START  TIMER 

455  GOSUB  6280  !  MOTOR  REV 

460  1=0  0  J=I 

470  IF  I=J  THEN  470 

475  GOSUB  5030 

480  Y (N1+N2+I ) ? Y=D  0  J=I 

490  GOSUB  770 

492  MOVE  N1+N2+I»Y(N1+N2+I) 

494  DRAW  N1  +  N2  +  I ,Y(N1+N2  +  I ) 

496  IF  I=N1  THEN  510 

500  GOTO  470 

510  !  STOP  ROUTINE 

520  GOSUB  6300  !  STOP  MOTOR 

530  OFF  TIMER*  1  0  PRINT  'RECEDING  REGION!' 

540  PRINT  USING  ' 10A » 6D . 3D  * / , 10A » 6D . 3D '  *  '  MEAN  =  '»MJ'  STDEV  =  '»SQ 
R  (  S 1 ) 

550  PRINT  0  COPY 

554  ON  KEY*  5> 'STORE*  GOTO  1000 

555  ON  KEY*  3, 'LIST'  GOTO  900 

556  ON  KEY*  4>'PL0T'  GOTO  560 

557  GOTO  37 

558  ON  KEY*  2 » ' READ'  GOTO  1900 

559  GOTO  37 

560  !  REPLOT  GRAPH  TO  SCALE 

561  GCLEAR 

562  DISP  'WHICH  REGION  DO  YOU  WANT  TO  PLOT  < A=ADVANCINGJ  S=STILLi  R=RECE 
DING)  ' 

563  INPUT  A1$0  N3=0 

564  IF  UPC$(A1*E1»13)='A'  THEN  N3=l  0  N4=N1  0  M$= ■ ADVANCING ' 

565  IF  UPC$(A1$C1?1])='S*  THEN  N3=N1+1  0  N4=N1+N2  0  M$='STILL* 


C-6 


566  IF  UPC$(A1$C1»13)=*R'  THEN  N3=N1+N2+1  0  N4=2*N1+N2  0  M$=  *  RECEDING  * 

567  IF  N3=0  THEN  562 
570  PRINT  USING  *  4/ *  f 
580  Y8  =  -l ♦ E99  0  Y9=1*E99 
590  FOR  I=N3  TO  N4 

600  Y8  =  MAX  <  Y  < I ) f  Y8 )  0  Y9=MIN ( Y ( I ) »  Y9 ) 

610  NEXT  I 

611  CLEAR  0  DISP  'MIN  Y  VALUE='iY9  0  DISP  'MAX  Y  VALUE='iY8  e  DISP  'DO  Y 
OU  WANT  TO  SET  SCALE< Y/N) *  i 

612  INPUT  A1$0  IF  UF‘C*(Al$i:i»i:i>  =  'Y'  THEN  613  ELSE  614 

613  DISP  'ENTER  Ymin>Ymax'  0  INPUT  Y9»Y80  GOTO  615 

614  Y9=0  0  Y8=D1 

615  LOCATE  30 » 120 r 20 > 80 

616  SCALE  N3-1 » N4  » Y9  r  Y8 

617  FXD  0,3  0  LAXES  - ( N4-N3+1 ) /10  » <  Y8-Y9 >710 » N3-1 » Y9  » 2 , 2 

618  SETGU  0  MOVE  70»3  0  LORG  4  0  LABEL  'TIME(see)' 

619  HOVE  8 » 45  0  LDIR  90  0  LABEL  ' FORCE < md ) '  0  LDIR  0  0  SETUU 
635  MOVE  <N3+N4)/2»Y8+< Y8-Y9>*3/20  0  LORG  4  0  LABEL  M$ 

640  MOVE  N3- 1 t 0 
650  FOR  I=N3  TO  N4 
655  MOVE  I » Y ( I ) 

660  LORG  5 

665  IF  N4-N3>40  THEN  PLOT  I*Y(I)  ELSE  LABEL  '+' 

670  NEXT  I 
680  COPY 

690  PRINT  USING  '7/'  i 
715  GOTO  37 
720  END 

770  !  RUNNING  AVERAGE  t  STD  DEV 

780  M=( <1-1 )*M+Y)/I 

790  Y1=Y1+Y 

800  Y2=Y2+Y*Y 

810  S1=<Y2-Y1*Y1/I)/I 

820  RETURN 

900  !  LIST 

902  CLEAR 

904  DISP  'WHICH  REGION  DO  YOU  WANT  TO  LIST  < A=ADVANCING J  S=STILL»  R=RECE 
DING) ' 

906  INPUT  A1*0  N3=0 

908  IF  UPC$(A1$C1»13)='A'  THEN  N3=l  0  N4=N1  0  M*= ' ADVANCI NG ' 

910  IF  UF’Cf  <A1$C1»11)  =  'S'  THEN  N3=N1  +  1  0  N4=N1+N2  0  M$='STILL' 

912  IF  UPC$<Al$Clf 13)='R'  THEN  N3=N1+N2+1  0  N4=2*N1+N2  0  M$='RECEDING' 

914  IF  N3=0  THEN  904 

915  PRINT  USING  '3/'  i  0  PRINT  M$  0  PRINT  'TIME(sec)  FORCE < m3) ' 

920  DISP  'ENTER  NUMBER  OF  DATA  POINTS  PER  PRINTED  VALUE'  0  INPUT  N5 
950  FOR  I=N3  TO  N4  STEP  N5 

952  F9=3 

954  IF  D2>20  THEN  F9=2 
956  IF  D2<=2  THEN  F9=4 

958  ON  F9  GOTO  959 » 959 » 960 » 965 

959  PRINT  USING  * DDDZ . DD r 4X r DDDZ « DD '  »  I*S2»Y<I)  0  GOTO  970 

960  PRINT  USING  ' DDDZ ♦ DD  r  4Xf DDDZ  *  DDD  *  J  I*S2»Y<I)  0  GOTO  970 
965  PRINT  USING  ' DDDZ ♦ DD » 4X » DDDZ ♦ DDDD '  i  I*S2»Y<I) 

970  NEXT  I 

980  PRINT  USING  '7/'  r 
990  GOTO  37 

1000  !  STORE  ACCORDING  TO  BSDM  FORMAT 
1010  01=2*N1+N2  0  N7=2  0  S7=3 


1020  U1$=*TIME  FORCE  1 

1030  S1*=*ADV«  STILL  RECED." 

1040  S2 ( 1 ) =1  8  S2( 2 ) =N1 +1  8  S2 ( 3 ) =N1+N2+1 
1045  REDIM  D ( N7  >  01 ) 

1050  FOR  1=1  TO  01 
1060  D ( 1 f I ) =S2*I 
1070  D(2rI)=Y(I) 

1080  NEXT  I 
1140  CLEAR 

1150  DISP  1  *  *  #  *  STORE  DATA  #  *  *  *  *8  PRINT 

1160  DISP  'File  names  must  be  <=  6  charac-  ters  Iona*  (Type  to  exit) 

■ 

1170  DISP  "Name  of  data  file  =  *58  INPUT  F$ 

1180  IF  LEN ( F$ ) >0  AND  LEN(F$)<7  THEN  1210 
1185  IF  F$= " E "  OR  NOT  LEN(F$)  THEN  1530 
1190  GOSUB  1590 
1200  GOTO  1140 

1205  IF  F$C1 » 13="E*  AND  LEN(F$)=1  THEN  1530 

1210  DISP  ' Is  data  medium  placed  in  tape  drive(Type  'E'  to  exit)*J8  IN 
PUT  N$ 

1220  ON  FNA ( N$ >  GOTO  1210 » 1230  *  1210 r 1530 
1230  CLEAR 

1240  ON  ERROR  GOSUB  1710 
1250  X=0 

1260  ASSIGN*  1  TO  F* 

1270  OFF  ERROR 

1280  IF  X=1  THEN  1360 

1290  DISP  "A  file  of  this  name  has  been  found*  Do  you  wish  to  store* 

1300  DISP  "present  data  set  under  this  name  and  destroy  old  data  set 

" ) 8  INPUT  N$ 

1310  ON  FNA ( N$ )  GOTO  1290 r 1370 » 1 1 40 f 1530 

1330  GOSUB  1590 

1340  CLEAR  8  GOTO  1290 

1360  CREATE  F*, (8+02*8)  DIO  700+2*700 

1370  ASSIGN*  1  TO  F$ 

1375  CLEAR 

1390  PRINT*  1*1  i  A$r01 * N7 * VI* * S7 *  SI* »S2 ( » ) 

1400  READ*  1*2 
1410  PRINT*  1  *  IK*) 

1420  ASSIGN*  1  TO  * 

1430  CLEAR 

1440  PRINTER  IS  2 

1450  PRINT  ********************************** i *  *  *  *TAB(32) *  *  *  * 

1460  PRINT  "*"iTAB(14)J " STORE "> TAB ( 32 ) J ** " 

1470  PRINT  *** 5TAB(32) »*#* } ********************************** 

1480  PRINT 

1490  PRINT  "Data  and  related  information  is  stored  in  "fFS?**" 

1500  PRINT 

1510  DISP  "Is  program  medium  replaced  in  device  *»8  INPUT  N* 

1520  ON  FNA ( N$ )  GOTO  1510* 1530* 1510* 1530 
1530  CLEAR 
1560  GOTO  37 
1590  BEEP 

1600  DISP  'Improper  response — please  try  again* ■ 

1610  WAIT  1500 
1620  RETURN 

1710  IF  ERRN=67  THEN  X=1 


1720  RETURN 
1730  DEF  FNA  ( A$  ) 

1740  IF  LEN(A*>=0  THEN  FNA=1  ELSE  FNA=P0S ( * YNE ’ » UPC* ( A$C 1 * 1 3 ) ) +1 
1750  FN  END 

1900  !  TRAP  FOR  UNINTENDED  LOSS  OF  DATA 
1905  CLEAR 

1910  DISP  ’Reading  date  from  tape  will  overwrite  current  memory* 

1920  DISP  ’Continue  <E  to  exit)'?!?  INPUT  N* 

1930  ON  FNA ( N$ >  GOTO  1920*2000* 554*554 
2000  !  READ  DATA  FROM  TAPE  BSDM  FORMAT 
2005  CLEAR 

2010  DISP  ’Name  of  data  file=’*0  INPUT  F* 

2020  DISP  ’Is  data  medium  placed  in  tape  drive(Type  'E'  to  exit)’*0  IN 

PUT  N* 

2030  ON  FNA ( N$ )  GOTO  2020*2040*2020*2190 
2040  CLEAR 

2050  ON  ERROR  GOTO  2080 
2060  ASSIGN#  1  TO  F$ 

2070  OFF  ERROR  0  GOTO  2090 

2080  DISP  ’File  ’*F$*’  not  found’  0  GOTO  2010 
2090  READ*  1*1  *  A# » 01  * N7 * VI* » S7 * Sl$ * S2 ( * > 

2100  REDIM  D ( N7  *  01 ) 

2110  READ*  1*2 
2120  READ*  1  *  D ( » ) 

2130  ASSIGN*  1  TO  * 

2140  PRINT  ’DATA  ENTERED  FROM  TAPE  FILE  ’>F*  0  PRINT 
2150  N1 =S2 <  2 ) - 1  0  N2  =  S2(3)-S2(2) 

2160  FOR  1=1  TO  01 
2170  Y<I)=D(2*I) 

2180  NEXT  I 

2185  S2=D(1*2)-D(1*1) 

2186  D1=D (2*01/2) 

2190  GOTO  554 
4900  ! 

4910  !  INTERRUPT  PROCESSING 
4920  ! 

5000  ENTER  3  »  Q3»Q4  !  READ  VALUE 

5001  IF  ABS ( Q3 ) <60000  THEN  5000 
5010  1=1+1 

5020  RETURN 

5021  ! 

5022  !  INPUT  DECODING 

5023  ! 

5030  Q1=Q3 

5031  02=04 

5032  IF  01*03  THEN  5030 

5040  IF  INT ( 02/10 ) =1  THEN  5045 
5042  DISP  ’Overran ge’ 

5045  Q2=RMD (02*10) 

5050  IF  Q2>0  AND  Q2<9  THEN  5100 
5060  ! 

5070  02=1 

5075  DISP  ’Exponent  Error* 

5080  GOTO  5190 

5090  ON  02  GOTO  5140*5160*5060*5180*5060*5060*5060*5185 
5100  ON  Q2  GOTO  5140*5160*5060*5180*5060*5060*5060*5185 
5140  02=  *  1 
5150  GOTO  5190 
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5160  02  =  *  01 
5170  GOTO  5190 
5180  Q2=  *  00 1 

5184  GOTO  5190 

5185  02= ♦ 0001 
5190  D=Q1*Q2 
5210  RETURN 
6000  ! 

6010  !  DATA  INITIALIZATION 
6020  ! 


6030 

RESET  3 

6040 

CONTROL 

3*1 

*  0  ! 

NO  INTERRUPTS 

6050 

CONTROL 

3*2 

*  0  ! 

CLEAR  OUTPUT  BITS 

6060 

CONTROL 

3*3 

*  5  ! 

5  DIGITS 

6070 

CONTROL 

3*4 

*  0  ! 

NO  EXPONENT 

6080 

CONTROL 

3*5 

*  2  ! 

2  FUNCTIONS 

6090 

CONTROL 

3*6 

*  0  ! 

NO  DECIMAL 

6100 

CONTROL 

3*7 

*  0  ! 

NO  HANDSHAKE 

6110 

CONTROL 

3*8 

*  0 

6120 

CONTROL 

3*9 

*  15 

6130 

CONTROL 

3*10 

*  0 

6140 

RETURN 

6200 

! 

6210 

!  MOTOR 

ON-FWD 

6220 

! 

6230 

CONTROL 

3,2 

*  1 

6240 

RETURN 

6250 

j 

6260 

!  MOTOR 

ON-REV 

6270 

! 

6280 

CONTROL 

3*2 

*  3 

6290 

RETURN 

6300 

i 

6310 

!  MOTOR 

OFF 

6320 

! 

6330 

CONTROL 

3*2 

*  0 

6340 

RETURN 

7000 

GOSUB  5000  ! 

READ 

START  VALUE 

7005 

GOSUB  5030  ! 

DECODE  VALUE 

7010 

D=49999*Q2  ! 

SET 

MAX  VALUE 

7020  ON  TIMER*  1*1000*S2  GOSUB  5000 
7030  RETURN 
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APPENDIX  D.  PROCEDURES  FOR  CLEAN  GLASSWARE  TO  PREVENT 
CONTAMINATION  OF  WETTING  LIQUIDS 


1.  Glassware  is  placed  in  boiling  1:1  diluted  nitric  acid  for 
five  minutes. 

2.  Then  it  is  rinsed  with  double  distilled  water. 

3.  It  is  next  steamed  with  double  distilled  water  for  thirty 
minutes . 

4.  It  is  then  vacuum  dried  at  150 °C  for  ten  hours. 

5.  Finally,  it  is  individually  wrapped  in  aluminum  foil  and 
stored  in  a  dessicator. 
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APPENDIX  E.  SINGLE  FILAMENT  GRAPHITE  FIBER  TENSILE  TEST 


A  representative  number  (generally  20)  of  single  filaments 
are  randomly  selected  from  a  tow  sample.  The  filaments  are  center- 
line  mounted  on  special  slotted  paper  tabs  as  shown  in  Figure  E-l. 
The  paper  tabs  are  gripped  such  that  the  test  specimen  is  uniaxially 
aligned  in  the  jaws  of  an  Instron.  An  Instron  constant  cross-head 
speed  of  0.02  inch/min  and  a  chart  speed  of  ten  inches/min  are 
employed.  The  procedure  of  mounting  the  graphite  fiber  on  the  tabs 
is  as  follows: 

1.  The  strand  bundle  is  placed  loosely  on  a  suitable  work  surface. 

2.  A  dental  pick  is  used  to  gently  separate  the  filaments  and  a 
suitable  test  specimen  is  selected. 

3.  A  tab  of  pre-punched  slot  is  used  as  the  mounting  medium. _ 

The  single  filament  is  centered  along  the  slot  of  a  mounting 
tab  and  a  small  strip  of  tape  is  used  to  anchor  one  end  to 
this  tab. 

4.  The  filament  is  lightly  stretched  and  the  opposite  end  is 
anchored  in  the  same  manner . 

5.  A  small  amount  of  sealing  wax  is  carefully  placed  on  top  of 
the  filament  at  each  edge  of  the  slot  to  secure  the  filament 
to  the  tab. 
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Cardboard  Tab 
Wax  Bead 

Graphite  Fiber 

Cut  Tab  Here 


Masking  Tape 


figure  E-l.  graphite  fiber  mounted  on  a  tab 


APPENDIX  F.  SINGLE  FILAMENT  TENSILE  STRENGTH 
AT  VARIOUS  GAUGE  LENGTHS 
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FIGURE  F-l .  T-6300  (STANDARD  SIZED  PRODUCT)  LOT  #1141- 


FIGURE  F-2.  AS-1  (STANDARD  UNSIZED  PRODUCT)  LOT  #156- 


FIGURE  F-3 .  C-6KV  (UNSIZED,  NONTREATED)  LOT  #7411 


9500 
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FIGURE  F-4 .  C-6K  (UNSIZED)  LOT  #9621 


FIGURE  F-5 .  C-6KE  (STANDARD  SIZE)  LOT  #0431 


FIGURE  F-7.  C-6K/ETP-10  LOT  #9Y31 


FIGURE  F-8.  C-6K/FRE-25  LOT  #9732 


p>  p»p»  tr  tr  cr  (dp 


APPENDIX  G.  TENSILE  FIXTURE  TO  STUDY  SINGLE  FILAMENT 
ADHESION  PROPERTIES 


c=l. 27x10  2m 

d-2.54xl0_2m 

-2 

e=1.9xl0  m 

f =6 . 4xl0-3m 

-3 

g=7.9xl0  m 
-2 

h=1.3xl0  m 

-4 

i=8xl0  m 


FIGURE  G-l .  TENSILE  FIXTURE  FOR  LEITZ  OPTICAL  MICROSCOPE 
DRAWN  TO  SIZE 


APPENDIX  H.  SINGLE  FILAMENT  INTERFACIAL  SHEAR 
STRENGTH  DETERMINATION  PROGRAM 


The  purpose  of  this  program  is  to  apply  a  Weibull  Distribution 
statistical  analysis  of  the  critical  aspect  ratio  of  the  fibers 
where  Weibull  parameters  as  well  as  the  average  shear  strength  and 
the  variance  of  the  shear  strength  are  determined. 
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10  !*  ***  INTERFACIAL  SHEAR  STRENGTH  PROGRAM  —  CRITICAL  LENGTH  DETERMINA 
TION  *** 

20  DIM  IK500)  fL(500)  »R(350) 

30  ON  KEY*  1 »  *  ENTER  *  GOTO  170 
40  ON  KEY*  2  t  '  STORE  *  GOTO  600 
50  ON  KEY*  3  t 1  DISPLAY  *  GOTO  760 
60  ON  KEY*  At  ' EDIT  *  GOTO  940 
70  ON  KEY*  5  r ' CALC  *  GOTO  1200 
80  !  ON  KEY*  6»'HIST0'  GOTO 
90  !  ON  KEY*  7»'WEIBULL'  GOTO 
100  !  ON  KEY*  8  r ' NORMAL  *  GOTO 
110  ! 

120  CLEAR  0  KEY  LABEL 

130  DISP  *  INTERFACIAL  SHEAR  STRENGTH  DETERMINATION!' 

140  DISP  'CRITICAL  LENGTH  MEASUREMENT.' 

150  DISP  'SELECT  OPTION* 

160  GOTO  160 

170  !  ****  DATA  ENTRY  ***# 

180  CLEAR  0  DISP  'ENTER  DATA  FROM  KEYBOARD  (K)  OR  TAPE  <T)?' 

190  INPUT  Y $ 

200  IF  Yf  =  *  K  *  THEN  250 
210  IF  Y$=  *  T  *  THEN  500 
220  GOTO  180 
230  ! 

240  ! 

250  !  ***  ENTER  DATA  FROM  KEYBOARD  *** 

260  ! 

270  CLEAR  0  DISP  'ENTER  SAMPLE  NAME?' 

280  INPUT  A$ 

290  DISP  'MEAN  VALUE  OF  FIBER  CRITICAL  LENGTH  (MICRON)  ISt'»@  INPUT  W1 
320  DISP  'TENSILE  STRENGTH  <PSI)  OF  FIBER  AT'iWIJ'  GAGE  LENGTH  IS:'i@  IN 

PUT  S 

350  DISP  'THE  TOTAL  NUMBER  OF  SPECIMEN  TESTED*.  *?@  INPUT  N1 

380  DISP  'THE  TOTAL  NUMBER  OF  BREAKS 5  ' » 0  INPUT  N 

410  DISP  'ENTER  FIBER  DIAMETER  AND  CRITICAL  LENGTH  (MICRON)' 

420  FOR  1=1  TO  N 

430  DISP  'D('m')»  L(MIJ')-* 

450  INPUT  D(I)»L(I) 

460  NEXT  I 
470  GOTO  120 
480  ! 

490  !  ****  DATA  ENTRY  **** 

500  ! 

505  DISP  'INSERT  DATA  FILE  TAPE!  (WAIT  FOR  THE  NEXT  COMMAND!' 

506  BEEP 

507  WAIT  3000 

510  DISP  'ENTER  FILE  NAME. ' »0  INPUT  A* 

520  ASSIGN*  1  TO  UPC$(A*) 

530  READ*  1  i  A$»W1»S»N1»N?D( ) »L( ) 

570  ASSIGN*  1  TO  * 

580  GOTO  120 
590  ! 

600  !  ***  STORE  DATA  ON  TAPE  *** 

610  ! 

620  DISP  'INSERT  DATA  FILE  TAPE!  CWAIT  FOR  THE  NEXT  COMMAND ! 3 ' 

621  BEEP 

622  WAIT  3000 

630  CLEAR  0  DISP  'ENTER  FILE  NAME "  J  0  INPUT  A$ 
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640  DISP  'NEW  OR  EXISTING  FILE  <N/E>?' 

650  INPUT  Q$8  IF  UPC$(Q*C1 f13)='N'  THEN  660  ELSE  GOTO  670 
660  CREATE  UPC$(A*)f50 
670  ASSIGN#  1  TO  UPC$(A*> 

680  PRINT*  1  f  A# fWIfSfNIfNfDOfLO 
720  ASSIGN*  1  TO  * 

730  DISP  '  DATA  ON ■ J UPC* ( A$ ) f ' STORED  ON  TAPE*  8  KEY  LABEL 
735  WAIT  3000 
740  GOTO  120 
750  ! 

760  !  ***  DISPLAY  DATA  *** 

770  ! 

780  CLEAR  8  DISP  'DATA  IN  MEMORY?  (Y  OR  N)'f8  INPUT  T* 

790  IF  T$=‘ Y'  THEN  810 
800  IF  T$='N*  THEN  120 

805  DISP  'ENTER  Y  OR  N' 

806  GOTO  780 

810  CLEAR  8  PRINT  'SAMPLE' fA$ 

820  PRINT  USING  830  f  'TENSILE  STRENGTH  (PSI)  OF  FIBER  AT'fUIf'GAGE  LENG 
TH  ISJ'rS 

830  IMAGE  34Af4D.3D»16AfD.3DE 

840  PRINT  USING  850  f  'TOTAL  NUMBER  OF  SPECIMENS  TESTED  IS'.'fNI 
850  IMAGE  42Af2D 

860  PRINT  USING  870  i  'TOTAL  NUMBER  OF  BREAKS  ISJ'fN 
870  IMAGE  32Af3D 
880  FOR  1=1  TO  N 

890  PRINT  USING  900  f  'DCfIf'If  L < ' f I f ' ) = ' fD ( I ) f L < I > 

900  IMAGE  2Af3Df5Af3Df4Af2D.2Df5D.2D 

910  NEXT  I 
920  GOTO  120 
930  ! 

940  !  EDIT  ROUTINE  ***** 

950  ! 

960  CLEAR  8  DISP  'ADD(A)f  CHANGE(C)f  END(E)'f8  INPUT  Z* 

970  IF  Z$= ' A '  THEN  1000 
980  IF  Z$='C  THEN  1110 
985  IF  Z$=  *  E '  THEN  120 
990  GOTO  960 
1000  !  ADD  ROUTINE 

1010  DISP  'ENTER  ADDITIONAL  NUMBER  OF  SPECIMENS  TESTED*. 'f8  INPUT  Ml 

1020  DISP  'ENTER  ADDITIONAL  NUMBER  OF  BREAKS: *f8  INPUT  M 

1030  N1=N1 +M1 

1040  FOR  I  =  N  +  1  TO  N+M 

1050  DISP  *D(*fIf')f  L('fIf*)=' 

1070  INPUT  D(I)fL(I) 

1080  NEXT  I 

1090  N=N+M 

1100  GOTO  960 

1110  !  CHANGE  ROUTINE 

1112  DISP  'CHANGE  TENSILE  STRENGTH  (S)f  MEAN  CRITICAL  LENGTH  (W)f  DIAMET 
ER  AND  LENGTH  (L)' 

1113  INPUT  D$ 

1114  IF  D$= 'S'  THEN  1117 

1115  IF  D*= ' L '  THEN  1120 

1116  IF  B*=  * W '  THEN  1181 

1117  DISP  'TENSILE  STRENGTH='fS 

1118  DISP  'CHANGE  TO  TENSILE  STRENGTH31'  8  INPUT  S 

1119  GOTO  1180 
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1120  DISP  *  ENTER  ITEM  *  TO  CHANGE’ 

1130  INPUT  I 

1140  DISP  ’OLD  VALUES:* 

1150  DISP  aD<a>Ira)r  L  <  *  r  I  »  *  >  =  *  »  D  <  I  >  » L  <  I  > 

1160  DISP  ’NEW  VALUESf 

1170  DISP  *DC»Ifa)f  L<a*I»a>*M0  INPUT  D(I)»L(I) 

1180  GOTO  960 

1181  DISP  ’MEAN  CRITICAL  LENGTH=*,W1 

1182  DISP  *  CHANGE  TO  MEAN  CRITICAL  LENGTH=*  0  INPUT  W1 

1183  GOTO  960 
1190  ! 

1200  »  ***  WEIBULL  DISTRIBUTION  PARAMETERS  *** 

1210  ! 

1220  S1»S2fS3»CtF2=0 

1230  !  INITIALIZE  ALPHA  TO  BE  10 

1240  A= 1 0 

1250  !  INITIALIZE  ITERATION  STEP  TO  BE  1 
1260  X=1 

1270  «  ASPECT  RATIO  DETERMINATION 
1280  FOR  1=1  TO  N 
1290  R(I)=L(I)/D(I) 

1295  NEXT  I 

1299  CLEAR 

1300  DISP  'STORE  DATA  FOR  HISTOGRAM/WEIBULL  PLOT  ?  (Y  OR  N)*»0  INPUT  Q* 
1310  IF  Q$  =  * Y *  THEN  1330 

1320  IF  Q$=  *  N  *  THEN  1460 
1325  GOTO  1300 
1330  ! 

1340  !  STORE  ASPECT  RATIO  FOR  H I STO/ WE I BULL /NORMAL  PLOT  *** 

1350  ! 

1360  CLEAR  0  DISP  "ENTER  FILE  NAME?  (END  WITH  H)*»0  INPUT  B$ 

1370  DISP  ’NEW  OR  EXISTING  FILE  (N/E)?* 

1380  INPUT  P$0  IF  UPC$(P$Clrl3)=*N*  THEN  1390  ELSE  GOTO  1400 
1390  CREATE  UPC$(B$)»20 
1400  ASSIGN#  1  TO  UPC*(B$) 

1420  PRINT#  1  t  N »  R  <  > 

1440  ASSIGN*  1  TO  * 

1450  DISP  'DATA  ON  *»UPC*(B$)J*  STORED  ON  TAPE** 

1460  DISP  ’DO  YOU  WANT  A  PRINT  OUT?  (Y/N)?*  0  INPUT  C$ 

1461  IF  C$  =  *  Y  *  THEN  1463 

1462  IF  C$=  *  N  *  THEN  1468 

1463  PRINT  USING  1464  f  *N*»*R* 

1464  IMAGE  A*8X>A 

1465  FOR  1=1  TO  N  0  PRINT  USING  1466  t  I»R(I>  0  NEXT  I 

1466  IMAGE  3Dr5X,30.3D 

1468  BEEP 

1469  CLEAR 

1470  SI »S2  » S3=0 
1480  FOR  1=1  TO  N 

1490  S1=S1+R(I) ~A#LOG ( R ( I > ) 

1500  S2=S2  +  R  < I ) "A 
1510  S3=S3+L0G ( R ( I ) ) 

1520  NEXT  I 

1525  DISP  *  A=  *  » A 

1526  DISP  " F2=  * f F2 
1530  ! 

1540  F=S1/S2-1/ A-S3/N 
1545  DISP  *  F=  *  » F 
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1550  IF  ABS<FK. 000001  THEN  1680 
1560  IF  C>0  THEN  1600 
1570  F2=F 
1580  C=1 
1590  GOTO  1650 
1600  IF  F2>0  THEN  1630 
1610  IF  F<0  THEN  1650 
1620  GOTO  1640 
1630  IF  F>0  THEN  1650 
1640  X=-X/10 
1650  A=A-X 
1660  F2=F 
1670  GOTO  1470 
1680  ! 

1690  B=(1/N*S2)''(1/A) 

1700  P=l- 1/A 
1710  GOSUB  1960 
1900  G2=G 
1910  T=S/2/B*G2 

1920  P=l-2/A 

1921  DISP  •?=•?? 

1930  GOSUB  1960 

1940  U=S"2/4/B~2*<G-G2''2> 

1941  y=SQR(V) 

1945  GOTO  2310 
1950  ! 

1960  !  GAMMA  FUNCTION  SUBROUTINE  *** 

1970  ! 

1980  IF  P - 5 7 < = 0  THEN  2020 

2000  G=1 « E75 

2001  DISP  * P>57 » OVERFLOW t G  SET  TO  1.E75’? 

2002  WAIT  3000 

2003  BEEP 
2010  RETURN 

2020  P1=P 
2030  E= ♦ 000001 
2040  C=0 
2050  G=1 

2060  IF  P 1 -2<=0  THEN  2110 

2070  IF  P 1 -2<=0  THEN  2220 

2080  P1=P1-1 

2090  G=G*P1 

2100  GOTO  2070 

2110  IF  P1-1C0  THEN  2140 

2120  IF  Pl-1=0  THEN  2270 

2130  IF  Pl-1>0  THEN  2220 

2140  IF  PI -E>0  THEN  2190 

2150  Z=SGN<P1 )*INT(ABS(Pl)-.5)-Pl 

2160  IF  ABS ( Z ) -E<=0  THEN  2280 

2170  IF  1-Z-E<=0  THEN  2280 

2180  IF  Pl-1>0  THEN  2220 

2190  G=G/P1 

2200  P1=P1+1 

2210  GOTO  2180 

2220  Z=P1-1 

2230  W=Z*(-.5684729+Z*( .2548205+Z*-. 0514993) ) 
2240  W=Z*( . 985854+Z* <-. 87642 18+Z*( ,8328212+U) ) ) 
2250  Gl=l+Z*<-»5771017+W> 


2260  G=G*G1 


2270  RETURN 

2280  DISP  "P  IS  WITHIN  l.E-6  OF  BEING  A  -VE  INTEGER' 

2281  WAIT  3000 

2282  BEEP 
2290  RETURN 
2300  ! 


2310  !  DATA  OUTPUT  *** 


2320  PRINT 
2330  IMAGE 
2340  PRINT 
2350  IMAGE 
2360  PRINT 
2370  IMAGE 
2380  PRINT 
2390  IMAGE 


USING  2330  f 
7A » 10A 
USING  2350  i 
26A » 3D 

USING  2370  i 
23A » 4D 

USING  2390  i 
27A » 3D « 3D 


•SAMPLE? 'fA$ 

•TOTAL  NUMBER  OF  SPECIMENS* * * N1 
•TOTAL  NUMBER  OF  BREAKS='»N 

• WEIBULL  DISTRIBUTION  ALPHA**  *  A 


2400  PRINT 
2410  IMAGE 
2420  PRINT 


USING  2410 
26A  1 4D « 3D 
USING  2430 


•WEIBULL  DISTRIBUTION  BETA**  »B 
•TENSILE  STRENGTH  (PSI)  OF  FIBER  AT*»W1» 


ENGTH=*  t S 


2430  IMAGE 
2440  PRINT 
2450  IMAGE 
2460  PRINT 
2470  IMAGE 


35A»4D.3D»13A>D.3DE 


USING  2450  i 
27Af D*3DE 
USING  2470  t 
31A»D«3DE 


•INTERFACIAL 
•VARIANCE  OF 


SHEAR  STRENGTH*  *  t T 
THE  SHEAR  STRENGTH**  » V 


2480  • 
2490  END 


GAGE  L 
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APPENDIX  I.  WEI BULL  DISTRIBUTION  OF  CRITICAL  LENGTH 
A  HISTOGRAM  PLOT 


The  purpose  of  this  program  is  to  plot  the  Weibull  Distribu¬ 
tion  curve  and  the  population  density  curve  of  the  critical  aspect 
ratio  of  the  fibers. 
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10  COM  Z9 

20  !  *****DATA  PLOT  PROGRAM*#*#*******VERSION  6.09********* 

30  DIM  P<2) *L<2>  *N1(2> *S2<2> »X< 2*350) f T*C80] » X$C803 f Y$C803 » S*C80  3 » LI ( 2 ) * 
Y1(2*350)*R(350) 

40  DIM  T1$C80] 

50  !  *****INITIALIZATION************ *********************** 

60  G*= '  AG  *  0  N$='NY'  0  K*='KTC“  0  L1(0)=0  9  S$=“ “  9  N  =  0 
70  N1=0  6  Z9=l  0  DEG 
80  ON  KEY#  1 * 'ENTER*  GOTO  180 
90  ON  KEY*  2r 'SETUP'  GOTO  530 
100  ON  KEY*  3 * *  HEADING  *  GOTO  830 
110  ON  KEY*  4 »  '  PLOT  *  GOTO  2130 
120  ON  KEY*  5  *  *  STORE  *  GOTO  390 
130  ON  KEY#  6  * *  EDIT  *  GOTO  1460 
140  ON  KEY#  7  *  '  LEGEND*  GOTO  1390 
150  ON  KEY#  8*'  RANGE*  GOTO  900 
160  CLEAR  0  KEY  LABEL 
170  GOTO  170 
180  ! 

190  !  *****DAT A  ENTRY************************************* 

200  CLEAR  0  DISP  'Enter  Date  from  Keyboard ( K)  or  Tape ( T)  or  from  Calcula 
t  i  o  n  ( C )  * 

210  INPUT  Q$0  Q2=P0S(K$»UPC*(Q$Cl»n) ) 

220  IF  Q2=0  THEN  200 
230  ON  02  GOTO  240*750*741 
240  FOR  I =N  + 1  TO  6 

250  CLEAR  0  DISP  'Enter  Data  for  Set  '515'  <Y/N)' 

260  INPUT  Q$@  Q2=P0S<N$fUPC$(Q$Clf 1]>) 

270  IF  02=0  THEN  250 
280  ON  02  GOTO  370*290 

290  DISP  “Enter  Data  Pai rs ( ENDLINE  ends)' 

300  FOR  J=1  TO  100 

310  DISP  'XCflJ'f'JJ*')*  YC5I*  '»'*  Jf ')  =  '*0  INPUT  Q$*R* 

320  IF  Q$  = "  THEN  350 

330  X(I* J)=VAL(Q$)  0  Y1 < I  * J)=VAL(R$> 

340  NEXT  J 

350  N1 ( I >  = J-l  0  N1=N1+N1(I) 

360  NEXT  I 
370  N=I-1 

380  ON  Z9  GOTO  530*160 

390  !  *****STORE  DATA  ON  TAPE*** 

400  Z9=2 

410  CLEAR  0  DISP  'Store  on  T ape ( Y/N > ■ 

420  GOSUB  3190  0  ON  09  GOTO  160*430 
430  DISP  “FILE  NAME' 50  INPUT  F$ 

440  DISP  “New  or  Existing  File  (N/E)' 

450  INPUT  Q$@  IF  UPC$ ( Q$C 1  *  13 )  =  • N ■  THEN  460  ELSE  GOTO  470 
460  CREATE  UPC$(F#)*50 
470  ASSIGN*  1  TO  UPC$<F*> 

480  PRINT*  1  »  N»N1< ) >X0»Y0»X1 *Y1 »X( * ) »Y1 <  * ) »L1< ) »P( ) »L( ) »S( ) »S$»T$»Y$*X 

$*F1*F2*Q1»X2*X3*Y2*Y3 

490  PRINT*  1  *  Tl$ 

500  ASSIGN*  1  TO  * 

505  CLEAR 

510  DISP  'DATA  SET  '5UPC$(F*>*'  STORED  ON  TAPE'  0  KEY  LABEL 
520  Z9=2  0  GOTO  170 

530  !  *****DETERMINE  COLOR*  LINE  AND  SYMBOL************** 

540  CLEAR 
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550  FOR  1=1  TO  N 


560  DISP 

570  DISP  *  ENTER  PEN  #  (1-4),*  8  DISP  'LINETYPE  *  (1-8)  AND  PL 

OTTING  SYMBOL  ♦  (1-6)" 


580  DISP  'FOR  DATA  SET*»I 
590  IF  1=1  THEN  600  ELSE  700 

600  PRINT  'PEN  COLOR  KEY*  8  PRINT  *1=  BLACK*  8  PRINT  *2  RED*  8  PRINT  *3 
GREEN*  8  PRINT  *4  BLUE* 

610  PRINT  8  PRINT  'LINETYPE  KEY*  8  PRINT  *1=  SOLID* 

620  PRINT  *2  BLANK*  e  PRINT  *3  DOT*  8  PRINT  *4  S*DASH* 

630  PRINT  *5  L.DASH*  8  PRINT  *6  DASH/DOT* 

640  PRINT  *7  L ♦ DASH/S ♦ DASH  *  8  PRINT  *8  DASH/DOT/DOT* 

650  PRINT  8  PRINT 

660  PRINT  'PLOTTING  SYMBOL  KEY* 

670  PRINT  *1=  DOT*  8  PRINT  *2  CIRCLE*  8  PRINT  *3  SQUARE* 

680  PRINT  *4  TRIANGLE*  0  PRINT  *5  PLUS  SIGN*  0  PRINT  *6  DIAMOND* 

690  PRINT  0  PRINT  8  PRINT  0  PRINT 
700  INPUT  P(Z)rL(I)rS(I) 

710  IF  P<IX1  OR  P  ( I )  >  4  THEN  700 

720  IF  L ( I ) < 1  OR  L ( I ) >8  THEN  700 

730  IF  S(IK1  OR  S  ( I )  >  6  THEN  700 

740  NEXT  I  0  ON  Z9  GOTO  830*160 

741  !  CALCULATION  FOR  WEIBULL  DISTRIBUTION 

742  DISP  ‘ENTER  FILE  NAME  WITH  RO  MATRIX  FOR  WEIBULL  DISTRIBUTION*  0  IN 
PUT  F $0  ASSIGN*  1  TO  F$ 


743  READ*  1  *  N5*R(  ) 

744  ASSIGN*  1  TO  * 

745  N=1  0  N 1 < 1 > =N5 

746  DISP  *  ENTER  VALUES*  0  INPUT  A  *B 

747  CLEAR  0  DISP  "CUMULATIVE  PLOT  (C)  OR  DENSITY  POPULATION  PLOT 
INPUT  E* 


( D )  ?  * 
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748  IF  E*=  * C *  THEN  750 

749  IF  E*=*D*  THEN  752  ELSE  747 

750  FOR  1  =  1  TO  N5  8  X< 1 » I >  =  <R< I >/B)"A  8  Y1 ( 1 » I ) =1-EXP ( -X ( 1  *  I ) )  e  X(1*I)  = 
R  < I )  0  NEXT  I 

751  Z9=2  8  GOTO  160 

752  FOR  1=1  TO  N5  e  X( 1 t I >=(R< I )/B)"A  8  Y1 ( 1 f I ) =A/B*(R< I >/B)~ ( A-l ) *EXP( - 
X ( 1  *  I ) )  8  X( 1  *  I )=R< I ) 

753  NEXT  I 

754  Z9=2  0  GOTO  160 


759  !  DATA  ENTRY*#**#**####*### 

760  DISP  'ENTER  FILE  NAME  *  *  0  INPUT  F* 

770  ASSIGN*  1  TO  UPC$(F$) 

780  READ*  1  *  N »N1 ( ) *X0 * Y0»X1 » Y1 »X ( » ) * Y1 ( t > »L1 < > »P( ) »L< ) »S( ) tS% * T$ » Y$ »X* 

*Fl*F2rQl»X2*X3*Y2»Y3 

790  READ*  1  J  Tl* 

800  ASSIGN*  1  TO  # 

810  Z9=2 
820  GOTO  160 


830  '  *****ENTER  HEADINGS****************************** 
840  CLEAR  0  DISP  'ENTER  TITLE  *  >  0  INPUT  T* 

850  DISP  'ENTER  SUBTITLE  *  » 0  INPUT  Tl$ 

860  DISP  'ENTER  X  axis  TITLE* ?0  INPUT  X* 

870  DISP  ’ENTER  Y  axis  TITLE* »8  INPUT  Y* 

880  P(0)=1 

890  ON  Z9  GOTO  900*160 

900  i  #####DEF I NE  GRAPH******************************** 

910  XO » Y0=10"99  0  XI » Yl=-10”99 
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1010  FOR  1=1  TO  N 
1020  FOR  J  =  1  TO  N1 ( I ) 

1030  XO=MIN(XO?X(I?J))  (?  X1=MAX(X1?X(I?J)> 

1040  Y0=MIN(Y0?Y1(I?J) )  @  Y1=MAX ( Y1 ? Y1 ( I ? J > > 

1050  NEXT  J 
1060  NEXT  I 

1070  DISP  "X  DATA  RANGE  IS  "JXO?"  TO  "?X1 

1080  DISP  "  Y  DATA  RANGE  IS  "?Y0?"  TO  "JY1 

1100  DISP  *  DO  YOU  WISH  TO  SET  GRAPH  LIMITS  "?0  GOSUB  3170 

1110  ON  Q9+1  GOTO  1 100 ? 1 160 ? 1 120 

1120  DISP  "ENTER  Xmin » Xmax  ?  Yinin  ?  Ymax  "? 

1140  INPUT  XO? XI ? YO?  Y1 

1150  IF  X1<X0  OR  Y1CY0  THEN  1120 

1160  DISP  "DO  YOU  WISH  AXES  OR  GRID  (A  or  G)"? 

1170  INPUT  Q$0  91=P0S(G*?Q$) 

1175  DISP  @  DISP  "NUMERICAL  LABELS: " 

1177  DISP  "on  horizontal  axis..." 

1180  ON  Gl+1  GOTO  1 160 ? 1 190 ? 1240 

1190  DISP  "Units  per  label?  units  per  tic  mark"? 

1200  INPUT  X3  ?  X20  X3  =  INT ( X3/X2+ . 51 )  0  IF  X3<0  THEN  1190 
1205  DISP  0  DISP  "on  vertical  axis..." 

1210  DISP  "Units  per  label?  units  per  tic  mark"? 

1220  INPUT  Y3 ? Y20  Y3  =  INT ( Y3/Y2 )  0  IF  Y3<0  THEN  1210 
1230  GOTO  1290 

1240  DISP  "units  per  GRID  LINE?units  Per  tic  mark"? 

1250  INPUT  X3  r X20  X3  =  INT < X3/X2+ ♦ 5 1 )  0  IF  X3<0  THEN  1240 
1255  DISP  0  DISP  "on  vertical  axis..." 

1260  DISP  "units  per  GRID  LINE?units  per  GRID  LINE"? 

1270  INPUT  Y3  ?  Y20  Y3=INT ( Y3/Y2+ .51 >  0  IF  Y3<0  THEN  1260 

1280  ON  ERROR  GOSUB  3220 

1290  IF  Y2<1  THEN  F2=2  ELSE  F2=0 

1300  IF  Y 2 < » 1  THEN  F2=4 

1310  IF  X2<1  THEN  Fl=2  ELSE  F1=0 

1320  IF  X2C.1  THEN  Fl=4 

1330  DISP  "INTEGER  LABELS  ON  X-AXIS"? 

1340  GOSUB  3170 

1350  ON  09  GOTO  1370?1360 

1360  F1=0 

1370  OFF  ERROR 

1380  ON  Z9  GOTO  1390»160 

1390  !  #**#*#*ENTER  LEGENDS****** 

1400  S$="  0  CLEAR 
1410  FOR  1=1  TO  N 

1420  DISP  "Enter  Legend  for  Group  *?I 
1430  INPUT  I$@  S*=S$&I$  0  LI ( I ) =LEN ( S$ ) 

1440  NEXT  I 

1450  ON  Z9  GOTO  390? 160 
1460  !  EDIT  ROUTINE 

1470  CLEAR  0  DISP  "Enter  data  set  number" 

1480  INPUT  I 

1490  DISP  "There  are"?Nl(I)?"  data  pairs  (X?Y)  in  this  set" 

1500  DISP  0  DISP  "Enter  starting  &  ending  item  numbers  for  list  (ENDLINE 

ends) " 

1510  INPUT  D$?E$ 

1520  IF  D*= " "  THEN  2122 

1530  FOR  J=VAL(D$)  TO  MIN (UAL (E* >  ?N1 ( I > ) 

1540  DISP  J?"="?X<I?J)?"?"?Y1(I?J) 

1550  NEXT  J 
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1560  DISP  *  Add ( A)  f  Delete(D)  FChan£Je(C>  f  Insert  ( I )  »End(E>  * 
1570  INPUT  Yl$ 

1580  IF  UPC*(Y1*C1  f13)=*A*  THEN  1640 
1590  IF  UPC$(Y1$C1»13)=*D*  THEN  1700 
1600  IF  UPC$(Y1$C1f13)='C*  THEN  2050 
1610  IF  UPC$( Y1$C1 f 13 )  =  * I  *  THEN  1850 
1620  IF  UPC$(Y1$C1f13)='E*  THEN  1500 
1630  GOTO  1500 
1640  !  ADD  ROUTINE 

1650  DISP  •XCJIJ'f  * JN1<I)+1J*>fY<*  JIJ * f * J N1 < I > +1 J * >  =' 
1660  N1(I)=N1(I)+1 

1670  INPUT  X(IfN1(I))fY1(IfN1<I>) 

1680  GOTO  1560 
1690  ‘ 

1700  !  DELETE  ROUTINE 

1710  DISP  'Enter  item  #  to  delete* 


1720  INPUT  J 

1730  DISP  *X( ' fit  * r * i Jf * ) »Y< * fI *  * f * i Jf * > 
1740  DISP  'Enter  D  to  delete* 


*JX(I»J)»*  *»Y1(I»J) 


1750  ! 

1760  INPUT  Yl$ 

1770  IF  UPC*<Y1$C1»13)*'D*  THEN  1560 
1780  N1(I)=N1(I>-1 
1790  FOR  J1=J  TO  N1 ( I ) 

1800  X<I» J1)=X<I» Jl+1) 

1810  ! 

1820  Y1(I» J1>=Y1(I» Jl+1) 

1830  NEXT  Jl 

1840  DISP  *Item*JJ»*  deleted*  Q  GOTO  1560 
1850  !  INSERT  ROUTINE 

1860  DISP  'Insert  item  before  item  ** 

1870  ! 

1880  INPUT  J 

1890  DISP  ■X<*m*»B>J»,>#Y<*»I»'»MJ**>  =  * 

1900  INPUT  G8fG9 

1910  N1 ( I ) =N1 ( I ) +1 

1920  FOR  J1=N1 ( I )  TO  J+l  STEP  -1 


1930  ! 

1940  XCI» Jl)*X(If Jl-1) 

1950  Yltlf Jl)*Yl(If Jl-1) 

1960  NEXT  Jl 
1970  X  < I f  J ) =68 
1980  Y1(IfJ)=G9 
1990  ! 

2000  GOTO  1560 
2010  ! 

2020  ! 

2030  • 

2040  ! 

2050  \  CHANGE  ROUTINE 

2060  DISP  'Enter  Item  #  to  ChansSe' 

2070  INPUT  J 

2080  DISP  'Old  Values:*  0  DISP 
2090  DISP  X(IfJ)f*f*fY1(IfJ) 

2100  DISP  0  DISP  'Enter  New  Values* 

2110  INPUT  X(IfJ)fY1(IfJ)0  DISP  0  DISP 
2120  GOTO  1560 

2122  CLEAR  0  PRINT  USING  *9/*  f  0  PRINT 


■Press  RANGE  User  Key  before 

1-5 


PLOT"  0  KEY  LABEL  0  GOTO  170 

2130  !  *#***GRAPH****************************************** 

2140  BISP  'PLOT  ON  ♦  * 

2150  DISP  ■  1.  HP  SCREEN* 

2160  DISP  *  2*  PLOTTER*  0  DISP  *  3*  POSTER* 

2170  INPUT  Z 

2180  IF  Z=1  THEN  2270 

2190  IF  Z=2  THEN  2230 

2200  IF  Z=3  THEN  2230 

2210  GOTO  2140 

2220  Z9=2 

2230  PLOTTER  IS  705 
2240  IF  Z  =  3  THEN  2280 
2250  LIMIT  0*  250* 15*195 
2260  GOTO  2280 
2270  PLOTTER  IS  1 
2280  LOCATE  26*120*20*80 
2290  SCALE  X0*X1*Y0*Y1 

2300  FXD  FI  *F2  0  CSIZE  4m6»0  0  LINETYPE  1  0  GCLEAR  0  PEN  1 
2310  ON  Q1  GOTO  2320>2340 
2320  LAXES  -X2*Y2*X0*Y0*X3*Y3 
2330  GOTO  2350 

2340  LGRID  -X2 * Y2 * XO * YO * X3 * Y3  0  FRAME 

2350  PEN  P(0)  0  MOVE  X0+ ( Xl-XO > /2 * Y1+ ( Yl-YO ) *3/20 

2360  LORG  4  0  CSIZE  2#Z+2*.7*0 

2370  LABEL  T$ 

2380  CSIZE  2#Z+1*.7 
2390  LABEL  Tl$ 

2400  SETGU 
2410  MOVE  70*3 
2420  LORG  4 
2430  LABEL  X$ 

2440  MOVE  8*45 

2450  LDIR  90  0  LABEL  Y$ 

2460  LDIR  0 
2470  SETUU 

2650  !  #*#**PLOT  DATA**#*#**#####**#*#*###*#****#****#*#* 

2660  FOR  1=1  TO  N 
2670  PEN  P(I) 

2680  MOVE  X(I*1) *Y1(I»1) 

2690  GOSUB  2890 

2700  MOVE  X(I*1) *Y1(I*1) 

2710  FOR  J=2  TO  NKI) 

2720  LINETYPE  L<I> 

2730  DRAW  X ( I  * J ) * Y1 < I > J ) 

2740  GOSUB  2890 

2750  MOVE  X(I»J)»Y1(I»J) 

2760  NEXT  J 
2770  GOSUB  2800 
2780  NEXT  I  0  PEN  0 

2782  IF  Z=1  THEN  CLEAR  0  DISP  USING  *7/*  »  0  DISP  'Press  KEY  LABEL  when 
ready*  0  WAIT  2000 

2784  IF  Z=1  THEN  CLEAR  0  GRAPH  0  GOTO  170 
2790'  GOTO  160 

2800  !  *****DRAW  KEY************************************** 

2810  SETGU  0  LINETYPE  L(I) 

2820  MOVE  110*96-3# ( 1-1 ) 

2830  IDRAW  5*0 
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2840  GOSUB  2890  0  SETGU 
2850  MOVE  117 f  96-3# ( 1-1 ) 

2860  LORG  2  0  CSIZE  3f»6  0  LDIR  0 
2870  LABEL  S*CL1 < 1-1 > +1 fLI ( I > 3 
2880  SETUU  0  RETURN 

2890  !  *****PLOT  SYMBOLS********************************** 
2900  SETGU  0  LINETYPE  1 

2910  ON  3(1)  GOTO  2920f2940f2990f3040f3080f3130 
2920  !  *###*DOT####*##*###***### 

2930  RPLOT  OfO  0  SETUU  0  RETURN 
2940  !  *##*#CIRCLE#**#*****#**** 

2950  IMOVE  ♦  15  f  0 

2960  FOR  A=0  TO  360  STEP  30 

2970  PDIR  A  0  RPLOT  .8f0 

2980  NEXT  A  0  SETUU  0  RETURN 

2990  !  ####*SQUARE*#*#####*##*** 

3000  IMOVE  - 1 f - 1 

3010  IDRAW  2  f  0  0  IDRAW  0f2 

3020  IDRAU  -2 f 0  0  IDRAW  Of-2 

3030  SETUU  0  RETURN 

3040  !  *##**TRIANGLE ************ 

3050  IMOVE  OfI 

3060  IDRAW  1 f -2  0  IDRAW  -2f0 
3070  IDRAW  1 f 2  0  SETUU  0  RETURN 
3080  !  *****PLUS  SIGN*********** 

3090  IDRAW  OfI  0  IDRAW  Of-2 

3100  IMOVE  OfI 

3110  IDRAW  IfO  0  IDRAW  -2f0 

3120  SETUU  0  RETURN 

3130  !  ###**D I AMOND# ************ 

3140  IMOVE  OfI  0  IDRAW  1f-1 
3150  IDRAW  -1 f -1  0  IDRAU  -1f1 
3160  IDRAW  1 f 1  0  SETUU  0  RETURN 

3170  !  *****YES  OR  NO************************************** 
3180  DISP  * (Y/N) * 

3190  INPUT  0*0  Q9=P0S ( N* f  0* ) 

3200  IF  Q9=0  THEN  3180 
3210  RETURN 

3220  !  *****CHECK  FOR  DIVISION  BY************ZERO******* 
3230  IF  ERRN=8  THEN  FI fF2=0 
3240  RETURN 


APPENDIX  J.  INSTRUMENTED  IMPACT  TESTER 


The  ETI-Dynatup  drop  dart  instrumented  impact  tester  monitors 
and  records  the  entire  impact  event,  from  initial  impact  and 
acceleration  from  rest  to  plastic  bending  to  fracture  initiation 
and  propagation  to  failure.  Also,  the  drop  dart  system  permits 
dynamic  tests  to  be  made  at  impact  energies  and  velocities  not 
available  with  other  testing  machines. 

The  drop  dart  impact  test  consists  of  a  drop  tower  and  instru¬ 
mentation  package  -  a  microprocessor.  The  drop  tower  is  a  gravity 
driven  impact  machine  equipped  with  remote  controls  for  release 
of  the  hammer  and  tup  assembly,  and  a  motorized  lift  mechanism 
for  easy  return  of  the  hammer  to  a  predetermined  drop  position. 

An  automatic  rebound  brake  is  furnished  as  standard. 

The  hammer  and  tup  assembly  is  designed  so  that  the  weight  of 
the  hammer  can  be  varied  and  the  tups  can  be  interchanged.  The 
base  plate  of  the  drop  tower  permits  interchangeability  of  anvil 
supports  to  accommodate  different  specimen  size  and  geometries. 

A  typical  output  is  shown  in  Figure  J-l.  Pj  and  Pp  are  the  loads 
correspond  to  the  initial  damage  of  the  composite  and  the  maximum 
load  required  for  punch  through.  The  corresponding  energy  ab¬ 
sorbed  at  each  stage  is  Ej  and  EF.  The  total  energy  absorbed 
during  the  whole  event  is  Erp.  In  this  report,  Pp  and  Eip  are  used 
to  characterize  the  impact  toughness  of  the  laminates. 
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